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1.0  INTRODUCTION 

B-K  Dynamics'  activities  during  the  third  quarter  (15  April  1975  - 
15  July  1975)  continued  to  focus  on  preparation  of  hardware  for  the 
interim  STINGER  simulation  and  conversion  of  the  STINGER  simula- 
tion to  the  ASC's  hybrid  computer  system.  This  report  describes  the 
STINGER  conversion  activities  and  presents  program  flow  charts  and 
coding  generated  in  this  task.  ij 

'I 
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2.  0 STING i:U  CONVERSION  TO  ASC  EQUIPMENT 


Converting  the  STINGER  real-time  simulation  from  the  interim 
system  to  the  ASC  equipment  has  required  replacing  IBM- 7094  soft- 
ware functions  with  equivalent  CDC-6600  functions,  and  IBM-UOS 
interface  operations  with  equivalent  DADIOS  ADC,  DAC  and  discrete 
word  liandling  capabilities.  In  Figure  2.  1 the  new  ASC  equipment 
utilized  in  phase  1 conversion  is  shown  in  relation  to  the  overall 
STINGEF^  simulation.  The  items  which  represent  new  ASC  equip- 
ment are  enclosed  in  dotted  lines. 

The  documentation  given  herein,  of  the  STINGER  simulation, 
reflects  the  rationale  and  realization  of  an  attempt  to  duplicate 
interim  STINGER  functions  on  ASC  equipment.  This  conversion 
took  place  within  the  constraints  of  available  documentation,  insta- 
' lation  of  new  hardware,  availability  of  new  hardware,  and  an  attempt 

i to  minimize  the  effect  of  phase  1 conversion  on  other  modules  of  the 

[ simulation.  The  phase  1 STINGER  conversion  presented  in  the 

! following  sections  represents  a one-on-one  conversion  of  software 

! and  hardware. 
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FIGURE  2.  1 PHASE  1 CONVERSION  OF  STINGER  SIMULATION 


2.  1 SOFTWARE  CONVERSION  OF  REAL-TIME  CODE 


i.f 

r 


Previously,  the  portion  of  the  interim  STINGER  simulation  im- 
plemented on  the  IBM/ 7094  consisted  of  approximately  2000  lines 
of  code  in  FORTRAN  and  500  lines  of  code  in  MAP  (IBM-7094 

assembly  language).  Converting  this  software  to  the  CDC-6600 

>!< 

required  minor  changes  to  the  FORTRAN  program  and  completely 
new  code  written  in  FORTRAN  hybrid  to  replace  the  IBM-7094 
assembly  code. 

This  new  code  performs  the  identical  functions  of  the  interim 
simulation.  These  functions  include  pre- real-time  data  processing, 
real-time  computations  and  post-real-time  data  analysis.  A 
summary  of  the  functions  performed  by  each  of  these  three  major 
program  segments  is  summarized  in  Appendix  G. 

In  Figure  2.2,  a flow  chart  of  the  real-time  portion  of  the  phase  1 
STINGER  program  is  presented.  An  explanation  of  the  symbols  used 
in  Figure  2.2  are  given  in  Table  2. 1.  Furthermore,  a limited  description 
of  the  computational  equations  are  given  in  Appendix  A.  These  equations 
represent  a one-on-one  cross  reference  between  program  symbols  and 
mathematical  notation. 

The  real-time  software  presented  here  has  been  tested  in  real-time 
using  simplified  test  functions.  The  simplified  test  functions  used  are 
presented  as  part  of  a FORTRAN  hybrid  real-time  subroutine  in  Appendix 

H.  These  simplified  test  functions  have  been  used  in  the  absence  of 
real  test  data  from  the  existing  simulation. 


*The  modifications  to  the  FORTRAN  program  were  performed  by 
MICOM  personnel. 
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FIGURE  2.  2 FLOW  CHART  OF  STINGER  REAL  TIME  CODE  ' i 
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PAUSE  2 


LOW 


iirrOLD 


PRINT  REAL 
TIME  STATU 


RETURN 


NPATH 


Pause  for  static  check 
ok  signal  then  give  DDS 
a GO  signal 


(Transfer  initial 
conditions 


Signal  ICS 
sent 


Enter  real  time  with 
GO  from  DDS 


READ 


Wait  for  ramp  up 
sijfnal 


If  ITl-NT  = 0 
table  has  been 
exceeded 


Check  to  be 
sure  ICR2>1 


Begin  interpolation 
of  XCOMP,  YCOMP, 
ZCOMP.YCandZC 


IT2  = 
ITl  + 1 


WMAN 


IND  = IND  + 1 


CALC 


MAX-IND. 


MAN(4*IND-3)=DT 


CALC 


Store  an  ADC 
if  range  table 
entry  is  reached 


DIV=TS(IT2)-TS{IT1) 

RATIO=(DT-TS(ITl))/DIV 

XCOMP=XDTGMS(ITl)+RATIO*(XDTGMS(IT2)-XDTGMS(ITl)) 

XC=XDM(ITl)+RATIO*(XDM(IT2)-XDM(ITl)) 

YCOMP=YDTGMS(ITl)+RATIO*(YDTGMS(IT2)-YDTGMS(ITl)) 

YC=YDM(ITl)+RATIO*(YDM(IT2)-YDM(ITl)) 

ZCOMP=ZDTGMS(ITl)+RATIO’i'(ZDTGMS(IT2)-ZDTGMS(ITl)) 

ZC=ZDM(IT1)+RATI0’!'(ZDM(IT2)-ZDM{IT1)) 


INTERP(TREAL.  TMA,  ZE,  NT,  NPX,  ZTA,  NERR) 


A1  = XTA^  + YTA^  + ZTA^ 


Begin  calcu- 
[ lation  of  COSE 
’rLB  and  A1 


= 0,  01745329  * 4.  637084242  ■:=  Al/(PHO  - VTI  * CLA) 


SA  = SIN(A1) 

CA=COS(Al) 

XXX=XXX=;‘20475.  O/SKK 
YYY=YYY=:=20475.  O/SKK 
ZZZ  = ZZZ=:--20475.  O/SKK 

RRR=XXX^+YYY^+ZZZ^ 

RFEET=RRR 

RRR=VTp:=-y/RRR 

F1  = S2=:=XXX-YYY*SPL+S3=:=ZZZ 

F2=XXX*S4+YYY*CPL+ZZZ*S5 

F3=CTL*ZZZ-STL-XXX 

SR=>/xC^+YC^ 

S1=CA+(SA/SR)=^ZC 

G1  = S1«XC 

G2  = S1*YC 

G3=ZC=:=CA-SA*SR 

E111=(F1*G1+F2*G2+F3*F3)/RRR 

COSE=l-Elll^ 

RLB=-y/COSE 

RLB=RLB*RLBK/SCALET 
COSE=Elll/l.  02375 


Begin  computation 
of  iris  ratio  no.  7 


NX  = 17 
XPX  = 2 

CALL  IXTERP(E111,  FCOS.  IIR7EF,  NX,  NPX,  HR7E.  NERR) 
HC7FT=1209.  675^:RIR7E-EXP(-2.  341*AL0G(RFEET/22965.  831)) 
P1IRSS=(2.  0 •=ATAN(0,  04884004884/RLB)) 

P2IRSS=230.  34375  ' RLB/RFEET 

AT7  = ((P2IRSS=:==:=2)-SlN(PlIRSS))/(2.  0*(COS(PlIRSS/2.  0))-'‘*2) 
PJTU1  = 0.  I’i'EXPd.  003258-ALC)G(AT7)) 

FI7  = IIC7FT/PJTU1 
XX  = 9 
NPX =2 

CALL  IXTERP(FI7.  FI7T.  FCI77,  NX.  NPX,  FCI7.  NERR) 

RN=FCI7 

Note:  for  MICOM  hybrid  applications  RN  is  a 
random  number.  See  Appendix  C for 
appropriate  code. 


RC1=Vf1^+F2^ 

RC1=F2/RC1 

RCB=ARCC0S(RC1) 

RCX=-RC1 

RCY  = -SIN(RCB) 


Begin  compu- 
tation of  scaled 
plume  rotation 
angle 


RCY= -RCY 
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I Output  DAC 


Transfer  to 
entry  point 
READ 


TABLE  2.  1 

DESCRIPTION  OF  PROGRAM  SYMBOLS 


VARIABLE  I 


PROGRAM 

SYAIBOL* 


TRP 


THETAL 


DESCPJPTION 


Scaled  plume  rotation  angle 

Plume  rotation  angle  (RAD.  ) 

Initial  elevation  angle  of  target 
(DEG. ) 

Real  time  (SEC.  ) 

Gravitational  constant 

Target  aerodynamic  accelera- 
tion table  in  earth  fixed  coor- 
dinates 


Target  angle  of  attack  (RAD.  ) 

Time  at  approximately  1000 
feet  to  go 

Aerodynamic  lift  coefficient 
due  to  angle  of  attack 

Target  inertial  velocity 

Apparent  plume  length  to 
breadth  ratio 


Subscript  which  denotes  earth 
fixed  coordinates 


I 


VARIABLE 

PROGRAM 

SYMBOL 

DESCRIPTION 

GN 

Subscript  which  denotes  gener- 
alized target  coordinates 

F 

Subscript  which  denotes  target 
fixed  coordinates 

G 

Subscript  which  denotes  guid- 
ance coordinates 

L 

Subscript  which  denotes  launch 
coordinates 

Rj.  R 

RI.  RFEET 

Initial  range  (FT, ) 

EDOT 

Scaled  7/t' 

GAM 

P 

PHO 

3 

Air  density  (Slugs/FT,  ) 

"''min 

GGG 

Minima  of  the  overload  function 

Scalet 

SCALET 

Scale  Factor  = 10.  2375 

K 

SKK 

Variable  scale  factor. 

K= 

^ 1+7  , t > t' 

XX 

Range  table  for  data  collection 

1 

XXS 

Scaled  range  table  (XX) 
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VARIABLE 

PROGRAM 

SYMBOL 

DESCRIPTION 

^r7 

RN 

For  MICOM  Hybrid,  a uniform 
random  number,  RN  f (-1,1) 

For  IRSS,  iris  ratio  number  7 

Cos  ( e ) 

COSE 

Cosine  of  angle  between  LOS 
and  center  line  of  target 

DX 

/ 

DY 

/ Scaled  xxx,  yyy,  zzz 

DZ 

’ 

t 

DT 

Real  time 

XXX 

x-Missile  minus  x~target 
position  (FT. ) 

^ml'^tl 

YYY 

y- Missile  minus  y- target 
position  (FT. ) 

ZZZ 

z- Missile  minus  z- target 
position  (FT. ) 

XPOT 

x-Missile  minus  x- target 
velocity  (FT./ SEC. ) 

« • 

1 

YDOT 

y- Missile  minus  y- target 
velocity  (FT.  / SEC. ) 

^ml'^tl 

ZDOT 

z-Missile  minus  z-target 
velocity  (FT.  / SEC.  ) 

1 I 


VA  RUBLE 

PROGRAM 

SYMBOL 

DESCRIPTION 

^TL 

XDTGMS 

Target  manuvers  and  manuver 
time 

^TL 

YDTGMS 

1 

^TL 

ZDTGMS 

e 

Apparent  plume  length  (FT. ) 

B 

Apparent  plume  breadth  (FT.  ) 

Cos  ^(L  / 7L2+L2) 

h 

ZALT 

Altitude  above  sea  level 

XCOMP 

Interpolated  value  of  XDTGMS 

YCOMP 

Interpolated  value  of  YDTGMS 

ZCOMP 

Interpolated  value  of  ZDTGMS 

^TE 

XC 

Interpolated  value  of  XDM 

Yte 

YC 

Interpolated  value  of  YDM 

^TE 

ZC 

Interpolated  value  of  ZDM 

. 
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VARIABLE 


PROGIiAM 

SYMBOL 


XDM 


YDM 


ZDM 


DESCRIPTION 


Target  velocity  table  selecteci 
is  stored  in  XDM,  YDM,  ZDM 
respectively 


TIME 

PPX 

PPY 

PPZ 

VMX 

VMY 

VMZ 

NPTS 

MISSED 

Till 


Storage  for  missile  position, 
velocity  and  time  in  the  region 
of  pre- specified  range  table 
entries 


\ 


Number  of  points  in  range  table 

An  array  of  values  representing 
miss'  conditions 


Plume  rotation  angle 


Time  = 0 on  target  trajectory 


VARIABLE 


PROGRAM 

SYMBOL 


DESCRIPTION 


X 


x-Missile  minus  x-target  (in 
launch  coordinates) 


y- Missile  minus  y-target  (in 
launch  coordinates) 


z-Missile  minus  z-target  (in 
launch  coordinates) 


VT 


VTI 


Target  inertial  velocity 


XTA 

YTA 

ZTA 


^Lx 

‘“Lz 


Interpolated  from  XE  at  TREAL 

Interpolated  from  YE  at  TREAL 

Interpolated  from  ZE  at  TREAL 

Longitudinal  center  line  of 
target 

Components  of  C 


Line  of  sight  vector 


VARL^BLE 


PROGRAM 

SYMBOL 


DESCRIPTION 


Components  of 


x-Missile  position  in  launch 
coordinates 

y- Missile  position  in  launch 
coordinates 

z- Missile  position  in  launch 
coordinates 

X- target  position  in  launch 
coordinates 

y- Target  position  in  launch 
coordinates 

z-Target  position  in  launch 
coordinates 

X- Position  of  target  in  inertial 
coordinates 

y- Position  of  target  in  inertial 
coordinates 

z-Position  of  target  in  inertial 
coordinates 


VARIABLE 


PROGRAM 

SYMBOL 

XDTGO 


DESCRIPTION 


YDTGO 

ZDTGO 


Initial  values  of  XDTGMS, 
YDTGMS.  ZDTGMS 


Commanded  irradiance  for 
plume 


f 


19.5 


Focal  length  of  projection 
lens  (CM) 


[1^=7000 

€ 


t7 


HR7E 

PlIRSS 

P2IRSS 

AT7 


Irradiance  at  7000  meters 
(W/CM^) 


Vertex  angle  of  plume 


Apparent  altitude  of  plume 
(IN.  ) 


Area  of  plume  transparency 
(IN?  ) 


J 


tul,  plume 


PJTUl 


Radiant  intensity  used  as 
function  of  plume  trans- 
parency area  taken  from 
calibration  in  IRSS 


T 


f,  plume 


0.  3 


Neutral  density  transmission 
factor  for  plume 


VARIABLE 

PROGIiAM 

SYMBOL 

DESCRIPTION 

DXG 

DYG 

j Dummy  storage  location  for 
/ current  value  of  DX,  DY,  DZ, 

1 XDOT,  YDOT  AND  ZDOT 

1 respectively.  These  variables 

I are  returned  to  the  main 

i 

1 

DZG 

1 program  for  diagnostic 
( purposes 

1 

" ' 1 

! 

XDO 

) 

YDO 

ZDO 

LAUNCH 

Array  of  DAC  variables 

A DIN 

Array  of  ADC  variables 

IWRITE 

Flag  for  writing  comments  if 
IWRITE  = 1 

1 

1 

I 

LEVEL 

Status  of  maneuver 
-7  implies  not  in  real  time 

0 implies  in  real  time 
+7  implies  target  traj  table 
exceeded 

1 

WMAN 

1 

Flag,  =1  implies  collect  data, 

= 7 implies  don't  collect  data 

1 

MAN 

Array  for  storing  target 
manuever  table 

MAX 

Maximum  points  collected 

2.27 
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VARIABLES 

PROGRAM 

SYMBOL 

DESCRIPTION 

FI.  F2.  F3 

Intermediate  variables  in  the 
computation  of  COS  (f  ) 

G1,G2.G3 

1 

RRR 

Range- to-go 

— 

IND 

Index  of  XCOMP,  YCOMP 
and  ZCOMP  arrays 

1 

INDEX 

Index  of  range  table  entries 

KCK 

Flag,  =1  implies  t>  t',  - -1 
implies  t S t' 

ICR2 

Counter,  initial  value  = length 
of  array  TS  minus  one 

NPTS 

Number  of  points  in  range 
table 

IDAS 

Logical  variable,  = -1  implies 
not  yet  in  3-loop,  =1  implies 
first  pass  in  3-loop  mode 

XMISS 

Array  to  save  ADCS  under 
miss  condition  (See  MISSED) 

In  some  cases,  program  symbols  are  used  for  intermediate 
storage  or  may  have  scale  factors  other  than  those  implied.  For 
further  information  see  Appendix  A. 
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2.  2 HEAL  TIME  I/O  CONSIDERATIONS 


Phase  1 conversion  of  the  STINGER  real  time  program  requires 
real  time  analog  and  discrete  communication  between  the  AD/4 
analog  computer  and  the  CDC/6600  digital  computer.  This  high 
speed  communications  link  is  provided  by  the  DADIOS.  Specific 
requirements  include; 

• 10-Analog  to  digital  channels 

• 11-Digital  to  analog  channels 

• 1-Discrete  control  line  (2  bits  required) 

• 1-Discrete  sense  line  (2  bits  required) 

In  Figure  2.3  this  arrangement  of  I/O  is  depicted  graphically. 

The  actual  assignment  of  these  analog  and  digital  signals  is 
described  in  Tables  2.  2 - 2,  5.  These  tables  provide  a cross  reference 
of  trunkline  assignment,  variable  name  and  scale  factors. 

Replacement  of  the  IBM/ DOS  interface  with  new  DADIOS  equipment 
required  scale  factor  modifications  to  both  the  ADC  and  DAC  channels. 
This  modification  was  implemented  in  the  digital  computer  code  as  a 
constant  scale  factor  of  1/1,  02375  multiplied  times  each  ADC  channel 
and  a constant  scale  factor  of  . 102375  multiplied  times  each  DAC 
channel.  The  rationale  for  this  scale  factor  change  is  presented  in 
Figures  2,4  and  2,5. 

In  Figure  2.  6 the  DADIOS  patching  required  for  the  STINGER  simu- 
lation is  presented.  This  figure  ndicates  that  a total  of  32  analog 
lines  are  available  although  fewer  lines  are  actually  used.  Figure  2.6 
similarly  depicts  a total  of  32  discrete  line  available,  where  in  fact 
only  four  lines  are  currently  utilized. 


f 
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Preliminary  checkout  of  the  DADIOS/CDC-6600/DDS/AD-4  system 
can  be  performed  in  an  open  loop  test  with  the  STINGER  real  time 
skeleton  test  program.  This  program  provides  a check  on  the  operational 
status  of  all  equipment  involved  in  the  phase  1 conversion  task,  A 
description  of  this  program  is  given  in  Appendix  D. 


The  DOS  provided  the  following  scaling 


IBM- 7094 


10. 2375 


AD/4 


where  -10<x<+10 


Now,  the  CDC-6600  must  provide  the  following  scaling 


CDC-6600 


added  scale 
factor 


Therefore,  digital  output  information  should  be  premultiplied 
by  . 102375  before  being  presented  to  the  DACs, 


FIGURE  2.4  DAC  SCALE  FACTOR  ADJUSTMENT 
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The  DOS  provided  the  following  scaling 


Now,  the  CDC-6600  must  provide  the  following  scaling 


added  scale  factor 


Therefore,  each  ADC  should  be  premultiplied  by  (1/1.  02375) 
before  being  used  by  the  STINGER  real  time  digital  program 


Figure  6,  ADC  scale  factor  adjustments 


FIGURE  2.5  ADC  SCALE  FACTOR  ADJUSTMENT 


2.  33 


TABLE  2.2 

PHASE  1 CONTROL  LINE  ASSIGNMENT 
(CDC  OUTPUT  DISCRETES) 


CDC-6600 

VARIABLE 


CTRLO- 

CTRL3 

CTRL4 


CTRL5, 

CTRL6 

CTRL? 


CTRL8- 

CTRL15 


CDC>6600 

BIT 

AD/4  # 2 
TRUNK 

LINE  i 

PURPOSE /ACTION 

0-3 

] 

TR57-TR54 

Not  in  use 

4 

TR53 

Is  set  high  after 

ICs  have  been  sent 
to  analog.  It  is 
also  the  signal  that 
starts  the  analog 
(SYS-OP)  for  MICOM 
hybrid 

5,6 

TR52,TR51 

Not  in  use 

7 

TR50 

i 

Is  set  high  when  the 
simulation  is  term- 
inated and  the  analog 
should  go  into  system 
hold.  Termination  can 
occur  for  "missed 
target"  or  "all  data 
collected" 

1 

8 - 15 

; TR37-TR30 

L 

Not  in  use 

f 


r 

1 

I 


TABLE  2.  3 

PI-IASE  1 SENSE  LINE  ASSIGNMENT 
(CDC  INPUT  DISCRETES) 


CDC-6600 

BIT 

AD/4  # 2 
TRUNK 

LINE 

PURPOSE /ACTION 

0-3 

TR47-TR44 

Not  in  use 

4 

TR43 

Not  in  use 

5 

TR42 

Set  high  by  the  static 
check  OK  switch  on 
AD/4.  This  but  must 
be  high  to  send  ICs 
to  analog 

6 

TR41 

Set  high  by  the  static 
check  OK  switch  on 
.the  AD/4  for  MICOM 
Hybrid,  or  set  high 
by  the  ramp  up  circuit 
for  IRSS  hybrid  when 
initial  IRSS  conditions 
are  attained 

7-15 

TR40  and 
TR27-TR20 

Not  in  use 

PHASE  1 ADC  ASSIGNMENT 
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ANAI.OG  AD/4  D!G]TA1>  CDC/6600  AD/4 

VARIABLE  ANLG.  FRACTION  ADC  TRUNK 

name  voltage  -i<n<+i  assignment  line 
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PHASE  1 DAC  ASSIGNMENT 
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DIGITAL  DIGITAL  I AD/4  CDC/6600  AD/4 


FOR  MICOM  HYBRID  APPLICATIONS 


3.0  ADC  CHECKOUT  PROG  HAM  FOR  AD- 4 #2 

3.  1 PURPOSE 

The  program  reads  16  ADC  channels  (220-363)  on  AD-4  number 
2 and  compares  the  values  read  to  either  of  two  references:  1.  stored 
references  entered  with  the  program  or  2.  references  obtained  by  a 
single  scan  of  the  ADC  channels.  Selection  of  this  option  is  under 
sense  switch  control.  When  a value  read  back  is  outside  a tolerance 
specified  by  the  user  (default  value  for  tolerance  is  1 bit)  the  program 
will  record  an  error,  nrte  its  magitude,  and  record  the  actual  number 
read  in  a high  or  low  failure  buffer.  The  program  prints  out  the 
accumulated  error  data  every  30  seconds.  A sample  printout  preceeds 
the  program  listing. 

3.2  USAGE 

A.  To  load  program  into  memory  do  the  following: 

1.  Put  Binary  program  deck  behind  universal  loader. 

2.  Ready  card  reader  and  printer. 

3.  On  the  console,  push  Idle  then  Reset. 

4.  Clear  register  lights. 

5.  Hold  down  Clear  and  Clear  Flags  for  1 second. 

6.  Set  Sense  Switch  4. 

7.  Press  Reset,  Run,  and  Cards. 

B.  To  execute  the  program  which  has  been  previously  loaded 
do  the  following: 

1.  Press  Idle,  Reset. 

2.  With  register  display  set  to  B,  enter  00160200. 

3.  Clear  AD-4  Interface. 

4.  Press  Run. 


'1 


C.  To  select  reference  mode,  set  SS2.  This  will  cause  the  program 
to  compare  values  read  to  those  stored  in  the  buffer  beginning  at 
location  65240.  Default  values  for  these  references  are  zero. 

D.  To  halt,  set  SS6.  To  continue,  turn  SS6  off  and  hit  Run. 

E.  To  suppress  print  out,  set  SSI. 

F.  To  modify  tolerance,  load  value  desired  in  location  65561. 

G.  To  modify  print  period,  load  value  in  RCOUNT  (e.  g.  20,  OOOg 
prints  out  every  half  hour)  at  location  65560. 


3.3  HARDWARE  USED 

• SDS  9300  computer,  printer  and  card  reader 

• AD-4/9300  Interface 

• AD-4  number  2 


FIGURE  3.  2 ADC  CHECK  PROGRAM  FLOW  CHART 

3.  6 


COMPUTE  2 X TOL 
STZ  PRINT  COUNT 
ZERO  OUT  ALL 
fiUFFERS  AND  FLAGS 


APPENDIX  A 


CROSS  REFERENCE  OF  REAL  TIME  CODE 
AND  RELATED  EQUATIONS 


DYNAMICS.  INC 


t 


COMPUTATION  OF  IRIS  RATIO  NO.  7 


RLB  = ( fi  / B)  / Scalet 


Elll  = COS  (e  ) 


RFEET  ^ R. 


HR7E  = H 


R=7000  (interpolated  from  table  HR7EE) 
f 


HC7FT 


rHR=700o1 

22965.331  [ e J ' 


plume 


H 

c7 

f,  plume 


- -1  . .5  . -1  , .04884004884 

PlIRSS  = Pj.  2tan  ( = 2tan  < ( g/s)  / o'oiiiT’ 

Note:  .04884004884  = .5/10.2375  and  B = 3, 0 

P2IRSS  = P„  = 22.  5 (10.  2375)  ( C«/B)  / 10.  2375)  / R (inches) 


Pg  sin  (P  ) 

AT7  = 2 2 ~ ■^t7 

2 cos^  (Pj/2)  ” 


. 1.003258  „ , 1. 0032581nA.„  . 

PJTUl-O.lA^'  =0.1e 


/V 

T , 
f,  plume 


RN=( 


iw 

T 

f,  plume 


-)  / J 


tul,  plume 


COMPUTATION  OF  SCALED  PLUME  ROTATION  ANGLE 


FI  = L 


F2  = L 


F3  = L 


G1  = (cos  (ft  ) + sin  (ft  ) 


G2  = (cos  ( ft  ) + sin  ( ft  ) 


G3  = Z^^  cos  (ft  ) - sin  ( ft  ) + Y^^ 


VTI  = V„ 


f 


RCl  = K I?  + l} 

X y 


RCl  = L ! \l}  + 1} 
y ' X y 


cos  ( P ) 


RCB  = cos"^  ( (3  ) = (3 


RCX  = -cos  ( 13  ) = C 


RCY  = -sin  ( 0 ) = C 


IF  (L  > 0)  . C = sin  (0  ) 

X ' y 

FIl  = (Ly  cos  («  ) -sin  (o  )V 


-L.  (Lc,  cos  (or  ) + sin  («  ) + Y^g)  ) / V, 


z '"TE 


A. 2 


J 


F22  = (L^  (cos  ( o ) + sin  ( a ) / Y 

■^x  ^ ^^^TE  ^TE  ^ ^T 


• f ‘ 2 ^ * 

F33  = (L  (cos  (o  ) + sin  (a  ) Z /y  X^^  + Y^^)  Y^^ 


~L^  (cos  { ot  ) +sin  ( or  ) / V ^ ^TE^  ^ ^T 


Note;  C 


Cj^  - ij,  X^j,  (cos  (o  ) + Z^g  sm  (ft  ) / Vx^g  + ^ ^ 


TE  TE'  ' T 


+3g  (cos  (ft  ) + Z^^  sin  (ft  ) / + Y^^)  / V^ 


+k„  (Z^„  cos  (ft  ) - Vx^E  "^TE  sin  (ft  ) ) / V, 


E TE 


Therefore: 


Fll  = (L  “ L C ) 
y Lz  z Ly 


F22  = (L  C,  - L ) 
z Lx  X Lz 

F33  = (L  - L C.  ) 

X Ly  y Lx 


FCR  = V ' Vlz  ■ * ‘L‘=Lx  ■ L'^Lz>  * ‘L'^Ly  ' Vlx’' 


Gil  = (L  C-  - L C ) C 
y Lz  z Ly  x 


G22  . 

(L  C.  _ - L_C,  C.,  + (L.Cl^  - 


^^y^Lz  ~ ^z^Ly^  ■'■^^z^Lx  “ ^x^Lz^  ''■  ^^x^Ly  ^y^Lx^ 


cos  ( e ) 


.3 


-1 


Till  = cos  (cos  e ) - 


^^os^^^  ZE  - COMP  ^ 


ZE  - COMP  < 0 and 


< Jr/2,  Tj^p  = d+  »/2 


^^os^  ZE  - COMP  < 0 and  6 >.  <2.  T^p  = fl-3  »/2 


SPO  = T _ / Scale  Factor 


COMPUTATION  OP^  , COS  (f  ) AND  R/g 

SKK  = K 

52  = cos  e cos 

SPL  = sin  4’  j 

Ij 

53  = sin  dj  cos  '(> 

Li  Li 

54  = cos  sin  ^ ^ 

CPL  = cos 

55  = sin  e sin  if- 

L(  Li 

CTL  = cos  6 

Li 

VTI  = VT 

A1  = (.01745329)  (4.637084242)  yfH  / ( VT  C^^) 

SA  = sin  ( a ) 

CA  = cos  ( ot  ) 

XC  = (200)  (102.  375)  /K 

T.  Jbi 
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YC  = (200)  (102.  375)  /K 

1 Xli 


ZC  = (200)  (102.  375)  /K 

Til 


XTA  = XE 


YTA  = YE 


ZTA  = ZE 


XXX  = X, 


YYY  = Y, 


ZZZ  = Z, 


interpolated  values 


RRR  = x2 


RRR  = V.J,\  x^+  Yj^  + Zq 


FI  = X^  cos  cos  sin  Z^  sin  0 ^ cos  4' 


L G 


F2  = X„  cos  sin  vt  ^ + Y^  cos  vt  . + Z^  sin  d sin  4'  , 
G L LG  LU  Li  Ij 


F3  = Z^  cos  - X^  sin 


SR  - V^e^’^TE 


SI  = cos  (a  ) + sin  ( or  ) + 

01=Xtj,  cos(«)  + sin(o)jc^jji^^/  ^ * '^TE 

G2  . Y^g  cos  («  ) + sin  (o  ) / V =^TE  * 4e 

G3  = Z^g  cos  (or  ) - sin  ( o ) \X^g  + V^g 

Elll  = |(Xq  cos  e ^ cos  “ Y^  sin  sin  0 cos  ^ j^) 

(X^j,  cos  ( a ) + sin  ( ot ) 

V^TE  '^TE 

+ (X_  cos  fi  . sin  T + Y/-,  cos  (f-  _ + Z^  sin  ® sin  >!>  ^ ) 

G Lj  LjG  L/G  \-i  Lt 

(Y^^  cos  ( or  ) + sin  ( or  ) Y,pgZ^^) 

V^TE  ^TE 

+ <Zg  cos  » ^ sin  « g)  (Z^g  cos  (o  ) -sin  («  ) V^E  * ^Te' 

/ C X^  + Y^+Z^l=  cos(o  ) 

COSE  = 1.  0 ~ Elll^ 

RLB  = 1.0-  Elll^ 

RLB  = ( \1.0  - Elll^)  (RLBK)  / Scale  Factor 
COSE  = Elll  / 1.02375 

A.  7 
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sin  ay  cos  tfig  y^  '^\  ' ^ 

T sin  ctj.  sin  V>^ 


FIGURE  A.  1 

TARGET  FIXED  COORDINATE. SYSTEM 


TARGET 


FIGURE  A.  2 
RANGE-TO-GO 

CONTROL  OF  A/D  DATA  COLLECTION 


SIMULATE  CONSTANT 

• ESTIMATE  SmiN 

BEARING  COURSE  AND 

CONVERGING  SEARCH 

ESTIMATE  f AT 

ROUTINE  AND 

1000  ft-TO-GO 

CALCULATE  T^IN 

TARGET 
AT  t - 1' 


y y I y 

MIN^ 


MISSILE  AT  t - t' 


*T'M1N 
TIME,  tac 


SEND  t'  AND  7m, N TO 

ANALOG  COMPUTER  TO 
CALCULATE  K 


FIGURE  A.  3 

SEQUENCE  OF  VARIABLE  SCALE  FACTOR  CALCULATION 
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I 


EARTH  FIXED  COORDINATE  SYSTEM  (Xg.  Yg.  Zg) 


FIGURE  A.  4 

EARTH  FIXED  AND  LAUNCH  COORDINATE  SYSTEMS 


2 3 4 S 6 

TIME,  t,  MC 


FIGURE  A.  5 

VARIABLE  ANALOG  SCALE  FACTOR 


A.  12 


APPENDIX  B 

AD/ 4 TRUNK  LINE  DISCRETE  PATCfflNG 


DYNAMICS.  INC. 


00 

10 

20 

30 

40 

1 

50 

60 

01 

11 

21 

31 

41 

51 

61 

02 

12 

22 

32 

42 

52 

62 

03 

13 

! 

23 

33 

43 

53 

] 

63 

04 

i — 

14 

24 

34 

44 

54 

64 

05 

15 

25 

35 

45 

55 

65 

06 

16 

26 

36 

46 

56 

66 

07 

17 

27 

37 

47 

57 

67 

VARIABLE  1 


VARIABLE  2 


Input  to  ^ ^ ^2  43  44  45  46  47  60  61  62  63  64  65  66  67 

CDC-6600  I 11 


Output  from  ^ ^2  53  54  55  56  57  70  71  72  73  74  75  76  77 

CDC-6600  r I I 


FIGURE  B.  1 AD/ 4 TRUNK  LINE  DISCRETES 
(FOR  AD/4  # 1) 


FOR 

(variable 

2 


TABLE  B.  1 

SUMRIARY  OF  POSSIBLE  DISCRETE  TRUNK  LINE  CONNECTIONS 


CDC-66Q0  INPUT  DISCRETES 

W50  --  INPUT  DISCRETE  1 --  COMMON /=^IDIS2/IIN 
\V51  " INPUT  DISCRETE  2 --  COMMON /’MDIS2/2,  IIN 
W52  --  INPUT  DISCRETE  3 --  COMMON  /*IDIS2/3,  IIN 
W53  --  INPUT  DISCRETE  4 --  COMMON /*IDIS2/ 4,  IIN 


POSSIBLE  AD/4  CONNECTIONS 

V50  --  AD/4  #1  LOGIC  TRUNKS  00-07  AND  20-27 
V52  --  AD/4  #1  LOGIC  TRUNKS  40-47  AND  60-67 
V60  --  AD/4  #2  LOGIC  TRUNKS  20-27  AND  40-47 


CDC-6600  OUTPUT  DISCRETES 

W60  --  OUTPUT  DISCRETE  1 --  COMMON  /*ODIS2/IOUT 
VV61  --  OUTPUT  DISCRETE  2 --  COMMON /*ODIS2/ 2,  lOUT 
\V62  --  OUTPUT  DISCRETE  3 --  COMMON /*ODIS2/ 3,  lOUT 
W63  --  OUTPUT  DISCRETE  4 --  COMMON  /^=ODlS2/4.  lOUT 


POSSIBLE  AD/ 4 CONNECTIONS 


V51  --  AD/4  #1  LOGIC  TRUNKS  10-17  AND  30-37 
V53  --  AD/4  #1  LOGIC  TRUNKS  50-57  AND  70-77 
V61  --  AD/4  #2  LOGIC  TRUNKS  30-37  AND  50-57 
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TABLE  B.  2 

SUMMARY  OF  POSSIBLE  ANALOG  TRUNK  LINE  CONNECTIONS 


CDC-6600  ADCS 


WOO  --  ADC  GROUP  1 — COMMON  /=:=ADC1/ADIN(16) 
WOl  — ADC  GROUP  2 --  COMMON /*ADC1  / 17.  ADIN(16) 
W02  --  ADC  GROUP  3 --  COMMON /"■'ADCl /33.  ADIN(16) 
W03  --  ADC  GROUP  4 — COMMON /*ADC1 /49,  ADIN(16) 


CDC-6600  DACS 


WIO  --  DAC  GROUP  1 — COMMON /*DACl/DAOUT{16) 
Wll  --  DAC  GROUP  2 — COMMON/*  DACl/17.  DAOUT(16) 
W12  --  DAC  GROUP  3 --  COMMON/*DACl/33.  DAOUT(16) 
W13  --  DAC  GROUP  4 --  COMMON/*DACl/49,  DAOUT(16) 


AD/4  ANALOG  CONNECTIONS  (ADCS  OR  DACS) 


voo 

“ AD/4 

#1 

FIELD 

0 

ANALOG 

TRUNKS 

10-17 

AND 

30-37 

VOl 

--  AD/4 

#1 

FIELD 

0 

ANALOG 

TRUNKS 

50-57 

AND 

70-77 

V02 

--  AD/4 

#1 

FIELD 

1 

ANALOG 

TRUNKS 

10-17 

AND 

30-37 

V03 

— AD/4 

#1 

FIELD 

1 

ANALOG 

TRUNKS 

50-57 

AND 

70-77 

V04 

--  AD/4 

#1 

FIELD 

2 

ANALOG 

TRUNKS 

10-17 

AND 

30-37 

VOS 

--  AD/4 

#1 

FIELD 

2 

ANALOG 

TRUNKS 

50-57 

AND 

70-77 

VO  6 

--  AD/4 

#1 

FIELD 

3 

ANALOG 

TRUNKS 

10-17 

AND 

30-37 

V07 

— AD/4 

#1 

FIELD 

3 

ANALOG 

TRUNKS 

50-57 

AND 

70-77 

U07 

— AD/4 

#2 

FIELD 

1 

ANALOG 

TRUNKS 

50-57 

AND 

70-77 

U06 

--  AD/4 

#2 

FIELD 

2 

ANALOG 

TRUNKS 

50-57 

AND 

70-77 

B.4 


APPENDIX  C 


COMPUTER  CODE  FOR  PRODUCING 
NORIVIALIZED  RANDOM  NUMBERS 


r 


«iiwcr:i,'TiNF  .xL . xu.Mij.sir-MA) 

nFFCwiPTION 

Tf'IS  ROuTTNF  PPOPwCff  NOPPALI7F.n  RANPn'^  NU^^pFRS  i"IThT'. 
A range  SPFCIFIEn  HY  THE  liSFW. 

INPUT 


1 

XL 

lower  lI'^it  of 

normal 

OUkvF 

? 

XU 

UPPFR  LIMIT  OF 

NORMAL 

OURVF 

T 

MU 

MEAN  OF  normal 

CURVF 

(TYPE  REAL) 

u 

SIG^«A 

STANDARD  DEVIATION  OF 

NORMAL  CURVE 

OUTPUT 

1 

RX 

THf  nopralizfo 

RANDOM 

NUMPFR 

REF:  General  Purpose  Computer  Subroutines 
Report  No.  TR-WS-75-2 
USAMC,  Redstone  Arsenal,  Alabama 
January  1975 


I 


J 


SUf-«OuTINF  NORMAL  ( PX.  XL»  XUt  MU,  SIG^^A  > 

niMFNsION  FZ ( 1?A) 

PF4L  MU 

DATA  oT.nTi ,NTP,NPT,SQ?Pi/.oi , .0^, i?f.,o. 

«;lf  = MU  - S.n  « sigma 
SU‘?  = MU  ♦ 5.n  » SIGMA 

IF(  XL  .LT.  SUS  .AND.  XU  .GT.  SL5  ) GO  TO  10 
PRINT  SOI.  Sl.S,  SU5,  XL,  XU,  MU,  SIGMA 
SOI  format  (*1LIMTTS  FOP  normal  OISTRIPUTION  SHOULO  «F  «E  T»-t- F N*F  1 0 . ? 
1 * AN0<»FlO.?  program  Tf RMlNATFO.o///*  XL  =*F10.?,SX 

? XU  =*  Fio.?,sx*  MU  =*>Fio.?,sx«  sigma  =*nn.?) 

STOP  1?T 

10  IF(MPT.NF.O)  GO  TO  ? 

T?=nTl 
F7 (1 ) =.S 
FC=F7( 1) 

T = 0. 

FO=S0?PI*FXP (-.S*T*T) 

PT?=nTo.S 

NPT=1 

1 T = T*r)T 

F = SO?PI»FXP<-.S*T»T> 

FCsFC  ♦ RT?«<F*FP) 

FP  = F 

IF(ARS(T-T?).GT..0001)  go  TO  1 
NPT=NPT  ♦ 1 
F7 (NPT)=FC 
T?=T?*DT1 

IF (NPT .LT.NTP)  GO  TO  1 
? CONTINUE 
PN  = PANF(R) 

P = PN 

IF(PN.LT..S)  P=1.-RN 
IF (P.GT.F7 (NPT) ) GO  TO  2 
■ ~-O0  3 1 = 1. NPT 

IF(P.GT.F7(I>  > GO  TO  3 
IF(I.FO.l)  GO  TO  4 
Xsl-l 

x=x*nT  1 
xi=x-nTi 

PX  = XI  ♦ (R-F7 (I-l ) > *nTl/(F7 (I)-F2 (I-) ) ) 

GO  TO  S 
3 CONTINUF 
PPIN'T  100 

lOO  format (13H0  ERROR  OISPN) 

STOP 

A CONTINUE 

py.=o. 

CONTINUF 

!F(PN.LT..S)  RXs-PX 

PXs  SIGMA*PX*MIJ 

IF  (PX.GT.XU)  GO  TO  2 

IF  (PX.LT.XL)  GO  TO  2 

PFTUPN 

FNO 


C.2 
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APPENDIX  D 

STINGER  SKELETON  REAL  TIME 
CHECKOUT  PROGRAM 


I 


SKELETON  REAL  TIME  PROGRAM 


EXECUTION  PROCEDURE: 


NOTE:  The  B-K  microprocessor  is  to  be  programmed  and  connected  to 
the  appropriate  trunks  on  the  AD/4  logic  board  as  specified  in  the 
section  "MICROPROCESSOR  REAL  TIME  SET-UP"  to  handle  the 
discrete  communication  with  the  digital  program.  The  symbol 
^ means  the  "SEND"  key  on  the  DDS  terminal. 


1.  Make  sure  all  the  input  discretes  from  the  microprocessor  are  low. 


2.  Catalog  the  SKELETON  REAL  TIME  program  by  submitting  the  deck 
with  the  following  cards  in  front: 


XXXGS,  T5,  CM6000. 

ACCT(PN=XXXXXXX,  PBC  = 7364G13R00,  CC  = 7300.  OP=A3,  JN=XXXX) 
REQUEST,  SAVE,  *Pp. 

COPY,  INPUT,  SAVE. 

CATALOG,  SAVE,  XXX GS,  ID=G1XXXX. 

EOR 


(NOTE:  X's  represent  characters  to  be  supplied  by  the  user. ) 


3.  Sign  on  the  DDS  terminal  as  follows: 

Screen  says  "ACTIVE  AND  NOT  BUSY" 

User  types  A 

Screen  says  "TYPE  USER  NM,  PASSWORD,  PBC,  CC,  PRG  NM,  JOB 


User  types  XXXXXXX,  XXXXXXX,  7364G13R00,  7300,  XXXXXXX,  XXXX  A 


Screen  says  "SELECT  A COMMAND  CODE  FROM  THE  LIST  BELOW" 


O Operator  Aid 
J Job  Control 
S File  Scan 

G File  Generation  and  Updating 
U File  Utilities 
P Permanent  Files 
T Terminate 


4.  User  types  P A to  get  permanent  file  command  set. 

5.  User  types  ATTACH(SKELRT.  XXXGS.  ID=G1XXXX)  A 

where  the  XXXGS  and  GIXXXX  are  thesame  as  on  the  catalog  card  in 
Step  2, 

Screen  says  "FUNCTION  SUCCESSFUL",  etc. 


6.  User  types  $J  /2k  to  get  job  control  command  set. 

7.  User  types  SKELRT,  nn^  where  nn  is  the  first  listed  free  control 
point  (01,  02,  etc. ). 


8.  Screen  says  "SELECT  STEP  MODE  IF  WANTED",  etc. 

User  types  A to  select  automatic  mode.  Program  will  compile, 
load,  and  begin  execution, 

9.  Screen  says  "PAUSE  l"  with  "PAUSED"  on  right  side.  At  this 
time,  all  output  discretes  should  be  high.  This  is  detected  by  the 
microprocessor  as  event  #1. 

10.  User  t3T3es  GO  ^ (a  second  GO  ^ may  be  necessary).  This 
will  clear  PAUSE  1. 

Screen  says  "PAUSE  2"  with  "PAUSED"  on  right  side.  At  this  time, 
all  output  discretes  should  be  low — detected  by  the  microprocessor 
as  event  #2. 


11.  Set  input  discrete  5 high  using  a switch  on  the  microprocessor.  User 
types  GO  A (a  second  GO  A may  be  necessary).  "PAUSE  2" 
will  stay  on  screen  but  "PAUSED"  on  right  side  should  disappear. 
Output  discrete  4 should  come  high  now — detected  by  microprocessor 
as  event  #3. 


D.  2 


A 


12.  After  a delay  of  approximately  15  seconds,  the  microprocessor  will 
automatically  set  input  discrete  6 high  (event  #4).  This  initiates  the 
ADC-DAC  conversion  loop  in  the  program. 

13.  When  all  conversions  are  complete,  the  program  will  set  output 
discrete  7 high  (microprocessor  event  #5)  and  terminate. 

Screen  says  "NO  MORE  CONTROL  CARDS". 

14.  User  types  $U  A to  get  file  utility  command  set. 

15.  User  types  1 ^ to  terminate  and  save  output  file. 

16.  User  types  L,  file  name,  XX,  AU  ^ in  order  to  list  the  output  file, 
where  the  file  name  is  as  listed  on  the  bottom  of  the  screen — 
RAIGSXX. 

17.  User  types  $T  A and  then^  to  sign  off  DDS. 
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MICROPROCESSOR  REAL  TIME  SET-UP 

The  software  checkout  chassis's  microprocessor  is  loaded  with  the 
program  below  using  the  front  panel  switches.  The  procedure  is  as 
follows; 

• Select  JAM  and  STEP  mode  (switches  up) 

• Load  000  into  high  and  low  address  registers  using  LAH  and 

8 

LAL  with  OOOg  set  on  switch  register 

• Load  program  using  DEPOSIT  switch  and  switch  register 

• When  program  is  loaded  set  005g  as  switch  register,  select 
RUN  and  RUN  on  mode  switches,  hit  INTERRUPT 

• Program  will  now  execute 

• At  completion  of  execution,  hit  INTERRUPT  to  re-execute 
The  AD/4  must  be  cabled  to  the  processor  as  specified  below;  • 


AD/4  TR  20 
21 
22 

23 

24 

25 

26 
27 

AD/4  TR  (SP) 


line  8*  Checkout  chassis  DO  X** 

7 6 

6 Connect  to  switch  on  CC 

5 X 

4 X 

3 X 

2 X 

1 X 


0 

1 

2 

3 

4 

5 

6 
7 


line  16  • Checkout  chassis  DI 


7 

5 
1 
4 

8 

6 
3 
2 


* This  cable  is  installed  between  the  AD/4  and  software  checkout  chassis 
**  X Denotes  "don't  care" 
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MICUOPKOCESSOU  UKAL  TIME  PROGIlAM 


LOCATION’ 


OP  CODE 


OPEllATOH 


CONNECTS 


000 

006 

MVl  A 

001 

000 

000 

002 

121 

A->P0 

003 

123 

A->P1 

004 

101 

0-»A 

005 

074 

CPI 

006 

377 

. 377 

007 

110 

JPZ 

010 

004 

oil 

000 

012 

006 

MVI  A 

013 

001 

01 

014 

121 

A->P0 

015 

101 

0— >A 

016 

074 

CPI 

017 

000 

ZERO 

020 

110 

JNZ 

021 

.015 

022 

000 

023 

006 

MVI  A “ 

024 

002 

002 

025 

121 

A-^PO. 

026 

101 

IPO  A 

027 

022 

RAL 

030 

022 

RAL 

031 

022 

RAL 

032 

022 

RAL 

033 

022- 

RAL 

034 

100 

JNC 

035 

026 

036 

000 

• • 

037 

006 

MVI  A • 

040 

003 

03 

041 

121 

A— >P0 

042 

106 

CALL 

013 

100 

DELAY 

044 

000 

• ' 

045 

006 

MVI  A 

046 

040 

040 

047  • 

123 

A-^POl' 

050 

006 

MVI  A ‘ 

051 

004 

04 

052 

121 

A— ^PO'i 

053 

101 

PO-^A 

054 

022 

RAL 

055 

022 

RAL 

056 

100 

JNC 

057 

053 

060 

000 

061 

006 

MVI  A 

062 

005 

05 

063 

121 

A-^PO 

064 

104 

JPU 

065 

064 

066 

000 

Check  for  all  ones 


Display  event  1 


Check  for  all  zeroes 


Display  event  2 


Check  for  bit  4 


Display  event  3 


Delay  15  seconds 


Raise  discrete  6 


Display  event  4 


Check  for  bit  7 


Display  event  5 


L. 
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FIGURE  D.  1 CHECKOUT  CHASSIS  FUNCTIONAL  BLOCK  DIAGRAM 
(DISCRETE  HANDLING  AND  DISPLAY) 
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PRCGf?^.''  TSja'5(  CUTPJT,HFILE,TAFE6=Cl,TPUT) 

USE  ONLY  Cn£  input  DISCRETE  AT  A TI^'£  TC  AVOIC  ERRORS  IK  LCGIC 


NOTE,  ALL  CUIFLT  DISCRETES  SHOULD  APPEAR  ON  DISPLAY  LIGHTS 
FRCGRAK  FCR  NCN  REAL  TIME  COMMUNICATION  OF  OAOIOS  VARIABLES 


PROGRAM  VARIABLES 
lERR  t 


ISTAT 

LEN 

PIMI) 

POUT(I) 

LARA 

N 

lODIS 

IICIS 


LCCF 


ERROR  CODE  FOR 
r=NOERRCR,  GT.O 
REAL  TINE  MODE 
0=IN  REAL  TINE, 
NUMBER  OF  TRUNK 
I-TH  TRUNCKING 
I-TH  TRUNCKING 
AN  ARRAY  OF  LOG 
THE  NLM3ER  CF  L 
DAOIOS  OUTPUT  N 
CAOICS  INPUT  Cl 
IF  LINE  1 HI 
IF  LINE  2 FI 
IF  LINE  3 HI 
IF  LINE  M HI 
NUMBER  CF  LCCFS 


RESERVATION 
= R£SERVATICN  ERROR 

ISTAT. GT.C  MCT  IN  REAL  TIME 
ING  STATION  CCN'NECTOR  FAIRS  USED 
STATION  CONNECTOR  INPUT 
STATION  CONNECTOR  OUTFIT 
ICAL  INTERRUPT  NUMBERS 
CGICAL  INTERRUPTS  IN  LARA 
ORO 

SCRETE  WORC 

GH  VARY  lOL'T  IN  BATCH  TR27  OR  TP67 
GH  GO  INTO  REAL  TIME  TR28  CF  TR66 
GH  RETURN  TC  BKGRNO  TR25  CR  TR65 
GH  STOP  7R2<*  CR  TR64 

IN  BACKGROUND 


CACIOS  PATCHING  REQUIREMENTS  tCNE  IIOIS  AND  ONE  lODIS) 


TRUNKING 
V-5C  TC  W-50 
V-52  TO  K-5C 
V-EC  TO  W-51 
V-E2  TC  W-51 

V-51  TC  W-6G 
V-53  TO  N-6C 
V-51  TC  W-61 
V-53  TC  W-61 


FORTRAN 

FOR  /IOIS2/l,IIDIS 
FCR  /I0IS2/1, IIDIS 
FOR  /IOIS2/2,IID  S 
FCR  /ICIS2/2, IIOIS 

FCR  /CCIS2/l,IODIS 
FCR  /CCISZ/i ,IODIS 
FOR  /ODIS2/2,IODIS 
FCR  /CDIS2/2 tIODIS 


AD/4  PATCH 

TR00-TR07  AND  TR20-TR27 
TR40-TR47  AND  TR6C-TR67 
TR00-TR07  AND  TR20-TR27 
TR4C-TR47  AND  TR6C-TR67 

TR10-TR17  AND  TR3C-TR37 
TR50-TR57  AND  TR7C-TR77 
TR10-TR17  AND  TR30-TR37 
TR50-TR57  AND  TR7fl-TR77 


THE  HIGH  CREEP  COC-6&QO  BIT  CCRRESFONCS  TO  TROX,  WHERE  X=0,2,4,6 

INTEGER  PIN(8),F0UT(8) 

DIMENSION  LARAI8) 

CCMMC^/1^TCCM/ICNT1 
interrupt (1=1, R=20 ,T=100000) 

CCMMON/*IDIS2/l, IIOIS 
JAM, CN 

CCMmCN/*ODIS2/1,IOOIS 
JAM, OFF 

INITIALIZATION 


LCCP-C 
LARA { 1 )=1 
N = 1 


OOO  ooo  C^O  oo  OO  Oo  oo  OOOO  o OOO  OOO 


c 


DC  1C  1=1, LEN 

10  WRIT£(6,8C20)PIN(I),POUT(I) 

JA'l  A CISCRETE  KORD  OUT  ONTO  AD-l*  TRUNK  LINE  IN  BATCH  NODE 

2C  I0UT=2**ICCUNT 
ICCUNT=IC  CLNT+1 
IF (I COUNT. EC. 16) ICCUNT=0 

' IC0IS=I0UT 


READ  BACK  A DISCRETE  WORD  FROM  AD-4  TRUNK  LINE  IN  BATCH  MODE 

3C  CALL  IMCNVT  (LARA,N) 

IN'=IIDIS 

PAUSE  AND  WAIT  FOR  APPROPRIATE  INPUT  CISCRETE  TO  EE  SET 
PAUSE 

TEST  LINE  1 OF  IIDIS 
IF(IN.£Q.1E)G0  TO  20 

TEST  LINE  2 OF  IIDIS 
IFdN.eO.  2E  )GC  TO  40 

TEST  LINE  4 OF  IIDIS 
IF  (iM.EQ.lCfi)GO  TO  50 

COUNT  NUMBER  OF  LOOPS  IN  BACKGROUND  MODE 
LOOP=LOO?  +1 
GO  TO  30 

RESERVE  EQUIPMENT  AND  ENTER  REAL  TIME 

40  CALL  STMRLN( ISTAT) 

V.'PITE  {t,2CCC)  ISTAT 

CALL  R£MARK(17M  JOS  IN  REAL  TIME) 

HCLC  GACKGFCUNC  UNTiLL  SIMHCLC  OR  SIMSTOP  IS  INCOUNTEREO,  ON  RETURN 
TO  BATCH  TURN  ALL  OUTPUT  CISCRETES  OFF 

E.6 


1 


8000  FORMAT  (/♦OTRUNK  LINE  PATCHING  IS  *»R3f*  TO  ■*,R3) 

ENO  ,1 


oo  ooo  o or>o 


SlG^^CUTIN'E  SUOl 

REAL  TIHE  SUBROUTINE 

CC'^MON  /INTCCM/ICMTl 

CC^'^'C^/♦IClS2/l,rIDIS 

ccy^'C^/»ccIS2/l,IocIS 

IN=IIDIS 

IF  IN  REAL  TI^E  TURN  ALL  CL'TPLT  CISCRETES  ON 


IOCIS=2**16-l 

TEST  LINE  3 OF  IICIS 
IF(IN.EO.i»B)CALL  SIMSTCP 
C 

C TEST  LINE  L OF  IICIS 

IF(IN.EQ.1CE)CALL  SIMSTOP 
C 

CALL  SIKIOLE 
END 


RT REE  TRJAN (C ) ,SUB1  (1  ) 
global  IMTCOM 

END 


RESERVATION  ERROR  CODE=  QO  0 S C C C C 0 C 0 C C C C D 0 000 


TRUNK  LINE  PATCHING  IS  V50  TO  W5C 


TRUNK  LINE  PATCHING  IS  V51  TO  W6C 


TRUNK  LINE  PATCHING  IS  V52  TO  KSl 


TRUNK  LINE  PATCHING  IS  V53  TO  W )F 
REAL  TINE  ST ATUS=0 0 Q 0 000 0 0 00 00 0 0 00 0 0 0 
TOTAL  LOOPS  IN  BACKGROUNO= 

PROGRAM  TERMINATED  NORMALLY  WITH  LINE  A HIGH 
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APPENDIX  F 

COMPUTER  LISTING  OF 
PHASE  1 STINGER  CONVERSION 

♦ • 
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I 

I 

I 

I 
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DYNAMICS.  INC. 

I 

• f 


MAIN 


74/74 


QPT  = 1 


FTMH  4.2475075 


07/16/75 


— Pr^OG'^A'i  M A I N ( I NPU  T , OUTPU  T » hF  I L T AP£6  = 0U  TPUT  ^ T APE5  ■=  INPUT  » T APE8> - 

CCMMLNT. 


COMMtNl , 
COMMENT. 

C*** 

C*** 

-c  NP.U 

C XT£ 

C QX 

C NPT 


hY5'?Il)  simulation  control 


NP.UNS  Is  [m£  number  of  runs  PER  FLIGHT  CONDITION.  MAXIMUM  OF  20 

XTERM  is  Th£  VALjE  OF  DELTA  X AT  WHICH  SIMULATION  IS  TERMINATED 
QX  IS  THE  DELTA  RANGE  BETWEEN  COLLECTED  DATA  POINTS 
NPTS  IS  THE  NUMUER  OF  DATA  POINTS  COLLECTED 

COMMCN/COMA/LEVEL ,IPTS,XXS(50) , XDT GO . Y OT GO , ZDT GO . R L 3 , COSE . SPO, RI , 
♦GAM.EDOT, THETAL.RN,  --  - - 

* P°X (50) ,PPY(50) ,PPZ(50) .TIME (50) , TM A S ( 30 ) , XD TGMS ( 3 0 ) 

* ,YurGMS(30) .ZUTGMS(30) ,XMAN{4,50) ,XMISS(7),NT 

* ,XCOMP,YCOMP,ZCCMP,TAMA(3  0) , DEL  TAR ( 30 ) , VM (3 O ) .G.GGG 

» V 'in  vnn  Jnn  n v r n v n n 7 n ^ c O C T c:  . y OM  / ■?  n \ . V HM  / ■?  n » - 7 OM  f fl  I . 


►.ECOSdZ)  .HR7EE(17)  , F I 7T  ( 9 ) , FC I 77  ( 9 ) 

COMMON  /C0M3/XD, YD.ZD, XD  TO . YDTG.ZOTG ,XDEF , YDEF  » ZOEF, ZZZ 
►.XJOE.YJUE.ZDDE 

COMMON/M ANEUV /I TT ,X OEM (30,45) , YOEM (30,45) ,Z DEM (30, 45) , 

* XDTGM(30) ,YDTGM(30) ,Z0TGM(3D) , XPRIME, YPRIME , 

* ZPRIME.TIMES  (50)  ,ITRAJ  (60  ,ITAU  (60)  , III,NAOJ 
CCMM0N/GrARG/FXA,FYA,FZA,FX3,FY3,FZB,FXC,FYC,FZC,AT,VT 
COMy.ON/COMP/ISl,I99,I83,CLMS(10)  , CL  SI  (10)  ,CLS2(10)  ,NSAM,Q1, 

* . UM  . - wu  - Wn  . Wf,  - TFT  . TPL  ANF.  TOON  . S2P1  .SlPl  , A2P1.  A 


CCMMCN/GrARG/FXA,FYA,FZA,FX3,FYU,yxci,(-XU,hYU,t-ilU,Ai,vi 

COMy.ON/COMP/ISl,I99,I83,CLMS(10)  , CL  SI  (10)  ,CLS2(10)  , NSAM,  Q1  , Q2  , Q3  , 
♦C4,a5,C3,Wl,W2,W3,W4,W5,W6  ,IFT,IPLANE,IC0N,S2P1,S1P1,A2P1, AlPl, 
♦■S2F2,Sl?2  ,A2P2  , A1P2,S2P3,S1P3,  A2P3  ,A1P3,S2P4,S1P4,  A2P4,A1P4,S2P5,- 
*S1P5, A2P5,AiP5,S2P6,SlP6,A2P6, A1P6 

- C0MM0N/EXTRA/NPTS,DX,XT£RM,KKK,AS,S£0,SC ,ASE0,R,SCALEL,VVVV»C1, " - 
♦XLOS, YLOS,ZLOS,CO,CX,CY,CLY,CLZ,TR3,XX (50) ,ThE 

cTM''rY  r'nccn.cTM‘CPr^_Pf^<‘Qf^O  - 


*S1?5, A2P5,AiP5,S2P6,SlP6,A2P6, A1P6 
- C0MM0N/EXTRA/NPTS,DX,XT£RM,KKK,AS,S£0,SC ,ASE0,R,SCALEL,VVVV»C1, " 
♦XLOS, YLOS,ZLOS,CO,CX,CY,CLY,CLZ,TR3,XX (5  0)  ,ThE 
COMMON  /ANG/  S INS 0 , COSSO , S INSCO , COSSCO 

C3:!M0\/Pi</XDH(22)  , Y J H ( 22  ) , ZJH  ( 22 ) , 8DHX  ( 2 2 ) , 3DHY  ( 22  ) , BDHZ  ( 2 2 ) 

—♦.RILL  (22)  ,P(4) ,SPK,SPHG,SPHS,SOM 

COMMON/M  I SD/XM ( 3 ) , YM  ( 3 ) , ZM  ( 3 ) , RM ( 3 ) , NOPT , I ON 
■ - CCMMON/JELT/PPXS  (50  ,PPYS(50)  , PPZS(50)  ,VMXS(50)  , VM  YS  (50)  . VMZS(50) 
COMMON/0  3MD/COSA, C05  3,C0SG,YA,  ZA 


REAL  NR 
u I y N u 1 0 N 
DIMENSION 


ZrE0(62)  ,TM(30  

AXM (20 ) ,AYM  (2  0)  ,AZM (20) , ARM (20)  , 
AXB(20) ,AY3{2  0)  ,AZ3(20) , ARB (20)  , 
•i7n(2Q)  .ARC  (20) 


* AXB(20),AY3(Z0),AZd(2U),AS;B(21J), 

* AXC(20),AYC(20),AZC(20),ARC(20) 

DIMENSION  ASC (62) ,RR(62) ,XMT (3) ,YMT (3) ,2MT(3) ,RMT(3)  - - 

DIMENSION  AXMT  (20 ) ,AYMT(20 ) ,AZMT (2  0 ) , ARMT ( 20 ) , AXBT ( 20 ) ,A YBT ( 20)  , 

♦AZ3T(2'j)  ,ARBT(20)  ,AXCT(20  , AYC  T ( 20  )“,  AZCT  ( 20 ),  ARC  T ( 20  ) 

INTERRUPT (I=1,R=50 ,T=80,P=10Q) 

CnMMnN/^.iiri/LAUMCHdl)  — 


liNic.rcfv.uri  li  — 

COMMON/^UACl/LAUNCHdl) 
C0MM0N/»ADC1/ADIN (10) 

COMMON/* IDIS2/IIN  - 

JAM, ON 

-C0MMCN/*00IS2/I0UT 

JAM,  OFF 


( 


NKr  CARO  IS  TAPE  SPECIFIC  FUR  LAUNCH  AT  WSMR 
DATA  ZAnOVE/-30S9./ 

UATA  RI  0,UTk/07. 2957795 ,0.01745329/  

UATAiTM(r),I=l,27)/0.,l.,2.,3.,4.,5.,6.,7.,3.f9.,10.fll.tl2.t 
* 13«,14.,15.,16.tl7.,18.,19.t20.,21.,22.,23.f24.»25.f26«/  — 

SCALE  = 1C2.375 

oCALtr  = 102.375  / lU.  ' 

SCALEL  = 200.  * 102.375 

SCALEV  = 20.*102.375  ^ 

SCALE2=b. ♦1.02375 

bCALED  = 20.*102.375  '•  

K.ASE=0 

RL0K  = 3.3333  

10  CONTINUE 

READ  (5,2000)  IS1,NSAM,IC0N,INL  — 

WRITE  (rt,S0Q) 

KKITl  (3,2002)  IS1,NSAM,IC0N,INL  

PRINT  300 

DO  11  1=1,10  ^ 

REAO  (5,620)  CLMS ( I ) , CLS 1 ( I) , CLS2 ( I ) 

PRINT  2005,  CLMS(I),CLS1(I),CLS2(I) 

11  CONTINUE 

REAO  (5,621)  Ql, 02,03,04,05,06  

PRINT  300 

. 2006,  01,02,03,04,05,06  

REAU  (5,621)  W1,W2,W3,W4,W5,W6 

PRINT  30  0 - 

PRINT  2007,  W1,W2,W3,W4,W5,W6 

- call  oclim  

READ  (5,902)  NRUNS 

— print  30  0 - - 

PRINT  2010,  NRUNS 

NR  = NR  U MS 

IF  (NRJNS.LE.20)  GO  TO  15 

PRINT  80  2,  NRUNS  . ^ 

NRUNS  =20 

j_5  CONTINUE  

READ  (5,903)  NPTS , DX , XTERM 

PRINT  30D - 

PRINT  2008,  NPTS, DX, XTERM 

LAUNCH  CONDITIONS.  NL  IS  NUMBER  OF  LAUNCH  SETS. 

READ  (5,902)  NL 

— PRINT  30  0 ^ - — - . - 

PRINT  2011,  NL  l! 

IF  (NL.LE.61)  GO-TO-  30'- r"~"r'~7r~ 

PRINT  323,  NL  ' ' ''  • ■ ■ - 

GO  TO  3 — 

30  CONTINUE  ji 

PRINT  30  0 - - - — !■: 

IF(ISl.NE.l)  PRINT  333  h 

C***  read  ANU  PRINT"L0ADE0-LA'JNCH  CONDITION-SETSi 'i 

NLL  = NL  + 1 

PRINT  30  0 - j 
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—  DO  40  1=1, NLL  

k£A3  (0,904)  ZTEO (I) , ASC (I ) ,RR(I) 

ZTE  = -ZTEOil)  ~ 

IF(ISl.NE.l)  PRINT  835,  I , ZTE , ASC ( I ) ,RR ( I ) 

40  CONTINUE  ' -■  — — 

3 CONTINUE 

—  PRINT  300  - 

NT  = 27 

iaci=o  — 

GO  TO  42 

— - C***  AUTOMATIC  SELECTION  OF  LAUNCH  CONDITIONS.  

41  INL=1NL*1 

C INDEX  for  AUTO  TEST  CASE  CHECK 

63  I99=u  ' 

U77=FLOAT (INL) 

— — - U73=  (U77-1.  )/5.  

T77=IFIX (U78) 

U79=FLCA r (177)  - ^ - 

U30=U73-U79 

— C test  for  every  fifth  LAUNCH  SET  ^ 

IF(U80.LT.C.01)  199=1 

C TEST  FOR  FIRST  LAUNCH'SET  - * 

IFdKL.EQ.l)  199=1 

—  IF  (INL-.ER.l)  IUC1  = 1 

C INDEX  FOR  OTHER  RUNS 

188=0  - '*  

C TEST  FOR  LAST  TEST  CASE  CHECK 

IF  (INL.GT.NL)  133  = 1 

C TEST  CASE  CHECK  REQUIRED  AFTER  LAST  LAUNCH  SET 

IFCINL.GT.NL)  199  = 1 

IF(I99,EQ.l)  GO  TO  62 

^ SZS-'KASE  = INL  

KA=KASE 

PRINT  918,  KA 

ITT=ITRAJ (INL) 

PRINT  916, ITT 

IT1=ITAU (INL) 

• PRINT  917, ITl - ' 

TAU=FL0AT(IT1) 

GO  TO  51  - . - 

42  READ  (8,2000)  III , IRSS ,KA, ITT, NCGO 

v:RITE-(8 ,800) 

>(RIT£  (8,2003)  III, IRSS, KA, ITT, NOGO' 

. AU  (8,  2001)  - TAU  

WRITE  (3,800) 

WRITE  (3,2004)-TAU  

IF(III.EQ.l)  GO  TO  41 

-50  CONTINUE 

5 PRINT  300 

23  continue — 

C***  LAUNCH  CONDITIONS  HAVE  BEEN  SELECTED. 

— GO  TO  51  - ' * 

62  CONTINUE 

PRINT  62  5 

TAU=3.0 

C INJi-X  FOR  TEST  CASE  DATA  CARD  OF  LAUNCH  CONDITIONS  — 
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C INUfc'X  FO,^  TtST  CASt  TARGET  MANEUVER 

C***  NEXT  CARJ  IS  TAPE  SPECIFIC  FOR  ITT=39  FOR  TEST  CASE  MANEUVER 
ITT=59 

C***  ADJUST  FOR  TAU  AnO  SCALE  TARGET  VELOCITY-TIME  MANEUVER  TABLE. 
51  NTAU  = INT(TAU)  + 1 

—  DO  25  I = NTAU,  NT 

X - I - NIAU  + 1 

TMA{i<)=TM(I)-TAU  _ - --  

TMAS(K)=  TMA (K)/SCALET 

— XQM(K)=X0EH(I,ITT)  

YUM(<)=YJLM(I,ITT) 

— ZDM(K>=ZDEM(I,ITT) 

25  CONTINUE  ^ 

KK=NT-INT  (TA'J) +1  — 

CO  2o  K=KK,NT 

— --  TMA<kJ=FLOAT(K)-1. — - — 

X0M(K)=XDEM(NT,ITT) 

YQM(K)=YD£M(NT,ITT) 

ZJM(K)  = Zl)EM(NT,ITT) 

— 26  CONTINUE  - • * 

generate  TARGET  AERODYNAMIC  ACCELERATION  TABLE 

— - XE(1)=XDN(2)-XDM(1)  

Y£(1)=YJM{2)-Y0M(1) 

ZE  (1)=Z0:1  (2) -ZOMt  11-32.174 — 

no  31  1=2, NT 

XE  (I) =XjM(I)-XOM(IV) 

{ YE  (I)=YCMm-YCM(IV)  ' 

ZE  (I)=Z0N{I)-ZDM(IV)-32. 174 

; — : 31~C0NTINUE 

I KASE=KA 

—  Z=ZTEOtKASE) ' " 

ZZZ=Z+ZABOVE 

I ZALT=2ZZ  

AS  = ASC(KAS£) 

I — R=RR(KASE) 

II  XD=X0M{1) 


d 
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GO  TO  t)J 

5 3 l!ETA  = 9a.  

GO  TO  5o 

54  aETA=270.  

GO  TO  5b 

- 55  u£TA=ATAN2(YD,XD)*RTD 

56  ASC0  = A3+ t'ETA-ldO. 

IF  (A:iCO.GT.O.O)  GO  TO  60  — - 

ASCO=obii . + ASCO 

60  CCNTll.lJt  

SC0=ASCJ*DTR 

— ARb  = -Z/R  

1 50-ASIN(ARG) 

1 1 riETAL  = GO*RTL)/SCALE  

ASO=SO*RTD 

--  - - X=R*COSCSO)  *005  (SCO)- 

Y=R*30G(  J.O)  •GIN(SCO) 

PRINT  805,  AjCO,ASO,R 

PRINT  'ai2,X,Y,2 

- — CALCULATE  CONSTANTS  FOR  COROINATE-TRANSFORMATION 

iI.N30  = sr  NTSO) 

COSSO=COS (SO)  

SINSCO=SIN(SCO) 

- COSSCO  = COS(SCO) — 

S2=CCSSO*COSSCO 

_ S3=SINSO*COSSCO  ’ 

S4=C0SS0»SINSC0 

S5=SINS0*SIN3C0  

bPL=SIMSCO 

_ CTL  = C0S30 

; ' STL=SINSO 

: CPL  = COSSCO  ' 

C**»  make  earth  TO  GUIDANCE  TRANSFORMATION 

; _ . CALL  GUia  

XDTGO  = XDTG 

YOTGC  =-YDTG 

ZDTGO  = ZDTG 

■ CAL CUlATE  SCALE-FACTOR -EQUATION 

R : = R 

j ..  ^.^LL  KSCALE,RETURNS(50) 

[ EDOT= (GGG/G) /SCALE2 

i;  READY  'TO  -fl AKE-FLIGHTS-ANO-COLLECT-POSTTION-POINT-S- 

^ SPR=Q. 

; . SPriG-Q.  - - 

; SPriS=0. 

' N111=MRUNS 

IF  iIRSS. EQ. 1)-NRUNS=IRSS 

5 IF  (IR3S. -Q.l)  NR=IR3S 

• IFdRSS.EQ.l)  199=0  -- 

IF(I99.Sa.l)  NRUNS=NSAM+1 

IF(I99.EQ.l)  NR  = NRUNS  

I'  DO  150  I-1,NRUN3 

' 1 call-ini  T 

IF  (ISl.NE.l)  PRINT  77, RN 
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C***  2t.RO  SIGNAGES  AREA 

DO  05  J=  1, 5 0 - - 

UO  85  K=  1 f 

aMAN(K,J)  = 0.  - 

65  CONTIUOE 

DO  90  0=1,50 

PPX(J)  = C.O 

— PPY  ( J)  = 0, 0 _ 

PP2(J)  = 0.0 
TIME( J)  = 0.0 
V M X ( J ) = 0.0 

- — \/'1Y(J)  = C-0  - 

VM2(J)  = 0.0 

90  CONTINOZ  - - 

IPTS=*JAOJ  _ _ , — 

GALL  FLIGHT  — 

- .--  -Tfriil.'i-^T.ir'GO  TO  91 

— - PRINT  99^'’,RI,GAM,ED0T,X0TG0,YDT50,Z0TG0,RN,RL3,SPD  

PRINT  m99,XC0MP,YCCMP, ZCOM?,G,GGG,XDO, YOO,ZDO,DXG,OYG 
PRINT  99''),DZG,S2,S3,S4  ,S5,RL3K,SCALEP,C0SE,F1,F2 
PRlNl  999  ,F3,G1,G2,G3, XC, YC,2C,S1, RRR,SR 
- — - PRINT  999,SPL  ,CTL,STL,CPl,A1,VTI,ZALT,CLA,PHO,ARG 

PRINT  99  9 ,AAA,SCALET,TREAL ,XTA ,YTA , ZTA , SCALEV ,QM,S A,CA 

PRINT  99&,LE7£L,IPTS,NT,NERR,NPX,NX 

print  999,TRP,XL0S,YL0S, ZLOS,CO,CX,CY,CLY,CLZ,VVVV 
PRINT  999,C1,AS,SE0,SC,ASE0,R,SCALEL,TH£TAL,THE  - 

91  continue 

NOPT  = NACJ  — - 

C***  SCALE  and  write  THE  POSITION  POINTS 

If  (ISI.Nl.1)  PRINT  8 00-  

IF  (XMI5S  (7)  .GE.O.)  PRINT  909,  XMISS 

IFdSl.Ea.l)  PRINT  809,ICNf 

IFdSl.NE.l)  PRINT  810, ICN  • 

-95  IF  dPTS.EU.O)  GO  TO  150  

JO  110  J=l,NAOJ 

- PPXS(J)  =-PPX(J)  * SCALEP 

PPYS(J)  = PPY(J)  • SCALEP 

- PP2:>(J)  = PPZ(J)  ♦ SCALEP—  — ^ — 

TIHES(J) =TIME(J)  * SCALET 

. - VMXS(  J)  = VMX  (J)  *SCALEL  — 

VMYS(J)=VMY(J> ’SCALEL 

VMZG(J)  = VMZ(J)*SCALEL--  

IFdSl.c.Q.l)  GO  TO  110 

PRINT  Hll,  PPXS(J)  ,?PYS{J)  ,PPZS{J),TIMES(J) 

lie  CO;.TI„JE 

C ChEG:<  digital  TIME  CRITICAL  DATA  COLLECTION  

!.C.<  = riADJ-l 

yj  330  ;-,  = l,NOK — ■ 

XCOHH=-XX(MI 

XCCLL  = PPXS(M)  — - - 

XCCOX=-XX (M+1) 

IF(XGaMM.LT.XCOLL.AND.XCOLL.LT.XCODX)  GO-TO-330 
GO  TO  329 

330  CONTINUE 

XCOMH=-XX (NAOJ) 

XCQLL=PPXS(NADJ)  


» 


J 


1 


L 
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- 329 

331 

83 

— -c**-*- 


111 


* SCALE! 


SCALEV 


“112 


114 


116 

C*** 


c*** 


120 
-C’^’^-SAVE 


— XCOOX=-XX  (NAOJ) +DX — - 

IF (XCJMM. LT.XCGLL. and. XCOLL.LT .XCODX) 

LiO  331  M = l,20  

PRINT  332 

CONTINUE  - - 

IF  (LEVEL. LT.O)  GO  TO  116 

SCALE  ANu  write  THE-MANEUVER-QATA 

IF(ISl.EU.l)  GO  TO  111 

PRINT  9U7  - — — 

CONTINUE 

DO  114  J = 1,NADJ  --  - 

XMAN(1,J)  = XMANdtJ) 

-uO  112  K = 2,4  - 

XMAN{i<,J)  = XMAN(!<tJ) 

continue  - ■ 

IF(ISl.EO.l)  GO  TO  114 
■print  903,  {XMAN(K,J) ,K=1,4) 

CONTINUE 

-IF  riSl.NE.D-PRiNT  -919vX00  ,YDO ,ZDO 

CONTINUE 

NEXT  CARO  IS  TAPE  SPECIFIC  FOR  ITT=1  FOR 
IFdTT.NE.l.AND.LEVEL.EQ 

IF  (LEVEL  .NE.O)  GO  TO  150 

CALCULATE  THE  M I SS -DI ST ANCE 

- CALL'  MISCOM  

PRINT  315 

PRINT  813  - - 

DO  120  K=l,3 

PRINT  314,-XM(K) ,YM{K) ,ZM(K) ,RM(K) 
CONTINUE 


GO  TO  83 


STATIONARY 
0. AND. VVVV.GE. 338.0)  CALL 


TARGET 

PRIME 


THE  MISS-DISTANCE-DATA- 


REFF'-rU^ 


AXM(ICN) 

= 

XM(1) 

A Y ri  ( I C N) 

— 

Y M ( 1 ) 

AZM (ICN) 

= 

ZM  ( 1) 

( I C w) 

— 

RM(1> 

AXJ(ICN) 

XM  (2) 

AY  3 (ICN) 

= 

YM (2) 

AZ  ;-,  (ICM) 

= 

ZM(2) 

ARE (ICN) 

= 

RM  ( 2 ) 

AXC(ICN) 

c 

XM  (3) 

AYC(ICN) 

= 

YM(3)--- 

AZC(ICN) 

ZM(3) 

— ARC (ICN) 

= 

RM(3LrFrr 

IF (VVVV. 

LT 

.338.0)  GO  TO  150 

JOEL. 


yjEL 


■/» 


yfJi) 


■ L , 


YMT(1)=YM(1)+19.5*YA 

ZNT(l)=ZM(l)+19.5-*-ZA 

XMT(2>=XM(2)-19.5 

YMT<2>=YM(2) 

ZMT  (2  ) =Z‘':  (2) 

XMT(3)=XM(3)-19.5 

YMT  (5)=YN(3) 

ZMT(3)=Zit(3) 

AXPT(ICN)=XMT(1) 

mYMT(ICN)=YMT(1) 

A2flT(ICN)=ZMT(l) 
AXBT(ICN)=XMT(2)  - 


U 


I 


FTNH  4.2  + 75075  07/16/7! 


! AYST  ( IC.\)  =YHr  (2)  

I AZliT  licrn  =ZMT  (2) 

AXCT{IC!i)=XHT  (3)  - . 

AYGT{ICN)=YMr(3» 

; Azcr  < u'.'M)  =ZMy  (3) 

' DO  121  :<-l,3 

- 121  RMT(K)^S0RT{XMT(K)**2  + Y‘1T{K)**2+2MT(K)  **2) 

; PRINl  3JJ 

i iRKT  ( ION)  =R^tT  (1)  - 

I ARDTaCN  ) =RMT  (2) 

i ARCT(  I3\)=RMr(3)  

I PRINT  ftla 

I PRINT  - 

I uO  122  R=1.3 

PRINT  314,  XMT  {<)  ,YMT(K)  ,ZMHK)  ,RMT  (K)  - 

I 122  G0NTINU2 

— call  LETH  — - 

: C***  GO  CACK  FOR  THE  NEXT  FLIGHT 

I 150  continue 

^ C***  AFTER  FLIGHTS  HAVE  BEEN  COMPLETED,  CALCULATE  MEAN  AND  STD  DEV 

PRINl  315  ' — ■ - - - 

PRINT  330 

CALL  MEAN  (AXM,AYM,A2M,ARM,NRUNS,XM{1) ,YM(1) ,ZM{1) ,RM(1)) 
PRINT  331' 

! CALL  MEAN  ( AX B , AY  3 , A Z3 , ARB ,NRU NS, X M ( 2 ) , YM ( 2 ) , ZM { 2 ) , RM ( 2) I 

—  PRINT  332  - — - - 

call  mean  {AXC,AYC,AZC,ARC,NRUNS,XM(3) ,YM(3) ,ZM(3) ,RM(3)) 

[ C***  GO  JACK  FUR  THE  NEXT  SET-OF  FLIGHT-  CONDITIONS 

IF  (VVVV. LT. 338.0)  GO  TO  197 

— PRINT  80  0 

P R I N T 3 1 6 


• PRINT  33  0 — 

IF(I99.EQ.l)  IPLANE  = 1 

— call  mean  (AXMT,AYMT,AZMT,ARMT,NRUNS,XMT(1) ,YMT{1) ,ZMT(1),RMT<1) > - 

PRINT  331 

IF  (199. EQ.l)  -IPLANE- =-2 

call  mean  (AXBT,AYBT,AZBT, ARBT,NRUNS,XMT (2),YMT(2) ,ZMT (2),RMT(2)) 

PRINT  332  - . 

IF (199. £0.1)  IPLANE  = 3 

— call  mean  (AXCT,AYCT,AZCT,ARCT,NRUNS,XMr{3),YMT(3),ZMT(3>,RMT(3)) 

AARUNS=SPK/NR 

AASOM  = SP;!G/NR 

AESUM=SPHS/NR 

AaRuNS  = SPM/NR  

PRINT  i02,ABRUNS,ABSUM,AASUM,AARUNS 

197  CONTINUE  --  

NR  UN'S  = Nlll 

NR  = Nill 

IF(IFT.NZ.l)  GO  TO  602 

DO  630  1=1,  10  

63G  PRINT  6C1 

IF (ICON. EQ.l)  GO  TO  635 

509  PAUSE - — 

IF(IQCl.FQ.l)  GO  TO  62 

IF  (laa.EQ.l)  GO  TO  -240  


0(,r<AM  MAIN 


/•♦//i*  0PT  = 1 


FTNH  4.2^75075 


07/16/7 


-633- 


1NL=I.NL- 
CO  TC  6J 
INNCL=NL 
INNCL=IN 
IFdNNCL 
138  = 0 
INL=INL- 
GO  TO  62 
IF  ( I'i'i.K 
IF ( lUCl. 
IF (133. £ 
GO  TO  63 
IOCl=i]  - 
GO  TO  66 
IF  ( ISiSS. 
GO  TC  41 

iaci=c 

IF (183. E 
199=0  - 
GO  TO  62 
GO  TC  60 
CONTINUE 
IF  (ilOGO 
GO  To  o 
PRINT  84 

FORMATd 

FORMATd 
3X,3H 
FORMATd 
FORMAT (3 
FORMATd 


NCL-5 

.GT.5)  GO  TO  241 


INNCL 


£.1)  GO  TO  627 
EC. 1)  GO  TO  626 
a.l)  GO  TO  50 
3 


NE.l)  GO  TO  50 


■.1)  GO  TO  50 


.tC.C)  GO  TO  300 


FORMAT  (2' 
FORMATd! 


FORMAT (of 
F0RMAT(6f 
FORMAT (1! 


FORMAT (1 1 


FORMAT (/- 
FORMATdl 


,10X,17H  RANDOM  NUMBER  IS,F10.4) 

(SH  =*Fo«4t3X^3H  =fF6.4^3Xf8H  “»F6#4t 

=,F6. 4) 

,37H  TARGET  IS  ASSUMED  TO  BE  A HELICOPTER)  — 

TIME  CRITICAL  DATA  COLLECTION  FAILURE.) 

,51H-IF  AUTOMATIC  MODE  SELECTION-IS  PLANNED,- SELECT-THE— 
4bH  DESIRED  STARTING  LAUNCH  SET  NUMBER  WITH  SENSE 

IQH  SWITCHES./  — - 

,52H  IF  MANUAL  MODE  SELECTION  IS  PLANNED,  SENSE  SWITCHES 

31H  WILL  BE  IGNORED  AT  THIS-PCINT./  - 

d2H  PUSH  START.) 

TEST  CASE  OUT  OF  LIMITS-//) 

,49H  IF  COMPLETE  PRINTOUT  IS  OESIRED,  SELECT  1 WITH 
16H  SENSE  SWITCHES./  - 

,51H  IF  only  essential  PRINTOUT  IS  DESIRED,  SELECT  0 

2QHWITH  SENSE  SWITCHES./  - - 

,12H  PUSH  START.) 

1.4) 

,49H  SELECT  TH£  QUALITY  CONTROL  TEST  CASE  SAMPLE  SIZE 
21H  WITH  SENSE  SWITCHES./ 

,43H  SAMPLF  Size  MUST  BE  IN  RANGE  2 THROUGH  ll.) 

,50H  IF  AUTOMATIC  MODE  SELECTION  IS  PLANNED,  AND  IT  IS 
49H-OESIREO  TO  CONTINUE  IF  A QUALITY  CONTROL  FAILURE- 
3H  OCCURS,/ 

,48H  Then  select  1 with  SENSE  SWITCHES.  OTHERWISE, 

13H  SELECT  2 .) 

1X,16H  TEST  CASE //> 

,43H  IF  THIS  IS  A MICOM  HY0RIO/IRSS  RUN,  SELECT 

23H  1 WITH-SENSe  SWITCHES./  “ • — 

,44H  IF  IT  IS  NOT,  SELECT  2 WITH  SENSE  SWITCHES./ 

,12H  PUSH  START. ) 


JG«AM  MAIN  74/7it  0PT  = 1 FTNH  4.24-75075  07/16/7 


710  f 0‘^MA  r ( 1 0 X , 4bh  itLEGT  TAPGtT  TRAJECTORY  WITH  SENSE  SWITCHES*/ 

♦ i:x,lSH  THEN  PRESS  START.) 

713  FORMAT{nx,20H  TARGET  TRAJcCTORY  »I2,15h  WAS  SELECTED./ 

* 10X,32H  IF  THIS  IS  CORRECT,  PUSH  START./ 

♦ iQX,S6H  IF  THIS  IS  WRONG,  TURN  SWITCHES  OFF 

♦ 16H  and  PUSH  START. ) 

720  FORMAT (10X,4HH  SELECT  POINT  TAU  ON  TARGET  TRAJECTORY  AT  TIME  OF 

* luH  MISSILE  LAUNCH. /lOX, 20H  USE  SENSE  SWITCHES.) 

725  FORMAT (10X,6H  TAU  = ,I2,39H  WAS  SELECTED.  TO  CONTINUE,  SELECT  1 

♦ 41H  WITH  SENSE  SWITCH  IF  TAU=0  WAS  SELECTED./ 

♦ 10X,52ri  IF  TAU=0  WAS  NOT  SELECTED,  LEAVE  SENSE  SWITCHES  AS 

• 9HTHEY  ARE./ 

♦ 1 OX, 12H  PUSH  START.  ) 

600  FORMAT  (IHO) 

801  FORMATdOX  50H  ScLECT  THE  LAUNCH  CONDITIONS  WITH  SENSE  SWITCHES./ 

* ICX  ISH  THEN  PRESS  START./ 

• lOX  50H  TO  END  RUN,  TURN  ALL  SWITCHES  ON  AND  PRESS  START  ) 

802  FORMAT  (IQX  15,  39H  IS  TOO  MANY  RUNS.  IT  IS  CHANGED  TO  20  > 

803  FORMAT (10X,16H  LAUNCH  CO ND IT 1 0 NS , 1 5 , 4X , 10 H WERE  SET./ 

* lOX  32H  IF  THIS  IS  CORRECT,  PUSH  START.  / 

♦ lOX  36H  IF  THIS  IS  WRONG,  TURN  SWITCHES  OFF 

* 16H  AND  PUSH  START.) 

805  FORMAT(li!X,13H  PSI  LAUNCH  =F7.2,4H  DEG,5X,15H  THETA  LAUNCH  =F6.2, 
•4H  D£G,5X,16H  INITIAL  RANGE  =F9.2,3H  FT) 

809  FORMAT (iaX,llH  FLIGHT  NO.  15) 

310  FORMAT  (5X  IIH  FLIGHT  NO.  I5,20X  IHX  19X  IHY  19X  IHZ  18X  5H  TIME) 
811  F0RMAT(26X  4(5X  F15.4)) 

613  FORMAT  (ICX  14ri  6X  3H  XM  12X  3H  YM  12X  3H  2M 

* 12X  3H  RM  ) 

614  FORMAT  ( 24X  3(5X  FIO . 2 ) , £15. 8) 

315-FCRMAT(1JX,23H  REFERENCED  TO  TAILPIPE)  

616  FORMAT (lOX, 17H  REFERENCED  TO  CG) 

623  FORMAT  {5X  15, 5X  20H  FOR  NL  IS  TOO  LARGE  ) 


630  FORMAT  ( IHO  lOX  2CH  Y-Z  PLANE  -----) 

331  FORMAT  (IHO  lOX  20H  X-Y  »LANE  - - - ) 

332  FORMAT  ( IHC  lOX  20H  X-Z  PLANE  -----) 


833-FORMAT  (IGX  27ri  TA9LE  OF  LAUNCH  CONDITIONS  ///)  - 

835  FORMATdOX, 5H  SET  I2,5X,11H  ALTITUDE  =F9.2,2X,3H  FT,2X, 

♦17H  CROSSING  ANGLE  =F7.2,2X,4H  DeG,2X,16H  INITIAL  RANGE  =F10.2, 
♦2X,3ri  FT) 

840  FORMAT  ( lOX  26H  END  OF  RUN,  SIGNING  OFF.  ) 

902  FORMATdD) 

--903  FORMAT  (1 10,2FiQ.O  > ■ 

904  FORMAT (3r 10.2) 

907  FORMAT  (IHO  4X  35H  TIME  XD  MANEUVER  YO  MANEUVER 

* 15H  ZD  MANEUVER) 

908  FORMATdX  FIO. 4,  3(3X  F8.2,4X)) 

909  FORMAT  (5H  X = F7.4,2X  5H  Y = F7.4,2X  5H  Z = F7.4,2X  8H  XDOT  = 

♦ F7.4,2X  3H  YOOT  = F7.4,2X  3H  ZDDT  = F7.4,8H  MISS=  E10.3> 

912  FCRMArdJX,7H  XTEO  =F10.2,3H  FT,5X,7H  YTEO  =F10.2,3H  FT,5X, 

•7h  ZTEC  =F10. 2, 3H  FT) 

913  FORMAT duX,62H  VERTICAL  TARGET  TRAJECTORY.  CROSSING  ANGLE  IS  NOT  A 
*PPLICADLE.) 

914  FORMATdOX, 53H  STATIONARY  TARGET.  CROSSING  ANGLE  IS  NOT  APPLICABLE 

*, ) 

915  FORMATdOX, 94H  PUSH  ALL  SENSE  SWITCHES  ON  FOR  AUTOMATIC  OPERATION 
•UR  ELSE  MANUAL  OPERATION  WILL  BE  SELECTED. /tQX,12H  PUSH  START.) 
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- •R16-F0RMAT(iax,l')H  TARGET  TRAJECTORY  tI2) 

517  FORMAT(10X,7H  TAJ  = .12) 

- 915  FORMAT(////10X,5H  SET  »I2) 

919  FOKMATdUX.GHXJG  = ,F  1 0 . 4 , 5 X , 5H  YDO  = » FI  0 . 4 , 5X  , 5HZD0  =,F10.4) 

996  F0RMAT(bI5)  _ - 

999  FO.RMAT(l2ilX,E11.4)) 

- 2000  FORMAT  (515)  - — . - . - - — 

2001  FORMAT  (r  10.  0) 

2002  FORMAT  (5X,6HISl  = , 15 ,5X , 7HNSAM  = , 15 , 5 X , 7HIC0N  = ,I5»5X,6HINL  = 
ltI5) 

2003  FORMAT  (5X,oHlII  = » 15 » 5 X , 7H IRSS  = »I5.5X»5HKA  = *I5,5X»6HITT  = ,I 
15  t 5X  » 7HriQGO  = »I5) 

20C4'F‘JRMAT  (5X,6mTAJ  = ,F10.5)--  ■ - 

2005  FORMAT  (5X,7HCLMS  = , F 1 0 . 5 , 5X , 7HCLS 1 = , FI  0 . 5 . 5 X , 7 HCLS2  = ,F10.5) 

- 2C0b  FORMAT  ('iX.SHQl  = , F 1 0 . 5 . 5 X , 5 HQ2  = , FI  0 . 5 » 5X  , 5HQ3  = , FlO  . 5 1 5X,  5HQ4 

1 = ,F10.5,5X,5H05  = ,F10.5 ,5X,5H06  = ,F10.5) 

_ 2007  FORMAT  (5X,5HWl  = ,F 1 0 . 5 . 5 X , 5 HK2  = ,F 1 0 . 5 » 5X » 5HW 3 = »F10 .5 f5X,5HW4  - - 

1 = ,Fl0.5.5X,5riW5  = , F 10 . 5 .5X . 5HW6  = ,F10.5) 

200  8'FORMAT- (5X, 7MNPT3  -=-.I5,5X,5HDX  - , FIO . 5 , 5X , 3HXTERM  ^ ,F tO .5) 

2009  FORMAT  (5X,7HZTE0  = . F 10 . 5 ,5X ♦ 6H ASC  = , F 1 0 . 5 , 5X , 5HRR  = ,F10.5) 

2010  FORMAT  (5X.8HNRUNS=  .15)  ' - - t - -- 

2011  FORMAT  (5X,5HNL  = .15) 

300  STOP  — - - - -- 

END 
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— 6L  uO K u«  T A CN  t - 

COMMJN/MANEjV/ITT  ,Xl)EN  (30,45)  ,Y3EM(3C,45)  ,ZOEM(30,45J  , 

* XOTG,*)  (30)  , YUTGH  (30)  , ZDTGM  (30)  ,XPRIME,YPRIME, 

• ZPR I ME,  TIMES  (5  0 ,rTRAJ(60)  , I TAU  ( 6 0 ) , 1 1 1 , NA  DJ 

DATAdTRAJd),  I = 1 , 5 9 ) / 3 » 4 , 5 , 3*  6 , 7 , 3 , 9 , 1 0 , 1 1 , 1 2 , 2*1 3 , 14 , 2 *1 5 , 16  , SERIES 

*17  , 2*  13,  12, 19, 20, 21, 22, 23, 24, 25, 21, 25, 27, 28, 29, 30, 31, 21, 32, 33, 34,  SERIES 


— *35,30,3/  ,33,39,40,41,39,42,38,43,34,44,5*45/  SERIES 

UATAdTAUd  ) ,I  = l,5  9)/5  9*0/  SERIES 

UATA( XDoM (I ,1 ) ,1  = 1 ,27) /27*0. 0/  - — -II 

JATA(Y^^M(I,1)  ,1=1,27)  727*  0.0/ 

()ATA(Zi.-'cM(:,l)  ,I  = i,27)/27*C.O/  - - - 

DAT  A (XDEM (1,2) ,1=1, 27)  727*0.  0/  JJ 

— DATA(YUEM(I,2) , I =1 , 2 7 ) /2 7* 168 . 9/ 

uArA(ZuEM (I, 2) ,1=1,27) 727*0.0/ 

DATA(XJEM(:, 3)  ,1  = 1, 27)  727*0.0/  - KK 

DATA (YOEH (I ,3) ,1  = 1,27)  727*337.  8/ 

OATAtZOEMd,  3)  , 1 = 1,27)  727*0.0/  - 

uATA(XuEM(I,4) ,1=1,27) 727*0.0/  LL 


DATA( YDEn d ,4) ,1=1,27) /27*5Q6.7/ 

DATA(Z0EMd,4),I=l,27)/2  7*0.D/ 

DATA (XDEM (I , 5), 1 = 1, 27) /O., 127. 9, 25 1.2, 355. 4, 4 66. 4, 55 0.5,  61 4. 7,  ^ A 

♦55 0.6,574.8, 56 8. 6, 63 3. 2, 58 4. 7, 5 10. ,416.9,303.7,189.4,63. 2, -65. 3, 

* *-191.5,-310.7,-413.7,-511.5,-585.8,-638.9,-658.9,-674.7,-656.1/ 

DA TA(YDEM  (1,5) , 1 = 1,27)  75 75 . 5 , 66 3. 4 , 627, i , 56 8. 2 , 438 . 8 , 3 91 . 6 , 280 . 3 , 

- *158.8,31 .6, -96. 7, -22 1.5, -3 33. 4, -443., -5 31. 5, -6 00. 9, -64 8. 5, -672, 6, 

♦-672.4 ,-547.9,-599.9,-530. 2,-441.3,-336.5,-219.5,-94.6,33. 8,161./ 

■ DATKZDEMd, 5)  ,1  = 1, 27)  727*0.  0/  - 

DATA (XDEM (1,6) ,1=1,27) 727*0.0/  MM 

DATA(YJc;M(I, 6)  ,1=1,27)727*675.6/  - — ^ 

DAT A (ZDEM (1,6) ,1=1,27) /27* 0.0/ 

DATA  (XDEM  (1, 7), 1 = 1, 27) / Q., -221. 1,-41 7.  8 , -5  68 . 5 , -65  6. 6 , -6  72 . 5 , 8 

•-6 14. 2,- 4 88. 4,-30 8. 7, -95. 1,129., 338. 9, 51 1.5, 627. 7, 674.9, 647. 7, 

*549,2,39  0.3,18  8.3,-34.4,-253.2,-444.3,-5  86.3,-68  3.9,-668.3 ,-599.1, 

*-464./ 

- DATA (Y DEM (I ,7 ) , I =1 , 27 ) /6 34 . 8 , 599. 9 , 498 . 9, 342. 9, 149 . 2 ,-60 . 9 ,-264. 3 , 

*-458.6,-564.7,-628.5,-623. 2, -5 49. 2 , -41 4 . 7 , -234 . 6 , - 28 . 7 ,180 . 4, 

*369.7,518.2 ,6  09.7^634. ,5  83  .5,4  78.3,315.3,117.7,-92.9,-29  3.3, 

*-461.4/ 

- CATA(ZDEMd, 7)  ,1  = 1, 27)  7231.1, 218. 3, 131. 6, 124. 8, 54.  3, -22.2, -96. 2,  - 

*-159. 7,-205.4,-228.3,-226. 8,-199.9,-151. ,-85. 4, -10 . 4 , 65. 7, 134.5, 

- *183.0,221.9,230.3,214.2,174.1,114.8,42.8,-33.8,-106.3,-167.9/  ' 

DATA (XDEM d ,3) , 1 = 1 , 27 ) /O . , -1 27. 9, - 25 1 . 2, -365 . 4 , -46 6, 4 , -550 . 5 , C 

* “614.7,  — 656.6,— 674.8,- 668.6,- 638.2,— 58  4'^  7,— 510.,  — 416.9,—  30  8.7, 
*-139,4,-63.2,65.3,191.5,310.7,418.7,511.5,585.8,638.9,668.9,674.7, 

*65  6.1/  - 

DATA( YDEM  d ,8) , I = 1 , 2 7 ) /o 75 .6 , 663 .4 , 6 27 . 1 , 5 68 . 2, 488 . 8, 391 . 6 , 280 . 3 , 
*158.8, ^1.6,— 96.7,— 221.6,— 338.4,— 443.,  — 531.6,  — 600.9,  — 648.5,— 672.6, 
*-672.4,-647.9,-599.9,-530.2,-441.3,-336.5,-219.5,-94.6,33.8,161./ 

uA1  A(2DEMd,3)  ,1=1,27)  727*0.0/ 

DA  FA (XDEM  (1,9)  ,1  = 1,27) /27*0.0/  NN 

— DATA(YDEM(I, 9)  ,1  = 1, 27)  727*634. 6/  

DA TA(ZQEM  1 1,9) ,1  = 1,27)  727*231.1/ 

^ArA(XDLM(I,10) ,I=l,27)/27*0.0/  - 00 

oATA(YDEM(I,10) ,1=1, 27) 727*714. 2/ 

DATACZDEMd  ,1  0 ) , I =1 , 27  ) / 27*259. 9/ 

OATA(XOEM(I, 11) ,1=1, 27)727*0.0/  . 0 

DATA (YDEM (1,11) , I =1, 27 ) /844.5 , 822. 5 , 757. 7, 65 3. 5, 515. 4, 35 0. 4, 167. 1, 


; _ .TA  OML 


7^/7i* 


opr=  1 


FTNH  U.2*7507'i 


07/16/7 


*— 2**«i1*~i2l5»*tt~'iSH,4»~553.8»~o83«87  — 773»3»~832«3»”843»»~809»8» 

^“7  — *“475*2f*304«6»*li8*2t  74«4  t2o3«  t43d*«390*3t711»df 

*79b.2/ 

CATA ( ZOlM (I , 11 ), 1=1. 27) /O., -191. 4, -3 72. 3, -5 34. 8, -6 69., -768. 4, 
♦-827.8,  — 844.1  .- 316.5.-746.4.— 637.6.— 495.9,-327.6,-142.7,49.6, 
*239.4.416.6,572.3,698. 1, 787.6. 836. 2,841. 2,302.4, 722. ,603.9,454.4, 
♦281.3/ 

nATA(XacM(I,12),I=1.27)/27*0.0/  Q 

DATA( YDLM(I.12) , I = 1 , 27 ) / 8 15 . 7 , 844 . .823.4,769.7,670.9,537.3,375.7, 
•194.5,3.  2 .-13  8.2, -36  9. 9.-5  32.3,-66  7. ,- 76 7. ,- 32 7.  1, -844 . 2 ,-81 7. 3, 

♦-74 3., -6 39. 7, -493. 1,-330. 6. -14 5. 9, 46. 4. 2 36. 3, 413. 8, 569. 9, 6 96. 3/ 
DATACZOtMd  ,12)  , I = 1 , 2 7 ) / 2 1 8 . 6 , 2 8 . , - 164  . ,-347.4,-512.8,-651.5, 

♦-756. 3,-821.8,-844.5 .-823. 2,-759.2,-655.6,-517 .9,-353.3,-170.3, 

♦21. 6. 212. 3, 39 2., 551. 3, 63 1.9, 777. 1,831. 8, 84 3. 2, 810.  7, 736-  1,623. 2, 
♦477.8/ 

LATA(XQ£H  (1 ,13) , I = 1 , 27 ) / 0 . , 159 . 9, 314 . , 456 . 8 , 5 83 . ,6 88 . 1 ,768 .4 , E 

♦820.8,843.5,835.8,797.8,730.9,637.5,521.2,385.9,236.7,78.9,-81.7, 
♦-239.4,-388.4,-523.3,-639.4,-732.2,-798.6,-336.2,-843.4,-820.2/ 

DATA(YDtM(I,13) , I =1 , 2 7 ) /844. 5 , 829. 2 , 78 3 . 9 , 7 10 . 3 , 61 0 . 9, 489. 5 , 35 0 ,4, 

♦193 -5, 39. 5,-12 0.9, -277., -4 23., -553. 8,-664. 5, -751. 1,-8 10. 6, -840. 8, 
♦-840.5,-809.3,-749.9,-662. 8, -55 1.7, -42 0.6, t2 74. 4,- 11 8. 2, 42. 3, 

♦201.2/ 

6'ATA(ZDtM(I,13),I  = l,27)/27*0.0/ 

DATACXOCM (I,14),I=l,27)/0. ,-64.3,-128.2,-191.4,-253.4,-314.,  F 

*“372. 8, -429. 4, -483. 5, -534. 8, -533. ,-627. 8,-669., -706. 3, -739.5, 

•-7  68.4,-792.8,-812,6,-82  7.  8 , -8 33 . 1 , -84 3. 5 , -844 . 1 , - 839. 8, -8 30 . 6 , 
♦-316.5,-797.8,-774.4/ 

DATA(YDbM (I, 14), 1=1, 27)/ 79 3. 5, 791. 2, 784. 3, 772. 9, 757. ,7 36. 6, 712., 

♦68  3.3,65  0.6 ,o 14. 1,574. 1,53  0. 7,484. 3,435. ,3  33.2,329.2,2  73.3,215.8, 
♦157., 97. 4, 37. 1,-2 3. 3, -83. 5, -14 3. 5, -2 02. 5, -26 0.3, -3 16. 6/ 

DATA  (ZL)£M  (1  , 14), 1=1, 27)/ 28  8. 8, 238.,  2 85. 5, 281. 3, 275. 5, 268. 1,259. 2, 

*243.7,236.6,223.5,209. , 1 93 . 2 , 1 76 . 3 , 158 . 3 , 139 . 5 , 1 19 . 8 , 99- 5, 78.5 , 
♦57.2,35.4,13.5,-8.5,-30.4,-52.2,-73.7,-94.7,-115.2/ 
DATA(XO£M(I,15>,I=1,27)/27*0.0/  C 

UATA(YD£M  (I ,15) , I = I , 27 ) / 34 4. 5 , 814. 6 , 72 7 . 2 , 5 33 , 3 , 40 7. 8 , 198 . 5 , -24 , 8 , 
♦-2  46. 4, -‘♦50. 6, -622. 9, -751,  1,-8  26. 3, -34  3. 0,-300.  I, -700. 7, -5  51. 7, 

♦—  363. 7, - 150., 74. 4, 293. 4, 491.7, 655. 3, 772. 5, 835.  1,838.6, 782, 9, 671, 8/ 

OAT A (ZO£M (I ,15) ,1=1, 27)/ 0. ,-222.6,-429.4,-605.8,-739.5,-820.8, 

♦-8 44. 1,-307. 7, -714. 2, -570. 2, -3 85. 9, -17 4. 3, 49. 6, 2 70., 471. 4, 639. 4, 

♦76 2. 1,831., 84 1.2, 791. 9, 6 86. 5, 5 3 2. 7, 341. 1,125. 5, -99., -3 16. 6, -51 1.7/ 
DATA(XO£^!<I,16),I  = 1,27)/0.,64.  3,128.2,191.4,25  3.4,314.,3  72.8,  H 

♦429. 4, 433. 5, 534. 8, 533,, 627. 8, 669. ,705.  3, 739. 5, 768.  4, 792.  8,812.6, 

♦■327.3,333.1,843.5,844.1,  33  9.  8, 830. 6 , B16. 5 , 797 . 8 , 774. 4/ 

„ATA(Yl)£M(I,16),I  = 1,27)/793.5,791.2,78  4.3,772.9,75  7.,736.6,712., 
♦533, 3, 65u. 6, 514. 1,574.  1,530.7, 484. 3, 435., 383, 2, 329, 2, 273.3, 215.8, 
♦157.,  97. 4, 37. 1,-23, 3, -8 3. 6, -14 3. 5, -2 32. 5, -260. 3, -3 16. 6/ 
LJATA(Z0£M(I,16)  , 1 = 1, 27  )/28  8. 8, 2 38.,  235. 5, 26 1,3, 2 75. 5, 263,1,259. 2, 
•248. 7, 23 6. 8, 223. 5, 209., 193. 2, 176. 3, 15 8. 3, 139. 5, 119. 8, 99. 5, 78. 5, 
♦57.2, 35.4,13.5,-3.5,-30.4,-52.2,-73.7,-94.7,-115.5/ 

DATA  (X J£ K (1, 17), 1 = 1, 27) /2 7 *0.0/  I 

LiATAlYDIMd  ,17)  , I = 1 , ?7  ) / 344. 5 , 839.  , 822 . 5 , 795.  3 , 757 . 7 , 7 10  .3  ,653 . 5 , 
♦583.3,515.4,435.7,350.4,260.4,167.1,71.6,-24.8,-120.9,-215.5, 
♦-307.2,-394.9,-477.4,-553. 8,-622.9,-683.8,-735.9,-778.3,-810.6, 
♦-832.3/ 

L'ATA  ( ZOFM  (I  ,17)  , I =1,  27  ) / O . ,-96 . 3 ,-19 1 . 4, -2  33. 9, -372 . 8 ,-456 . 8, 

♦-5  34.8  ,-u05. 3, -669., -7 23. 4,-76 8. 4, -803. 3,- 827. 8, -84 1.4, -844.1, 
♦-835.8,-816.5 ,-786.5,-746.5,-696.5,-637.6,-570.2,-495.5,-414.3, 


u->TA  one  74/7^4  0-^r  = l FTNH  4.2  + 75075  Z7/lb/7 


— ♦-327.0,-230. 7, -142.7/  

L.-\IA(X0£M(I,13),I  = l,27)/0.t-159.9,-314.,-45o.a»-58  3.,-68  8.1»  J 

♦-7o.-i, 4, -820. 3,-84  3. 5, -83  5,  3, -7  97, 8, -73  0.9, -63  7.  6 , -521. 2, -3  35. 9, 

♦-2 36. 7, -73. 9, 81. 7, 239. 4, 338. 4, 5 23. 3, 639. 4, 732, 2, 79 8. 6, 83 6. 2, 84 3. 4, 

♦320.2/  

'JATA{YDEM(I  ,13)  , I = 1 , 27  ) / 34  4. 5 , 829.  2 , 78  3 . 9 , 7 1 0 . 3 , 61  0.  9 , 48  9.  5 , 35  0 . 4 , 

- — ♦198.5,39.5,-120.9,-277. ,-423. ,-553, 8, -664. 5, -751.1  ,-810.6,-640.8, 

*-340,5,-809.8 ,-7h9. 9, -66 2. 8,-551.7,-420.6,-274.4,-118.2,42.3, 

*201.2/ 

UATA(ZOEM (I ,18) ,I=l,27)/27*Q.O/ 

i:ATA(XD£M<I,l9),I=l,27)/C.,128.3,254.1,375.1,488.3  ,593.2,6  86.3,  R 

♦766. 2, 33 1.3, 360. 6, 91 3., 927. 9, 9 25., 904. 4, 86 6. 4, 811. 9,741.8,657.5, 

— *56  0.5, 45  2. 9, 3 3 6. ,213.  7, 86.3,-41.7,-16  9.5,-294.  ,-412.9/  

CATA(YOEM (I ,19) ,1=1, 27) /872. 9, 364. 5, 339. 6,798. 6, 74 2. 3, 671. 7,588.3, 

♦493. 0,33 9. 4, 277. a, 160. 9, 40. 3, -80., -199. 2, -3 14. 7, -4 24. 1,-525. 4, 
♦-616.6,-696 .1,-76 2. 1,-813. 6, -8 49. 5, -86 9. 1,-372., -8 53. 2, -82 8., 

•-731.9/ 

DATAtZDEM  (I ,19)  ,1  = 1,27) /317.  7,314.  7,305.6,290.7,2/0.2,244.5,214.1, 

♦179.7,14 1.7,101.1,53.6,14.9,-29.1,-72.5,-114.5,-154.4,-191.2,  

♦-224.4,-253.3 ,-27 7.4, -29o. 1,-309.2,-316.3,-317.4,-312.4,-301.4, 

•-284.6/  

DATA  (XD::-*  (I  ,20  ) ,1  = 1,27  )/D.  ,-96 . 3 , - 1 9 1.  6 , - 2 84 . 9 , - 37  5.  1 , -4  61 . 2 , K 

•-542.3,-617.6,-686.3,-747.5,-800.7,-845.2,-880.6,-906.5,-922.7, 

♦ — 923.8,  — 925.,— 911.2,— 887 .6,  — 854.3,— 811.9,  — 760.7 ,-701.3,-634.3, 

*-560.5,-480.7,-395.6/  

DAT  A tYDEM (I ,20 ) , I = 1 , 27 ) / 92 3. 9 , 923. 9 , 9 0 3 . 9 , 8 34. 1 , 34 9. 8 , 80 6.  3 , 754. 2, 

* 69  3. 3, 62 o., 55 1.5, 4 71. ,335.4,295.6,202.  7,107.6,11.3,-85.1,-180.6, 
♦-274.1,-36‘*. 7,  — 451.3,-533.1,  — 609.1,— 678.6,  — 740.7,  — 794. 9, —840. 4/ 

DATA(ZDZMCI,20) ,I=1,27)/27*0.0/ - 

END 


DATA  THO 


/‘t/Y'-t 


0PT  = 1 


FTMH  '♦.2+7507S 


07/16/ 


- i^LOCK  uATA  TWO  

CCM;iON/.'^ANtUV/ITT,xa£;M  tYDEMt  33,451  , ZDEM  ( 30 , 45  ) , 

♦ K0TG3 (30) , YDTGM ( 30) , ZOTCM (30 ) , XPRIME, YPRIME, 

* ZPRIf-'EtTIMPG  (50),ITRAJ(60),ITAU{60),III,NADJ 

i;ATA(XJ- P(I  ,2i)  ,I  = l,27)/27*0-0/  S 

DA  TA(YO£M  (I, 21), 1-1,27) /sea.  9, 91)1. 7, 321. 3, 6-3  3. a,  52  5. 3, 32  6.  ,107.6, 

♦-117.1,-334.9,-533.1  ,-700.2,-826.2,-304.  ,-923.9,-899.4  ,-817.3,  

*-o3 7. 5 ,-517.4 ,-317 .,-98. 1 , 126. 5 , 34 3 . S , 54 0 . 9 , 7 0 6 . 4, 83 C . 6, 90 6, 1, 

*923.7/ 

DATA(ZD£M  (I,  21 ), 1 = 1,  27) /O., -223., -4  3 3.  , -5 1 7.  6 , - 766 . 2 ,- 86  9.  8 , 

*-9 22 . 7 ,- S 21 . 5 ,- 356 . 4 ,- 7 6 0 . 7 ,-6 1 0 .5 , -424 . 5 , -2 1 3 . 7 , 9 . 6 , 2 32 • 3 , 441 . 4 , - — 

♦624.7,771.5,373.2,923. 7,920.3,853. , 755 . 2 , 6 03. 2 ,4 16 . , 2 0 4.  4 , -19. 1/ 

DATA(XnEv(I, 22), 1=1, 27)/0. ,160.1,315.3,461.2,593.2,707.5,800.7, 1 

♦36 9, 8, 913., 92  8. 3, 916. 9, 977. 5, 8 11. 9, 7 22., 61 0.5, 48  0.  7, 336. 5,  182. 3, 

♦22, 6,-137. 7, -294. ,-441. 4, -575. 6, -692. 7, -789. ,-661. 6, -908. 6/  — 

JAT A ( YDEK  (1 ,22) , 1=  1, 27 )/92 3. 9, 915., 373. 8, 806. 3, 714. 8, 6 01. 9, 471., 

♦326. , 171 . 2,11. 3,-143.9,-304. 7, -451. 3,-584.5,-700.2 ,-794. 9,-865.8, 

♦-910.8,-92 8. 6, -918. 6, -831. 2, -a  17. 3, -72 9., -61 8. 9, -4 90. 3, -34 7., 

*-193.4/  * 

l:ATA(ZjEM  (I,22),I  = l,27)/27*0.0/ 

DATA  (XGE  (I, 23), 1 = 1, 27  )/C., -96. 3, -191, 7, -284. 9, -37  5. 1,-461. 2,  U - - 

*-542. 3, -617. 5, -686. 3, -747.  5,- 3 00. 7, -345. 2-, -83  0,6,-906. 5, -9  22.7, 

*-923. 3, -925., -911. 2, -3 87. 6, -3 54. 3, -311. 9, -760. 7, -7 01. 3, -634. 3,  

*-560.5,-480.7,-395.5/ 

DATA  lYuEM  (I  ,23),I  = l,  27)/804.5, 8 0 0 , 1 , 78  7 . 2 , 765 . 7 , 73  6 . , 698 , 3 , 653.  I', 

♦6  0 0.9,5^+2.2 ,477.6,407.  9,333.  8, 256.  , 175 . 6 , 9 3.  2 , 9.  3, -73.  7, -1 56. 4 , 

♦-237.4,-315.3,-390.9,-461,7,-527.5,-587.7,-641.5,-688.4,-727.8/  

0ATA(2DEM(I,23),I=l,27)/464.5,462. , 454 . 5 , 44 2 . 1 , 424 , 9 , 4 03 . 2 , 377. 1 , 

- *346. 9, 313. ,275. 7, 235. 5, 192. 7, 147. a, 101. 4, 53. 8, 5. 7, -42. 5, -90. 3,  — 

♦-1 3 7, 1 ,-162. 3,-225. 7,-266. 6,-304.6,-339.3,-370.4,-397.4,-420,2/ 

DATA (XDEH  (I, 24)  , 1 = 1, 27 )/C,, -22 3., -4 33. , -61 7 . 6 , - 766 . 2 ,- 86 9.  8 , V 

*-922, 7, -921. 5, -8 66. 4, -760. 7, -6 10. 5, -42 4. 3,-213.7,9.6,232.3,441.4, 

*624. 7, 771. 5, 873. 2, 923. 7, 920. 2, 853. ,755. 2, 603. 2, 416. ,204.4, -19.1/  ' 

DATA(YDEr')(I,24)  , I = 1 , 27  ) /92  8 . 9 , 9 0 1.  7 , 321 . 8 , 693 . 8 , 52  5. 3 , 326.  ,107.6, 

*-117.1  ,-334,9,-533,1  , - 70 0 . 2 , -8 26. 2 , -90 4 . ,-928.9,-899.4,-817.3,  

*-5  37.  *+,-517. 4, -3 17., -98.  1,  12  6.  5 , 34  3 . 3 , 54  0 . 9 , 7 06. 4,  830.6,90  6.1, 

*923.7/  - ■ 

DATA (ZOEM (I ,24) ,1=1, 27)/27*0,Q/ 

DATA (XDE9  (I  ,25)  ,1  = 1,27) /O. ,-191.6,-375.1,-542.3,-686.3,-800.7,  L 

*-8  3 0.6, -922. 7, -925., -88 7. 6, -81 1.9, -7 01. 3, -56 0.5 ,-3 95. 6, -21 3, 7, 

*-22.6,169.5 ,354.3,523.8,670.3,789. , 873 . 2 , 919 . 8 , 926 . 8 , 894 . , 322. 6 , 

*715.9/ 

...  oaTA  (YjE^:(I,25)  , 1 = 1, 27  )/92  8. 9, 90  8. 9, 849. 3, 754. 2, 626.,  4 71.,  295. 6, 

*10  7.6,-85.1,-274.1,-451. 3, -6 C 9 . 1 , - 74 0 . 7 , -84 0 . 4 , -90 4 . ,-92 8. 6 , 

*-913. 3,-858.7,-767.1 ,-642. 6,-490.3,-315. ,-130. ,62. 5, 252-4,431.4,  

*591.9/ 

- DATA(ZutM(I,25) ,I=1,27)/27*0.0/  

DATA(XD£t^(I,26),I  = l,27)/0,,223.,433.,617.6,766.2,369.8,922.7,  M 

*921.5,56  6.4,760.7 ,610.5,4  24.5,213. 7,-9.6,-2  32.3,-4  41.4,-624.7, 

♦-7 7 1.5, -5 73. 2, -92 3. 7, -92 0.2, -8 63., -755. 2, -60 3. 2, -4 16., -2 04. 4, 19.1/ 

DATA(YDEM{I, 26) ,1=1, 27)/928. 9, 931. 7 ,321.8,693.8,525.3,326. ,107.6,  

*-117.1,-334.9,-533.1,-700. 2,-826.2,-904. ,-928.9,-899.4,-817.3, 

4, -517, 4, -317, ,-98. 1.126. 5, 343. 3. 540. 9, 706. 4, 830. 6, 906.1, 

*928.7/ 

DATA(ZaE:3(I,26)  ,I  = l,27)/27*0.0/ 

DA TA ( XDEM (I ,27) ,1=1, 27)/ C. ,191.6,375.1,542. 3,686.3,800.7,880.6,  N 

*922.7,925., 837.6, 811-9,701.3,560.5, 395.6,213. 7,22. 6,-169.5,-354.3, 


BBeanasa 


uATA  TWO  74/7^  0PT  = 1 FTNH  i*.2<-75075  07/16/7 


♦-523. S,-670.8,-739. ,-873. 2 , -9 19. 8 , -926 . 8 , -694 . ♦ -82 2 . 6 ,-7 15 . 9/ 
DATA(YnEM(I,27) , I =1,27) 7928.9,908. 9,349. 8, 754,2,626.,471.,295.6f 
♦10 7.0  ,-8  5 .1,-2  74. 1,-451. 3, -6 09. 1,-740. 7, -84 0.4,-904., -928.  6, 

♦-9 13. 3, -358. 7, -7b 7. 1,-04 2. fa, -490. 3, -31 7., -130., 62. 5. 252. 4, 
♦431.4,591.9/ 

DATA(ZQ£M(I,27),I=l,27)/27*0.0/ 

uATA(Xl}£M{I,28)  ,1  = 1,27) / 27*0.0/  PP 

DATA(Yu£M(I,28),l=l,27)/27»928.9/ 

DA TA(ZDEM (I, 28), 1=1, 27)727*0.0/ 

DA  TA(XD£!8(1  , 29)  ,1  = 1,27  )727*0.0/  W 

UAI A(YD£M(I,29) , I = 1 , 27 ) / 92 8. 9 , 90 8. 9 , 849 . « , 754 . 2 , 62 6. ,4 71 . , 295. 6, 
*107.6,-85.1,-274. 1,-451.3,-609.1,-740.7,-840.4,-904. ,-928.6, 
♦-913.3,-85  8.7,-76  7.1,-64  2.6,-490.3,-317. ,-130 . ,62. 5, 252.4, 431.4, 
♦591.9/ 

OAT A <ZOEM (I, 29) ,1=1, 27)/ 0. ,-191.7,-375.1,-542.3,-686.3,-800.7, 
*-880.6,-922.7 ,-925. ,-887.6 ,-811.9,-701.3,-560.5,-395.6,-213.7, 
♦-22.6,159.5,354.3,523.3,670.3,739. ,873.2,919. 8, 925. 8, 894., 822. 6, 
*715.9/ 

DATA!  XOE  .9(1, 30),  1 = 1, 27) / O.  ,96.3,191.7,234.9,  375.1,  461.2,542.3,  0 

*0  17.6,68  6.3,74  7.5,80  0.  7 , 845. 2 , 8 8 0.  6 , 90 6. 5 , 922. 7 , 92 8. 8 , 92 5.  ,911.2, 
♦837. fa,  85 4. 3, 311. 9, 760. 7,  70 1.3, 6 34. 3, 56 0.5, 480. 7, 39 5. 6/ 

DATA(YO£M.  (1,30)  ,I  = 1 , 2 7 ) / 8 0 4 . 5 , 8 0 0. 1 , 78  7 . 2 , 765 . 7, 73 6.  , 6 98 .3 , 653 .1 , 

♦fa  0 0.9, 542. 2 ,47  7.6,40  7, 9, 33  3.8, 256. , 1 75 . 5 , 93. 2 , 9. 8, -73.  7, -1 56.4 , 
♦-237.4,-315. 3,-390.9,-461. 7,-527.5 ,-537. 7,-641.5,-683.4,-727.8/ 

— DATA (ZOEM (1,30) , I =1 , 2 7 ) / 46 4. 5 , 462. , 454 . 5 , 442 . 1 , 424 . 9 , 4 03 . 2 , 3 77. 1 , 
♦346. 9, 31 3., 2 75. 7, 2 35. 5, 192. 7, 147. 8, 101. 4 ,53.3,5.7,-42.5,-90.3, 
*-137.1,-182.3,-225.7,-266.6,-304,6,-339. 3,-370.4,-397.4,-420.2/ 

DATA (XuEM (1 ,31) ,1  = 1,27) /O, ,- 16 0 . 1 , - 3 15 . 3 ,-461. 2 , -5 9 3. 2 ,- 70  7.5,  X 

♦-8 00.7,-869.8,-913. ,-92 8. 8, -91 6. 9, -8 77. 5, -811. 9, -7 22., -6 10 .5, 

*-48 0.7, -336. 5, -18 2. 3, -22. 6, 137. 7, 294., 441. 4, 575. 6, 692. 7, 789., 

- *861.6,90  8.6/ 

DATA(YD£H (I ,31) , I = 1 , 27 ) /92 8 . 9 , 915. , 8 73 . 8 , 8 05 . 3, 714 . 8 , 6 01 . 9 ,471. , 

*326. ,171 .2,11. 3,-143.9,-304.7,-451. 3,-584.5,-700.2,-794.9,-865.8, 
*-910.8,-928.6,-918.6,-881.2,-817.3,-729. ,-613.9,-490.3,-347., 

♦-193.4/ 

OATA(ZDE:1{I,31),I=1,27)/27*0.0/ 

UATA(XOEM(I,32),I=1,27)/0, ,-128.3,-254.1,-375.1,-488.6,-593.2,  Y 

♦-6 36. 3, -766. 2, -83 1.3,- 8 3 0.6, -91 3., -927. 9, -925. ,-904.4,-866.4, 
♦-311.9,-741.8,-657.5,-560.5,-452.9,-336.5,-213.7,-86.8,41.7,169.5, 
♦294. ,412. 9/ 

0ATA(YD£M(I,32) , I = 1 , 2 7 ) /8 7 2. 9, 864. 5 , 83 9. 5 , 79 3 . 6 , 74 2. 3 , 67 1. 7 , 58 8. 3 , 
♦493.6,339.4,277,3 ,160.9,40 .3,-80., -199. 2, -3 14. 7,-424. 1,-525.4, 

1,-76  2.  1,-813.6,-84  9.5,-869.1,-872.  ,-85  8. 2 , 82  8.  ,-781 . 9/ 

DATA (ZOEM (I ,32) ,1=1,27) /31 7.7, 314. 7 , 30 5. 6 , 290 . 7, 27 0. 2 , 244. 5 , 214. 1, 
♦179.7,141.7,101.1 ,58.6,14.9,-29.1,-72. 5,-114.5,-154.4,-191 .2, 
♦-224.4,-253.3,-277.4,-296. 1,-309.2,-316.3,-317.4,-312.4,-301.4, 
♦-284.6/ 

nATA(XDEM(I,33),I=l,27)/27*0.0/  P 

0ATA( YUEM (1,33) , I = 1 , 2 7 ) / 10 1 3. 4 , 1 0 1 1 . 3, 1 0 05 . 2 , 995 . , 980.9, 962.7, 
♦940.7,914.9,885.4,852.3, 31 =.8, 776. , 733 . 1 , 6 87 . 3 , 6 38 . 6 , 5 87 .4 , 533. 8 , 

♦47 8. 1,420. 4, 36 1.1, 30 0.3, 238. 2, 175. 2, 111. 6, 47. 4, -16. 9, -81. 2/ 
DATA(ZD£M(I,33),I=l,27)/a. ,-64.3,-128.4,-191.9,-254.9,-316.4, 
♦-376,8,-435.7,-492.9,-548. 1 , -6 01 . 1 , -55 1 . 7 ,-6 99, 6 ,- 744 , 7, -7 86 . 8 , 

♦-8  25.8,-861.4,-89  3.5,-92  2. , -94 6. 9, -967. 9 , -985 . ,-998. 1 , -1 00 7. 2, 

♦-1012. 3, -1013. 2,-1010. 1/ 

UATA(XDEM(I,.5  4),I  = 1,27)/27*Q.0/  QQ 

3ATA(YDul-1(I,34),I  = l,27)/?7*9  78.3/ 
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- UATA<ZQ:m(I,34> ,I=l,27)/27+262.3/  - 

jAlA(X^i.”(I,35),I  = l,27)/27  + 0.0/  Z 

DAT AtYJEM (I t35 ) , I = 1 , 27 ) / 10 13 . 4 t 988 . 4 , 9 14 . 9 , 796 . 3 , 6 3 8 . 6 ,445 . 5 , 

•238.2,15.3,-208.5,-421.9 ,-614.6,-777.1,-901.4,-981.3,-1012.9, 
•-994.7,-927.0,-314.9,-652. ,-476, 6 , -267 . 8 , -45 . 8 , 178 . 5 , 394 . , 59 0. 1 , 
•757.2,88/./ 

- — UA  fA  (ZOEf!  (1,35)  ,1  = 1, 27) /O.  ,- 22  3, 4 , -4  35 . 7 , -626 . 7 , -7  86. 8 ,-90  8.  2, 

•-9 85., -1013. 2, -99 1.7, -921. 4, -8 05. 7, -65 0.4, -463. 1,-25 3., -30. 5, 

•193. 5, 40 8., 602. 4, 767. 2, 8 94. 3, 9 77. 3, 1012. 3, 99 7. 5, 93 3. 6, 82 3.  8, 

♦o73.5,4  )0,/ 

OATA{XD-7:(I,3d)  ,I  = 1,27  )/27*Q,0/  AA 

DATA  (Y OEM  (I, 36), 1 = 1, 27) /1013. 4, 1008. 8, 995., 972. 3, 940. 7, 900. 6, 

- *85 2. 3, 7 3 0.3, 7 3 3. 1,65 3, 3, 5 3 7. 4, 5 06, 2, 42 0.4, 33 0.8, 23 8. 2, 14 3. 5,  

♦4  7.4,-49 .1,-14  5. 1,-2  39.8,-332.4,-421.9,-5  0 7.0,-5  88.7,-664.5, 

- - ♦-734.3,-797.4/ 

DATA (ZDEM (I ,36) ,I=l,27)/0. ,-96.4,-191.9,-285.6,-376.8,-464.6, 

- ♦-548. 1,-626. 7, -699. 6, -760. 1,-825. 8 , -877, 9, -922. ,-957. 8, -985., 

*-1003.2,-1012. 3,-1012.2,-1002.9,-934.6,-957.3,-921.4,-877. ,-824.8, 

♦-765.1,-698.4,-625.4/  

DATA{XDEM{I,37),I=l,27)/0. , 1 28 . 4, 2 54 . 6 , 376 . 3 , 492. 9 , 6 0 1. 1 , 6 99. 6 , BB 

*73 6. 8, 86 1.4, 922. ,967. 9, 993. 1,1012. 2, 1010. 1,991. 7, 957. 3, 90 7. 5, 
♦843.1,765.1,074.7,573.5,463.1,345.2,221.8,94.7,-33.8,-161.8/ 

DATA(YDE.M{r  ,37)  , I = 1 , 27 ) / 10  13.  4 , 1 00 5 . 2 , 930 . 8 , 94 0 . 7,  885. 4,  81 5. 8, 

*733.1,633.6,533.8,420.4,300.2,175.2,47.4,-81.2,-208.4,-332.4, 

*-451.  ,-5  5 2.  3,  - 654. 5, -756. 1, -83  5.4  f-901T4-, -952. 8, -9  8 8. 8, -10  0 8. 9, 

*-1012.8,-1000.4/ 

_ „ . . DATA(ZU£MII,37) ,I=l,27)/27*0.0/  - 

DATA (XOEM  (1,38  ) ,1  = 1, 27) /27*0.0/  CC 

DATA(  Y'u£H(I,38)  ,1  = 1,27)  /lO  13.4  ,995  . ,940.7,352.3,733.1,587.4,420.4,  ’ 

♦238.2,47.4,-145.1,-352.4,-507.6,-664.5,-797.4,-901.4,-972.7, 

“ *-10  03.9,  -10  08. 6,-971. 8,-899.  8,-795-;  3, -662. -,-5 04V8, -329. 3,- 141. 8, 

♦50.7,241.4/ 

0ATA(ZDEK(I,38) ,I=l,27)/0. ,-191.9,-376.8,-548.1,-699.6,-825.8,  

♦-922. ,-985. ,-1012.3,-1002.9,-957.3,-377. ,-765.1,-625.4,-463.1, 

♦_234.i,_94. 7, 98. ,287, 2, 466. ,628. ,767. 2, 878. 7, 958. 4, 1003.4, 1012.1, — 

♦934.2/ 

- CATA(XOEM(I  ,39)  , I = 1 , 2 7 ) /27  *0 . 0 / RR- 

UATA(YD£M  (I  ,39)  , I = 1 , 27 ) / 27* 10 1 3. 4/ 

DATA(ZDEM{I,39),I=1,27)/27*0.0/  - 

OATA(XDEM(I,40),I=1,27)/0. ,-64 . 3,- 128. 4 ,-191. 9, -25 4 . 6 , -316 .4 , OD 

_ ,-376.3, -435. 7, -492. 9, -543. 1,-601. 1,-651. 7, -699. 6, -744. 7,-786. 8, 

♦-325.3,-361.4,-393.5,-922. ,-946. 9, -967. 9, -985., -99 8. 1,-1 00 7. 2, 

*-1012.3,-1013.2,-1010.1/ 

CATA (YDEM (I ,40) ,1=1, 27) /37 7. 6, 875.8,870. 5,861.7,349.4,833. 7,814.7, 

*792.3, 766.8,738.1, 706.5,672.1,643.9,595. 2,553.1,508.7,462. 3,414. , 

♦364. 1,312. 7, 260., 2 06. 3, 151. 8, 96. 6, 41. 1,-14. 6, -70. 3/ 

D AT A ( ZD£ M ( 1,4 0 ) ,1=1, 27) /50 6. 7, 505. 7,50 2.6, 497. 5, 49 0.4, 481. 4, 470.4, 

*45 7. 4, 44 2. 7, 42 6. 2, 40 7. 9, 38 8., 366. 6, 343. 6, 3 19. 3, 293. 7,266.9,239., 

*210. 2, liJC. 5, 150. 1 ,119. 1,87. 6, 55. 8, 23. 7, -3. 5, -40. 6/  

uATA(XL)EM(I,41),I  = l,27)/0.,64.3,12  8.4,191.g,2  54.6,316,4,37  6.8,  EE 

*435.7,492.9,548.1,601.1,651.7,699.6,744.7,786.8,325.8,861.4,893.5, 

♦922. ,946.9,967.9,985. ,993. 1,1007.2,1012.3,1013.2,1010.1/ 

DATA (YDEM (I, 41) , I = 1 , 27 ) /87 7. 6 , 875 . 8 , 87 0 . 5 , 861 . 7 , 34 9 . 4 , 83 3 . 7 , 814 . 7 , 
*792. 3, 76 6. 8, 738. 1,706. 5, 6 72. 1,634. 9, 595. 2, 553. 1,508. 7, 462. 3,414., 

*364. 1,312. 7, 260. ,206. 3, 151. 8, 96. 6, 41. 1 ,-14. 6, -70. 3/ — 

DA TA( ZDEM (I, 41), 1=1, 27)/506. 7, 505. 7, 502. 6, 49 7. 5, 49 0.4, 48 1.4, 47 0.4, 
*457. 4, 442. 7, 426. 2, 407. 9, 388., 366. 6, 343. 6, 319. 3, 293. 7, 266.9, 239., 
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— — eLOCK  DATA  THRtE  - — - • 

CU«,'^QN/JOMA/l£\/EL,IPTS,XXS  (50)  , XDT  :;o , YDT  GO  ♦ Z3TG0»  RLB , COSE.  SPO.RI  , 
*GAM,£DOT .THETAL.RN , 

* PPX (50 ) ,PPY (50 ),PPZ(50),TIME(50) , TMA S ( 3 0 ) . XDTGMS ( 30) 

1 - * ,YOTGMS(3'0)  ,ZDTGMS(3fl)  .XMAn1(4.5Q)  ,XMISS(7),NT 

* ,XCOMP,YCOMP,ZCOMP,rA'(A(30)  ,D£LrAR(30),VM(30)  .G.GGC 

— ^*,XOO,Yl)O,ZDO,DXG,DYG,DZG,S2,S3,S4,S5,XU^(30)  ,YDM{30),ZDM(30)  , 

*RL0<.SCALEP,F1,F2.F3.G1,G2,G3.XC.YC,ZC,S1,RRR,SR.SPL.CTL.STL, 

— - +CPL,Al,uTI,XE(3Q)  ,YE(30)  ,ZE(30)  .ZALT.NERR.CLA.NPX.NX.CLAAdO), - 

^PHO  . ARo  . AAA  ALET,rR£AL,TMA(30).XrA,YTA,ZTA. SC ALE  V, QMM ( 10 > , QM 

- •,SA.CA,V.'1X(5C)  ,VMY  (50)  ,VHZ  (50)  

* .EGGS (17 ) ,HR7EE (17) , F I 7T  ( 9 ) . FCI77 ( 9 ) 

—  COMMCJ'J/PK/XOH  (22)  . YOH(  22)  , 2QH(  22)  .3DHX  (22)  ,3DHY  (22)  ,8DHZ(2  2) ^ 

».KILL(22) ,P(4) .SPK.SPHG, SPHS.SPM 
COMMQN/MAN£Ut//lTT.X3£.M{30,45).YO£M(30,45)»ZDEM(30t45), 

* XDTGM(3C),YDTGM(30) .ZDTGM(30) , XPRIME, YPRI ME , 

— ♦ ZPRIHE.TIMES (50) tlTRAj (60) ,ITAU (60) tlll.NAOJ 

CArA(XiJEii(I,42).I  = l,27)/0.  ,-160.2,-316.4,-454.6,-6  01.1,-722.5,  FF 

■ ♦-325.8,-906.2,-967.9,-10  03.2,-1013. 2,-997. 8,-957.3,-892.7,-805.7-, 

♦-693. 4,-573.5,-434.2,-284. 0,-126.7, 33. 8,193.5,348. 3,494.4,628., 

I—  *74  5.8,34  4.9/  • ---  --- 

l;ATA(Y0£M  (I, 42), 1 = 1, 27)/ 10  13. 4, 1000. 6, 952. 7, 900. 6, 815. 8,  710. 6, 

— ■ *537.4,449.5,300.2,143.5,-16.9,-176.9,-332.4,-479.5,-614.6,-734.3, 

*-8  35,4,-915.6, 972.7,-1005.4,-1012,  3,-994.7,-951.6,-884.6,-795.3, 

*-686.  ,-559. 5/  - 

iJArA(ZG£M  (I  ,42)  ,I  = i,27)/27*0.0/ 

— - GATA(XUEM (I ,43) ,1=1, 27)/a, ,191.9, 376.8,548.1,699.5,825.8 ,922., 985. -GG 

*,1012.3, 100 2. 9, 957, 3, 877., 765. 1,625. 4,463. 1,2 84.1, 94, 7, -98., -287. 2 

+,-466. ,-623. ,-767. 2, -878. 7 ,-953.4,-1003.4,-1012.1,-984.2/  

UATA(YU£M(I,43) , I = 1 , 2 7 ) / 97 8. 8 , 96 1 . 1 , 9 0 8 . 6 , 323 . 3 , 70 8. 1 , 567, 4 ,4 06 , 1 , 

— *230.1,45. 3, -140. 2,-321.1,-490. 3,-641.9,-770.2,-870.6,-939.6,-974.5 

* ,-9  74. 2, -933. 7,-86  9. 2,-768.2,-639. 5,-437.6,-315. ,-137.  ,49,  ,233.2/ 

DATA(2D£M(I ,43) ,I=1,27)/262.3,257.5,243.5,220.6,189.7,152.,108.8, 

*61. 7, 12. 3, -37. 6, -86., -13 1.4, -172. ,-206. 4, -233. 3, -2 51, 8, -26 1.1, 

* -2  b 1 • , -2 5 1 . 5 , - 2 32 . 9, -2  0 5 . 8 , - 17 1 . 3 , — 1 3 0 • 5 , — 35 . 2 , — 36 . 7 , 1 3.  1 , 62 .5 / - 

DATA(XO£M(I,44) ,I=l,27)/27*0.0/  HH 

DATA  (YDEMd  ,44)  , 1 = 1 , 27  ) /952. 3 , 1 0 05 . 2 , 1 0 0 3.  3 , 962 . 7,  869 . 2, 73  3 , , 

*56 0.7, 3b  a. 9, 14 3, 3, -81. 4, -302. ,-5  07,  3, -63 8. 6, -8 35. 5, -94 1.4, -1000. 9, 

*-1011.2,-371. 8,-884.5,-753.7,-585.9,-389.2,-173.4,50.9,272.7, 

*481.2,655.9/ 

— DATA(ZO£M(I,44> ,I=l,27)/346.6,128.2,-96.6,-3l6.5,-520.9,-699.7,  

* — 844.1,  — 946.9,- 1003.2,— lOlC.l,— 967.  3,  — 877.,  — 743.5,— 573.4,  — 375.1, 

— — *-153.4,66.1,287.4,494.5,677.3, 826. -8,935.6,^98.4, 1012,1,976V,  891. 9, 

♦763.9/ 

eArA(XO£M(I,45)  ,I  = l,27)/27*0.0/  - - - SS 

DATA{YDEM(I,45),I=1,27)/27*1124.8/ 

CATA(ZOEM(I,45),I  = 1,27)/27*0.0/  — _ 

DATA  (GL AA (I) ,1=1,8) /O. 043,0.0435,0.0455,0.0485,0.0515,0.059, 

* - - 0.0477,  0.033/  - - — - 

DATA  (3M«(I),I=1,8)/0.0,0.3,0.6,0.3,0.9,1.05,1.5,2.2/ 
CATA(TAMA(I),I=1,28)/0.,.25,1.,1.7,2.,2.5,3.,4.,4.5,5.,6.,7., 

* 7.5, 8. ,9., 9. 5, 10. ,11. ,11. 5, 12. ,13. ,14.,l5.,l&. ,17., Id., 19., 

- - - * 20./ 

DATA (DEL  tar (I ) , 1 = 1, 28) /O., 31. 2, 445. 6, 14 18. 5, 19 78, 1, 295 8. 1,  3945.6, 
♦ 6 8 70.6,6  3 05.6,7728.1,955  8. 1,11363.  1,1226  3.  , 131 15.  , 1458  0 . , _ _ 

* 1518  8.,  15  73  3. ,16  753. , 1 72 18 . , 1766 0 . , 1 8 5 0 5 . , 1 9295. , 2 00 38 . , 

- * 20745. ,21420. ,22060. ,22665. ,23240./  
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I 

j - jAIAtVrtd) ,1=1, 28)/125., 125. »980., 1800. .1930. ,1990. tigeo. *1890. ♦ — 

I • 1840. ,1340. ,1820. ,1800. ,1780. *1630. ,1300. ,1130. ,10 70. ,960., 

♦ 9 CO. ,370. ,820. ,760. , 72 5. ,6 90. ,660. ,620., 590. ,5  60./ 

CAT A (XOH (I) ,1  = 1,22) /1 3. 5 , 0 . , 9 . , - 1 1 . 3 , -16 . 8 , -1 3 . 2 , -1 . 5 ,-4. 2 ,-16. 3 , 
— ♦ -18. 2, -2 0.3, -6. 8, -19. 4,- 19. 4, -19. 8, 3, 4, -1.6, -5. 6, -1.6, -16. 7,-1 2. 2, 

•-16.3/ 

DATA  (YDHd)  ,I=1,22)/12^0.  ,-4.  ,4.,Q.,0.,.9,0.,-,9,3*0./ 

DATA  (ZOM(I), 1=1, 22)/.!,. 2, .6,0. ,3. 8, 2. 8,  — .3,— .3, .3, .4, 5. 7,— .3, .3, 
•.3,0.,— 1.6,— 1.6,— 1.6,— l.D,— .2,— 1.9,— .1/ 
data  ( 30 HX (I), 1 = 1, 22 )/6. 5, 12. 5, 6., 6. 5, 2. 7, 2.  3, 6. 8,  3., 1.8, 1.5, 2. 5, 
•2.6, 2.2, 2. 2, .2, 1.3, 1.3, 1.3, 1.3, 3. 1,3. 6, 2./ 

DATA  ( CD  DY(I) ,l  = l,22)/i.7,2.,1.9,2. ,.3, .4, 7. 5, 9. 5, 3. 2, .3, .1,11. 9, 

♦2. 2, 2. 2,1. ,1.6, .6, 1.6, .6, 1.9, .3, 2.1/-  - 

DATA  (3DhZ(I) ,1=1,221/2. ,2.5,2.8,2.3,2.7,1.&,.3,.3*.2*2.3*2.3,.2, 
•.1,.1,1* ,.9, .6, .9, .6, 1.9, 1.4, 2.1/ 

DATA  (KILL(I),I=1,22)/4»1,0,1,0,0,1,5»0,1,5*0,1,1/ 

- DATA  (ECnSd)  ,1  = 1, 17)/-!., -.9848, -.9397, -.836, -.7071, -.5, -.342, 

♦-.  1736,0.  0,.  17 35,. 5,. 7 C71,. 866,. 939 7,. 984 8,. 996 2,1./ 

0ATA-(HR7£E(I) ,1=1,17) /5 . 3 568E- 1 0 , 6 . 43 85E-10 , 5 . 7 85 4E-1 0, 3. 7351 E-10 

♦,1.7935E-10,6.3327E-11,4.0449E-11,2.3362E-11,2.0376E-11,2. 0276-11, 
•1.3443E-11,1.5&3E-11,1.184E-11,8.7423E-12,5;2709E-12,3.4869E-12,1- 
♦3586E-12/ 

DATA(FI7T(I) ,1=1, 9)/. 18, .332, .612, .737, .833, .92, .97, .99999, 1.0000/ 

DATA  (F3I77 (I) , I = 1 , 9 ) / 0 . 0 , .1 049 1819 , . 3 14 7545 8 , . 4 1967277, .5 24590 96 , 

•.  6295091 6,.  734  42735,.  7 85  89 133, -.976  8 0098/ 

END 
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CO^!MC^;/IXTRA/NPTS,DX,xrERM,KKK»  AS,  SE  0 , SC  , A SEO  , R , SD  ALEL  ,V  V VV  ,Cl  , 
♦XLOS,YLOS,ZLOS,CO,CX,CY,CLY,CLZ,TRP,XX(50),THE 
COMMON/COMA/LEVlL  , IPTS,XXS  (5  0)  ,XDTG0,YDTG0,2IjTG0,RLB,C0SE,  spo,ri, 
*GAM,tQOT,THETAL,RN, 

♦ PPX (50) ,PPY(50) ,PPZ(50) ,TIME(5Q) , TMAS ( 30 > , XDTGMS ( 30) 

— * - ,YjrGMS(30) ,ZDTGMS(30) ,XMAN{4,50) ,XMISS(7),NT  

♦ ,XCO’-!P,YCOMP,ZCOMP,TAMA(30)  , 0 ELT  AR  { 30  ) , VM  (3 0 ) ,G,GGG 

♦,XCO,YDO,2OO,DXG,OYG,DZG,S2,S3,SA,S5,XDM{30),YDM(30),ZDM{30), 

♦RL?,<,SCALEP  ,F1,F2,F3,G1,G2,G3,XC,YC,ZC,S1,RRR,SR,SPL,CTL  ,STL, 

— — ♦CPL,A1,7TI,XE(30) ,YE(30) ,ZE(30) ,ZALT,NERR,CLA,NPX,NX,CLAA(10), 

*PH3,ARG, AAA,SCALET,TREAL,TMA (3 0 ) , X T A , Y T A , Z T A , SC ALE V , QMM ( 10 ) ,QM 

—  ♦,SA,CA,VHX(50)  ,VHY(5C)  ,VMZ(50)  - 

♦,ECCS( 17 ) , riRZEE (17) , FI 7T (9) ,FCI77( 9) 

— COMMON  /COM3/XO,YU,ZO,XDrG,YDTG,ZaTG,XDEF,YOEF,ZDEF,ZZZ 

* fXh^JE  ,YfcjiJ£,  Zl3GE 

— - COMMON/COMP/I  SI, 199, I 3 3, CL  MS (10) ,CLS1(10) ,CLS2(10)  , NSAM, Q1 , 02 , Q3 , 

•UA,a5,a6,i;i,.')2,W3,WA,N5,.'J6  ,IFT,IPLANE,ICCN,S2P1,S1P1,A2P1,  A1P1» 

*S2P2,S1P2 ,A2P2,A1P2,S2P3,S1P3, A2P3, A1P3,S2P4,S1P4, A2P4,A1P4,S2P5,- 

♦S1P5, A2P5,A1P5 , S2 P6, S 1P& , A 2P6 , A 1P6 

— COMMON/MAN£UV/ITT,XOEM (30,45) ,YOEM (30,45) ,ZOEM(30, 45) , 

* XCTGM(30) ,YDTGM(30) ,ZDTGM(30) , XPRIME,YPRIME, 

— — * ZPRIKE, TIMES (50) ,ITRAJ(60) ,ITAU(60) , III, NADJ 

C0MM0N/GTARG/FXA,FYA,FZA,FX3,FY3,FZ3,FXC,FYC,FZC,AT,VT 

CATA  RTD,OTR/57.2957795,0.01745329/-  

C***  GENERATE  THE  DISTANCE  TABLE  FOR  COLLECTING  DATA  POINTS 

—  bCALEP  = 2v;.»  (102.375) /(I. ♦GGG)  - — - 

STPN=NPTS 

OIF=  (XTERMt  (STPN-1.  )*DX) -(  (20.-*102. 375)/ (l.  + GGG) ) 

IF  (DIF.GE.O.)  GO  TO  47 

—  NAOJ=NPTS  - — 

GO  TO  43 

47  ARG=DIF/DX  - . 

NADJ=NPTS-IFIX (ARG)-l 

48  CONTINUE  — ' 

DO  20  I=l,NAOJ 

_ IN0£X=NADJ-I  + 1 

XI  = I - 1 

XX  (INDEX)  = XTERM  + (XI* OX) ""  “ 

XXS ( INDEX) =-XX (INDEX) /SC ALEP 

— — 20  CONTINUE-  

C***  WRITE  OUT  RANGE  TABLE 

IF (ISl.NE.l)  PkINT-820, -NPTS,DX,XTERM — ' 

IF (ISl.EQ.l)  GO  TO  21 

PRINT  30  0 

PRINT  700 

—  PRINT  80  0 - - • - “ 

PRINT  705 

PRINT  3 JO — 

PRINT  322  , (XX  (I) ,XXS(I) ,I=l,NAOJ) 

— --  print  3J0  - 

21  CONTiriUE 

XO  = XDM(l)  - ■ 

YJ=YUM(1 ) 

—  ZD=ZDM  (1)  — ■ 

XOX=XD 

- - YUY  = YJ  _ - 
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call  GUIj 

Xl)TGO  = \D  tg  - - — 

rDrG0=v3rG 
iL)fGO  — 2jfG 

PrJINT  8LG,XD,Y0,ZD,X0TGt  YDTG.ZOTG 

—  ARC  = (XO  + »-2  + Y0++2) 

IF(KK<.£Q.l)  GO  TO  65 

— i£0  = AVhN2  (Zl)*ARG)  - 

65  SC=AS*DTR 

— - ASEO  = SEO  ♦ RTO  

IFdTT.NF.l)  PRINT  807, AS, ASEO 

— — C^^-^-SCAlE  INITIAL  LAUNCH' CONDITIONS 

kFEET  = RI 

RI  = R/SCALtL  - ^ 

GAM=G/SCALET 

—  — XUTGQ  = XiJ  TGO/SCALEV 

ydtgo=ydtgo/scalev 

Z0rG0  = ZlirGC/SCALEV — 

TRANSFORM  MANEUVER  VECTORS  TO  LAUNCH  COORDINATES  AND  SCALE  THEM 

IFCISl.EQ.l)  GO  TO  66'  ' - 

PRINT  800 

- PRIM7  910  - ■ " 

PRINT  800 

— — PRINT  907 

66  CONTINUE 

00  70  I = 1,NT  

XD=XDM(r ) 

Z0=Z3M(I) 

' CALL  GUID  

XOTGM(I)  = XDTG 

■ YDTGH(I)  = YDTG  

ZOTGN(I)  = ZDTG 

- XOTGMSd)  = XDTGM(I)  / SCALEV 

YUTGMS(I)  = YDTGM(I)  / SCALEV 

- ZUfGMSdJ  =“ZUTGHd)~/  SCALEV 

IFdSl.'JE.l)  PRINT  908,  TM  A ( I ) , XOM  (I ) , YDM  (I ) , ZDM  (I  ) 

70  CONTINUE  

XD=XUX 

YD  = YUY  - - “ - 

ZD=ZOZ 

IF  dSl.EC.D'GO  TO-76 ^ 

PRINT  Oil 

print  duo  • 

PRINT  907 

PRINT  8D0  

00  75  1=1, NT 

print  938,  TMA(I)  .XQTGMd)  ,YOTGM(I)  ,ZDTGM(I> 

75  continue 

76  CONTINUE  — - 

C CALCULATE  INITIAL  TARGET  ANGLE  OF  ATTACK 

AAA  = lll6.8g  + 0,0Q3894*ZZZ  ' - 

ARC  = 0.i}0  003  + ZZ2 

PhO=C.00  2 3 769  2+EXP(ARG) 

QM=VVVV/AAA 

NPX  = 2 ■ “ . - 
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CALL  iNfLRP  (QM,QMM,CLAA,NX,NPX,CLA,N£RR) 

■hT  = li5i}J.C 
ST=243.0 
Xt'U=XtL(l » 

VEL=1'L(1  ) 

ZEU=ZE(1 ) -- 

AT=2.*M*SQRT  (XEU»*2  + Y£U*»2*ZEU**2)  / (PHD*  ( VVVV^*2)  ♦ST»CLA*32.174» 

AT=AT*DTR 

C CALCULATE  INITIAL  COSE 

Ci  = COStAT>  — 

Sl-SIN(AT) 

S9  = SURT(  XD*XU  + YO*yO)  - 

CQSE=  <STL*(Z0*Cl-S9*Sl)-(Cl*(ZD*Sl/Sg) ) ♦ (XD*S2+Y0*S4) ) /VVVV 

= COSE  --  - --  --  ■ ■ 

C0^C0SE»*2 

- C0SE=C05E/1. 02375  — 

C CALCULATE  INITIAL  APPARENT  PLUME  LENGTH-TO-BRE ADTH  RATIO 

xLy  = PLO<*SQRT(l.-CO) 

RL J=RL j/SCALET 

C CALCULATE  INITIAL  LOS ^ ■ - 

XL0S--R*S2 

- YLCS  = -k'»S4  - 

ZLOS=R*STL 

C- CALCULATE  HORIZON  REFERENCE  VECTOR 

ARG  = YLCS/SURT  {XL0S*‘*2  + YL0S**2) 

DET=ACOS ( ARG)  ~ 

IFtXLOS.LT.O.O)  CY=SIN(B£T) 

IF JXLOS.GE.O.C)  CY=-SIN(CET)  ' 

CX=-COS(BET) 

- C CALCULATt:  IN  IT  I AL-T  ARGE  T - CENTERLI NE  UNI  T - VECTOR 

CL Y=YD»C1/VVVV 

CLZ=-SOR T (Y0**2) ♦Sl/VVVV  — 

C CALCULATE  INITIAL  PLUME  ROTATION  ANGLE 

- ARbl=  (yLOS*CLZ-ZLOS*-CLY)-*CX*  lZLOS»CLX-XLOS*CLZ)*CY  

ARG2=S0RT  {{YL0S»CLZ-ZL0S*CLY)**2+(ZL0S*CLX-XL0S*CLZ>**2 

♦*(XLCS»CLY-YL0S*CLX)*»2) 

ARG=ARG1/ARG2 

TriE=ACOS  ( ARG)  - - - ' r— 

ARG  = XLC5*-CLY-YL0S*CLX 

- - IF(ARG.LT.O.a)  GO  TO  80 - 

IF (THE.LE.1.57079632S)  GO  TO  85 

- 1, RP  = THE-4. 71238898 — 

GO  TO  92 

aQ  TRP=1.570796326-THE  * 

GO  TO  92 

85  1RP=1.570795326*THE 
92  SP0=TRP/SCALET 

NPX=2 

CALL  INI £RP(£lll,EC05,HR7EE,NX,NPX,HR7E,NERR)  

,iC7FT=12a9.b75*HR7E*EXP(-2.341*AL0G (RFEET/22965.831) ) 
P1IRSS=(2.0*ATANIC.04884004884/RL6) ) - - 

P2IRSS=230.34375*RLB/RFEET 

AT7  = ( (P2IRSS**2)*SIN(P1IRSS) ) /(2.0* (COS (PllRSS/2.0) )**2> 

PJTJ1  = 0.  l*EXP (1.0  0 3258*ALOG(AT7) ) 

FI 7=hC7F T/PJTUl  ~ 
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NPX  = 2 

CALL  INTERP  ( F 1 7 , F 17 T , FC 17 7 , NX , NPX , FCI 7 , NERR) 

RN=FCI7 

70C  FORMAT(luX,33H  RANGE-TO-GO  TABLE  FOR  DATA  COLLECTION)  

7C5  FORMAT(22X»gH  R AN G£ , FT , 5 X , 13ri  SCALED  RANGE) 

800  FORMAT  (IHO) 

806  FORMAT(10X,3H  XOTEO  =F10.2,7H  FT/SEC, 5X,3H  YDTEO  =F10.2,7H  FT/SEC, 
»5X,3H  Z-^TEO  =F10.2,7H  FT /SEC  , /I  OX , 3 H XDTLO  =F10.2,7H  FT/SEC»5X, 

♦3H  YDTLO  =F10.2,7H  rT/SEC,5XT8H  ZDTLO  =F10.2,7H  FT/SEC) 

8C7  FORMAT (1 GX t 17h  CROSSING  ANGLE  =F7«2,4H  DEG,5X,13H  DIVE  ANGLE  =- 

♦F7.2,4H  DEG) 

820  F0RMAT(1QX,7H  NPTS  = 1 3 , 4 X 1 5H  DX -=F5 . 2 , 3H-F Ti 4X , 8H-XTERM-=F7. 2, 

*3H  FT) 

822  FOkMAT  (2CX  Fl0.2,5X  E15.8)  - ' 

907  FORMAT  ( IHO  4X  35ri  TIME  XO  MANEUVER  YO  MANEUVER 

• -15H  ZD  MANEUVER) - 

908  FCRMATdX  F10.4,  3(3X  F3.2,4X)) 

910  FORMAT (10X,24H-EARTH  FIXED  COORDINATES) 

911  FORMAT (10X,1DH  LAUNCH  COORDINATES) 

return  ' " ■ ■ ■ 

END 
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— S'J  ',Pj‘jTINE  L:.TH  — - - 

COVMOM/P-C/XDH  (22)  , Y3H(22)  , 2DH(22)  ,3QHX  ( 2 2 ) , 33HY  ( 22  ) t BDHZ  ( 2 2) 

•,<ILL  (22)  »P  (4)  ,SPK,SPh&,SPHS,SPii  

CJ’^MON/MISD/XM  (3)  ,YM  (3)  ,2M  (3)  ,RM(3)  ,M0PT,1CN 

COMMJN/MAnEUV/ITT  ,XDE!^  (3a»45)  »YD£M  ( 30,45)  ,ZO£M  (30»‘45)  » 

♦ XDTGM(3  0)  ,YDTGM  (30)  ,ZDTG!1(30)  , XPRI M E , YPRI ME  , 

—  ♦ — ZPRIME,TI«ES(50)  ,ITR4J{60)  ,ITAU{60),III,NAOJ  

C0MMeN/L'£LT/PPXS(5D)  , PP  Y S ( 50  ) , PPZS  ( 5 0 ) , VM  XS  ( 5 0 ) ,VMYS(50)  ,VMZS(50) 

- C0MMCN/GTARG/FXA,FYA,FZA ,FX3,FY3,FZ3,FXC ,FYC,FZC,AT,VT  

COMMO.“-l/3  0MO/COSA,COSa,COSG,YA,ZA 

—  COMMON/GCMP/ISi,I95,I3d,CLMS(10) ,CLS1(10) ,CLS2(10) ,NSAM,Q1,Q2,Q3,  

♦04,05,20 ,W1,W 2, W3,W4,W5,W6 ,IFT,IPLAN£, ICON , S2P1 , SI PI , A2P1, AlPl , 

— ♦■S2P2,S1P2  ,A2P2,  A1P2,S2P3,S1P3,  A2P3  , A1P3,S2P4,S1P4,  A2P4,A1P4,S2P5, 

♦S1P5, A2P5, A1P5 ,S2Pb,SlP5,A2P6, A1P6 

COSA= (1. /ZA)/SQRT( (l./ZA)»»2+<YA/ZA)*»2+l.)  ' 

C0S3= (YA/ZA)/SQRT ( (l./2A)**2+ ( YA/ZA)»*2+1. J 

—  COSG=l./'CaRT  ( (l./ZA)**2+ (YA/ZA)»»2*1.)  

SSX=0. 0 

; SSY  = Q.C  - 

; ssz=a.Q 

: — DO  4 J = 1,NAUJ  ' - - 

SSX=SSX+VMXS( J) 

j __  SSY  = SSY4-VMYS(J) — — ■ ' ■ ■ 

! SSZ=SS2+VMZS(J) 

! V-CONTINDE  ^ 

I ADJ=FLOAT (NAOJ) 

— - z***  AVERAGE  LAUNCH  COORDINATES  OF  RELATIVE  VELOCITY  VECTORS  OVER  DATA  ' - 

C COLLECTION  RANGE 

! - C00L5=  (S3X/ACJ) - 

C00L6= (S SY/AOJ) 

, — COOL7=  (SSZ/A3J)  

; Z***  GENERALIZED  TARGET  COORDINATES  OF  AVERAGE  RELATIVE  VELOCITY  VECTOR 

i C ■ REFERtNCED  TC  ORIGIN  AT  TAILPIPE  

I CUUL5=FX  A*COOL5+FYA»COOL6  + FZA»COOL7 

j — — C'JUL&  = FXD*C00L5<-FY8*C00L6+FZe»C00L7  — - " 

' CU’JL7  = FXC*CCOL5  + FYC»COOL6  + FZC*COOL7 

! — C1  = C0S(AT)“ 

i S1=SIN(AT) 

1 C»**~TARG£T  FIXED  COORDI NA TES-OF  AVERAGE-RELATIVE  VELOCITY  VECTOR 

) C referenced  TO  ORIGIN  AT  TAILPIPE 

I _ CARL5=C’JDL7^CH-CUUL5»S1 - - - 

i CARLb=-CJUL6 

i CARL7  = CD-JL5*C1-CL'UL7*S1 — — 

] C»**  TARGET  FIXED  COORDINATES  OF  AVERAGE  INERTI,.L  MISSILE  VELOCITY 

I C VECTOR  - - ■ 

I C REFERENCED  TO  ORIGIN  AT  TAILPIPE 

VM1=  {CuJL7*VT)*C14-CUUL5*S1 

' VM2=-CUUL6 

; VM3  = CUGL5»C1-(CUUL7  + VT)^S1  ' - - 

; c^**  Generalized  target  coordinates  of  position  vector  of  last  data 

' - c POINT  --  --  - - - -- 

I C REFERENCED  TO  ORIGIN  AT  TAILPIPE 

I XAI  = FXA»PPXS(NADJ)  <-FYA^PPYS(NADJ)  ♦FZA*PPZS(NAOJ) - 

I ZAI=FXC»P?XS(NAOJ) +FYC*PPYS{ NAD J)+FZC* PPZS (NAD J) 

— YAI  = FX3»PPXS(NADJ)  <-FYa»PPYS{NADJ)  ♦FZS^PPZS  (NAOJ) : 

C***  TARoET  FIXED  COORDINATES  OF  POSITION  VECTOR  OF  LAST  DATA  POINT 
C REFERENCEj  TO  ORIGIN  AT-TAILPIPE  . - 
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~ >.ilMJ:ZA  I^Cl*XAI*bl  

YAI>'J  = -Y.'.  1 

:ai:-!0  = xai  »ci-ZAi*si 

C***  01STA.\Ca  r>FTft'£EN  LAST  DATA  POINT  AND  XF-YF  PLANE  PIERCING  POINT 

C r<:;FLR:.NSc  u)  TO  ORIGIN  AT  TAILPIPE  - 

D1  = S.SRT  ( ( XAIRa*ZR  ( 1)  / 2A)  **2*  (YAIMD-YM  ( 1)  * ( Y A/2A ) *2  M ( 1 ) )»*2 

--  *-*ZAIR0**2)  — ---  

Z***  CCNPONENTS  of  unit  VECTOR  IN  DIRECTION  FROM  XF-YF  PLANE  PIERCING 

C POINT  TO  position  OF  LAST  DATA  POINT  - 

C REFcRENCZD  to  ORIGIN  AT  TAILPIPE 

AA=  (X  AIM j + ZM  (1 ) /ZA  ) /ul  

■;3=  (YAIM  0-YM{l)  * ( YA/ZA)*ZM  (1)  ) /01 

— CC  = ZAli-!J/01  - -•  

ir  (COS  A/ ASS (COSA) .NE. AA/ABS(AA) ) COSA=-COSA 

IF (COSR/AdS (0033) .NE.S3/A3S{SB) ) C0S3=-C0S3  — 

IF  iC0SG/A3S (OCSG) . NE . CC/ ABS ( CC ) ) C0SG  = -C0SG 

C***  COSINE  OF  THEfA  

C0ST-AA*C0SA+3B*C0S3+CC*C0SG 

— c***  OISTANCE  between  PK  model  aim  POINT-AND -XF-YF  PLANE-PIERCING-POINT 

C RcFERENCEO  TO  ORIGIN  AT  TAILPIPE 

D2  = Ji*COST  

C***  target  fixed  COORDINATES  OF  PK  MODEL  AIM  POINT 

C REFERENCED  TO  ORIGIN  AT  CG  - 

XAIM=02»CCSA-ZM (1) /ZA-19.5 

—  YAiM=J2»CCS3*YM(l  J-(YA/ZA)*ZM{1) 

ZAIM=D2»CCSG 

IFdSi.EJ.l)  GO  TO  2 00  

PRINT  999,XAIMD.YAIMD,ZAIMD,D1 ,AA,B3»CC»C0ST,02 

— print  999,C0SA,C0SD,C03G  

300  CONTINUE 

—  — NELIPS  = 22 - 

VSnELL=SQRT (VM1**2+VM2**2+VM3»*2) 

- tiUX  = VMl/VSHELL  - - 

GUY=VM2/ySHELL 

SUZ  = VM3/VSrlELL  

IFdSl.EC.l)  GO  TO  301 

PRINT  999,C00L5,C00L3,C00L7,CUUL5,CUUL6.CUUL7tCl,Sl 

PRINT  99  9,CARL5.CARL6»CARL7,VMl,VM2,VN3,XAI,YAI»ZftI 

PRINT  999»XAIM,YAIM,ZAIM,VSHELL,BUX,BUY,3UZ  ; 

301  CONTINUE 

— T5TAB  = ia  J.  


MP£N=C 

—  DO  4m4  M = 1,NELIPS — 

CAY1=(XAIM-XDH(M) )/30HX(M) 

— CAY3= (YAIM-YUH(M) )/30HY(M)  — 

CAY5= (ZAIM-ZDH{M) )/B0H2{M) 

CAY2=CARL5/30HX (M)  — 

CAY4=CAXL6/BDHY (M) 

—  CA Y5=CARL7/3DrtZ{M)  — 

CAY7  = CAY  Z»*2*CAY4**2  + CAY6*»2 

CAYd  = 2.*  (CAY1*CAY2*CAY3»CAY4  + CAY5*CAY6) 

CAY9=CAYl**2*CAY3**2*CAy5**2-l. 

— TEMP1=(CAY8*0.5)/CAY7  - - 

1EMP2=CAY9/CAY7 

— — Ul  SCRI=T£MP1**2-TEMP2 

IF (Cl SCRI) 444. 445, 445 


445  TTRIAL=- (EM PI-SORT (DISCR I) 
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- IF(TrSIAL-TSTAt3)446,444,444 

44b  TSlAa=TTRIAL 

Mp£,SjrM 

44  4 CCiMlNUf 

'JO  4ua  M=l,4 
4b0  P(M)=a. 

- - IF  (MPHf-i)  461  *461 ,462 

461  P(  1) =1. 

GO  TO  600 

462  IF(KiLL(MPtN) >471,471,472 
471  P(2)=l. 

GO  TO  50 : 

472-TT>  = fbTA;3 

M=MPcN 

xs=>;aim+carl5*ttt  ' 

YS=YAIM+CARL6*TTT 

2S=ZAIt^  + CARL7+TTr  

AN=2.* IXS-XOH (M) ) / (BDHX (H) **2) 

t3N  = 2,*(YS-YC!H(M)  ) /(30HY(M)+*2)  

CN=2.*<ZS-ZDH(M) ) /(BDHZCM) **Z) 

CSANG=(-AN+3JX-BN  + Bay»CN  + aUZ)/SQRT  {AN»*2+.BN++2  + CN**2)  

IF  {CGANG-.174)4  81,43l,48  2 

481  P(3)=l.  

GO  TO  503 

432-p(4)  = l.  - 

500  CONTINUE 

- --  PRINT  121,  P(l)  ,P<2)  ,P(3)  ,P(4) - — 

GP.N  = SPi<+°  (4  ) 

SPriG=SPHG+P (3)  

SPHS=SPHS+P (2) 

SPN  = SPM+P  (1)  

101  FORMAT<1HO,7H  = , F&. 4 , 3X , 1 3H  = ,F6.4  ,3X, 16H 

• =,F6.4,3X,7H =,F6.4) 

999  FORMAT (10 (IX, Ell. 4) ) 

- RETURN  - 

END 
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SU'.<POUTIilE  uOLIM  ' 

CO.-IMON/^UMP/ISl  ,I99,I88,CLMS  (10)  , CL  S 1 ( 1 0 ) , CL3  2 ( 1 0 ) ,NSAM,  Q1  »Q2,Q3, 

*34  ,^6  ,X1  ,^2,W3,W4,W5,W&,IFT,IPLANE7IC0N,S2P1,S1P1,A2P1,  AlPl,  

♦o2->2,31P2,A2P2,A1P2,32P3,S1P3»  A2P3,A1P3,S2P4,S1P4,A2P4,A1P4,S2P5, 

*S1P5, A2?5,AlP5,S2P6,SlPo,A2P6,AlP6  — 

C calculate  C0:<TR0L  LI'IITS 

i2Pl=;-.l»CLC2  (NSAX)  - 

SlPl=.-il*CLSl  (N3AM) 

A2?l  = ai  + .sl*  CLMS  (NSAM)  

AlPl^ai- W1*CLMS (NSAM) 

i2P2-l'i2+CLC2(NSAM)  - 

SlP2=.<2*Ci.31  (Ni.A'^ ) 

A2P2  = a2+W2  + CLKS (NSAM)-  ^ 

AlP2=a2-W2*CLMS (NSAH) 

i2P3=W3<^CL32  ( MSAM  ) ' — - 

S1P3=W3*CLS1  (NSAP,) 

A2P3  = Q3+W3  + ClMS(NSAM)  - - - - 

AlP3=23-r(5  + CLMS  (''ISAM) 

^2P4=;'(4*CLS2  (NSAM)  

S1P4  = '.I4*CLS1  (NSAM  ) 

A2P4  = C4+'.^4+CLMS  (NSAM)  - - 

A1P^  = C4-^I4♦CLMS  (NSAM) 

--  - j2P5=W5+CLS2(MSAM)  - “ — 

S1P5=W5*CLS1 (NSAM ) 

A2P5  = C'5+W5*CLMS  (NSAM) 

AlP5=Q5-w5+CLMS (NSAM) 

i2Pa=li6+CLS2(NsAM)  

SlP6=Wb+CLSl (NSAM) 

A2P5  = ab+W6  + CLMS  (NSAM)  ^ 

AlP6=ib-W&+CLMS (NSAM) 

— RE  T ukr; 

lNl. 


n CT 


A 
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c**  This  Su3kjgtin-:  transforms  earth  fixed  to  launch  coordinates. 

COMMON  /ANG/SINSO,GOSSa,SINSCO ,COSSCO 

common  /CCMB/X0,YD,ZD,XDTG,YDTG,2DTG,XDEFty0EF,ZDEF,2Z2 
♦,XUDE,YD0E,2UD£ 

c*+» 

■ XJTG  = XU  ■*  COSSO  ■»  COSSCO  + Y D'+COSSO* S I NSCO  - ZD*SINSO 

TG  = -XD^SINSCO  + YD*COSSCO 

TG  = X U^SINSC*CCSSCO  + Y 0*SI NSO* SI MSCO  + ZO*COSSO  


o o n o j I o o o o 
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t.U•5^iOuTI''JL  KSCALE«R£TURNS(Z1»  - — 

C***  THIS  ROUTINE  GIVES  ThE  LSTIMATED  THE  AT  1000  FT-TD-GO  (T-PRIME) 

C***  ANJ  Inc.  V \LUE  OH  GAMMA  THAT  WILL  NOT  ALLOW  ANALOG  OVERLOAD  (GAMA).  — - 
THEY  ..LlrRHINC  THE  VARIABLE  SCALE  FACTOR  EQUATION.  THE  RETURNED 
C***  values  are  G for  T-PRIKE  and  GGG  for  GAMA. 

CO.M.MGN/CQMA/LEVEL  ,IPTS  ,XXS  (5  0)  ,XDTGO,YnTGO,ZDTGO,RLB,COSE,  SPO,RI, 

•oA'-.EjCI  , THETAL.RN,  

• PPX(5C),PPY(50) ,P?Z(50),TIME(50) , TM AS ( 3 0 ) , XO TGMS (30 ) 

• ,YDTGMS(3G) ,ZDTGHS(30),XMAN(4,5C),XMISS(7),NT 

• .XCOMP. YCOMP , ZCOMP, TAN  A (30 ) ,DELTAR{30 ) , VM (3  0)  ,GfGGC 

• ,XJO,  Y.jJ  ,ZUa,UXG,UYG,UZGtS2,S3,S4,S5,XDM(3C>,YDM(3  0),ZDH(3  0)» 


^.■xL-'StSoALEPtFl  iht?.F3*GA*G2»o3*XCtYCtZC*ol  fRRRtSR^^PL^CTLfSTLf 

♦CPl , M ,V  1 I,X£ ( 30)  , YE (30)  ,ZE(30 ) , Z AL T , N ERR , CL A , NPX , NX , CLA A ( 1 0 ) , 

*r  mO ,ARG, AAA .SCAlET ,TR£AL,TNA (30) . X T A , Y T A , 2T A , SC ALE V t QMM ( 10 ) ,QM 
♦ ,S  A,CA  ,V  MX(50 ) ,VMY  (30)  ,VMZ  (50) 

•,EC3S(1Z)  ,HR7EE(17),F17T('3)  ,FCI77('3) 

CCMMCN/J3NP/IS1 , lOP, I3a,CLMS(10) ,CLS1( 10 ) ,CLS2(1C) » NS AM, Q1 , Q2, Q3 , 
♦I-,,U3,ab,!'.'l,W2,N3,W4,W5,W6,IFT,IPLANE,IC0N,S2Pl,SlPl,A2Pl,AlPl, 
'►S2P2,S1P2,A2P2,  A1PZ,S2P3,S1P3,  A2P3,  A1P3,S2P4,S1P4,  A2P4,A1P4,S2P5,  - - 
•.ilP5,  A2P5,AiP5,S2P6,SlP6  ,A2P6,A1P6 


BLGIN  CONVERGING  SEARCH  FOR  T-PRIME  AT  lOCO  FT-TO-GO. 

DG  = 1. 

AN=U. 

r.=  o. 

10  b=U*'UG 

ARC1=G»SQRT (YDTGO**2+ZOTGO**2) 

A.RG2-R  + U'^XUTGO-10  00,  

T he T A 2=A TAN2 (A RGl ,ARG2) 

NPX  = 2 - 

IF  (G.GT.2Q.)  GO  TO  20 

j _ 2 

CALL  INTcRP (G, TAMA, JELTAR, N0,NPX,0R,NERR) 

VALUE  = R-DR*C0S  (THETA2)  ♦G*XDTGO-100U. 

A=AJS ( VALUE) 

IF  (A.Lf.2a.)  GO  TO  40 
IF  (VALUE.GT.O.)  go  to  10 

AN  = ANtl.  — 

C=G-DG 

GO  TO  10 

20  P-JlNT  30  

30  FORMAT(10X,35H  T-PRiME  IS  GREATER  THAN  20  SECONDS) 

RETURN  .11  — 

40  continue 

IF(lSl.NE.l)  PRINT  50,G,UR - - 

50  FCPMATdGX,  lOH  T-PRIME  =F10.4,4H  SEC,5X,26H  DELTAR  AT  1000  FT  TO  G 
*0  =Flfl.2,3H  FT) 


CcGIN  CONVERGING  SEARCH  FOR  GAMA. 

w J - 1 • 

A.\  = U . 


) 


)UIIN£  .< SCALE 
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60  GG-Go+nL? 

— - A=-xGrGj*x'orGo 

CALL  INiE R?(GG.TAMA,UELTAR.N3,N?XtDR.NERR) 

CALL  INT::RP  {GG,TAMA,VM,NG,NPX,V»r)ERR) 

0=2.  “X  or  GC''-  (20  0C0.-R  + DR*COS{THETA2)  ) -2  0 0 0 0 . ^COS  ( THETA2 ) 

C=(R-nR*C0S(  rH£TA2  ) ) + ( 20  0 0 0 . -R+DR* C OS  ( THET  A2 ) ) 

VAlo£=A*  GG’^GG  + B'^GG  + C 

AA  = AbS(VALUE)  — — 

IF  (AA.LT.  10  000  , ) GO  TO  70 

- IF  (VALUE.GT. 0, ) GO  TO  60  

A.\  = AN  + l. 


DG=.1/AN 

GO  TO  60  ■-  

70  CONTINUE 

IF(ISl.NE.i)  PRINT  «Q,GG  

80  FORMAT (1GX,25H  T AT  GAMA  MIN  OVERLOAD  =F10.4,4H  SEC) 

-IF  lISl.NE.l)  PRINT  90, DR — - 

90  FORMAT (1  OX, 30H  DELTAR  AT  GAMA  MIN  OVERLOAD  =F10.1,3H  FT) 

IF(lSl.ME.l)  PRINT  95, V — - : 

95  FORMAT {10X,26H  VM  AT  GAMA  MIN  OVERLOAD  =F10.1,7H  FT/SEC) 

A1  = 2GGOJ./(OR»COS(TH£TA2)-R-XDTGO*GG)-^  

G&G=- (1. +A1 ) *G/GG 

- -IF  t ISI.Nl.1  ) PRINT  96, GGG 

96  FORMAT  (inx,3H  GAMA  = F10.4) 

A2  = G.43/ (l.+GGG)  

PRINT  97, A2 

97  FORMAT  (10 X,  31H-  SMALLEST-BIT-A/D  RESOLUTION^- F10.4 , 2X , 3H-FT) - 
RETURN 
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cijn'i-.UTI'.t  - — 

:0>->:0N/0OMA/LF.VE;l  ,IPTS,XXS  (5  0)  ,xdtgq,ydtgo,zdtgo,^lb,cose,  spo,ri, 
"^GAMj-DOT  ,Th£TAL»RNt 

♦ PPX (50) ,PPY (50) ,PPZ(50) ,TIME(50) , THAS ( 30 ) , XDTGMS ( 30 ) 

— * ,YDTGMS(30) ^ZDTGMS(3C) »XMAN(4,50) ,XMISS(7),NT 

♦ ,XCOMP,YCOMP,ZCOMP,TAMA(30) ,D£LTAR (30 ) tVM (30) ,G»GGG 

*,XDO,YD0,Z0O,DXr,,OYG,0ZG,S2,S3,S4,S5,XnM(30),YDM{30),ZDM(30),  

+ JK,SCALLP,Fl,F2  ,F3,G1-,G2  ,G3,XC,  YC,ZC,Sl,RRR,SRtSPL,CTL,STL, 
*CPL,A1,VTI,XE (30) ,YE(30) ,Z£(30 ) , ZALT , NERR , CL A , NPX, NX,CLA A ( 10 ) , 
»PHJ,ARG,  AAA,3CALET,TREAL  ,TMA(30)  ,XTA,YTA,ZTA,SCALEV,GMM(10  ) ,(JM 
♦,SA,CA,\/,XX(50)  ,VMY  (50)  ,(/HZ  (50)  - 

♦ ,ECOp  ( 17 ) ,HR7EE (17) , FI 7T (9) ,FCI77 (9) 

COr.MOfl/MANEUV/ITT,XJEM(3  0,45  ) , YDEM  ( 3 0 » 45  ) , ZDEM  ( 3 0 » 45 ) , 

♦ XOTGM{3C) ,YDTGM(3C) ,ZOTGM{3a) , X PRIM E,Y PRIME, 

♦ Z“RIME,  TIME5  (50  ,ITRA  J (60)  ,ITAU  (60)  , III,  NADJ  — 
COMMON  /ANG/  S I NSO , COSSO , S INSCO , CO SSCO 

COMMON/gELT/PPXS  (5  0 ,PPY3(50)  , PPZS  ( 5 0 ) , VMXS  (5  0 ) , VM  YS  (5  0)  ,VMZS(50) 

COMMON  /C0M3/XD,YD,ZD,XDTG,YDTG,ZDTG,XDEF,YDEF,ZDEF,2ZZ 

■•,XODE,YDDE,ZOD£  - - - 

COMMON  /MLV  C0SMU1,3INMU1,C0SMU2,SINMU2,C0SMU3,SINMU3 

COMMON/MISD/XM  (3) , YM (3 ) , ZM  (3) ,RM(3 ) ,NCPT , ION 

C0MM0N/CCMP/IS1,I99, I83,CLMS(10),CL31(10),CLS2(10) , NSAM, Q1 , Q2, Q3 , 

*04,05 ,06 ,W1 ,W2 ,H3 ,W4 ,W5 ,W6 ,IFT , IPL A NE , I CON , S 2P1 , SlPl , A2P 1, AlPl , 

*S2P2,SiP2 , A2P2 , A1 P2 , S2P3 , S 1P3 , A2 P3 , AlP 3 , S2 P4 , S1P4, A2P4  ,A1P4,S2P5, 

*SiP5, A2?5,A1P5,32P6,31P6,A2P6,A1P6  

REAL  HU1,MU2,MU3 

— - UATA  Rri),DTR/57.2957795,0.Q1745329/  ^ ^ 

50  CONTINUE 

TLAST  = TIMES(NAOJ)  

TFIRGT  = TIMES(I) 

TSPAN-=  TLAST--  TFIRST 

THALF  = IFIRUT  + TSPAN  / 2, 

C***  INTERP  returns  THESE  VALUES  FROM  TAJ  AOJJSTEO  TARGET- VELOCITY-TIME - 
C***  MANEUVER  TAUL'E  IN  EARTH  FIXED  COORDINATES. 

CALL  INTER?  ( TL AS T , T M A , X DM , NT , NPX , X DEL , N E RR ) - - 

CALL  INTERP  ( TF IRST , T M A , XDM , NT ,NPX , XDEF, NERR) 

XODE  = (XDEL  - XDEF)  / TSPA.N  

CALL  INTERP  (TLAST, TMA,YDM, NT, N?X,YDEL,Nt:RR) 

- - CALL  INTER?  ( TF I RS T , T M A , Y DM , NT , NPX , YDEF , NERR)  “ 

YJjE  = (YDEL  - YDEF)  / TSPAN 

--  CALL  INTERP  (TLAST, TMA,ZDM, NT, NPX,ZDEL, NERR)  

CALL  INTERP  ( TFIRST , TM A, ZDM, N T ,NPX , ZOEF, NERR) 

ZOOE=  (ZUfL-ZOEF) /TSPAN  -32.174  v - - 

C***  COMPONENTS  OF  AERODYNAMIC  ACCELERATION  IN  EARTH  FIXED  COORDINATES 


XJ 

= XDDE 

YD 

= YuDE 

ZD 

= ZDJE 

transform  these  to  launch  coordinates 

C***  components  of  AERODYNAMIC  ACCELERATION  IN  LAUNCH  COORDINATES 


XDDTG 

XDTG  — ■ 



YDDTG 

= 

YDTG 

ZDDTG 

= 

ZDTG 

r/- 

C***  DEFINE 

X- 

Y-  AND  Z0TG5 

CALL  INTERP  (THALF, TMAtXDTGM, NT, NPX,XDTG5\NERR)-  . -- 

C***  : .TERP  RETURNS  THESE  VALUES  FROM  TAJ  AOJJSTEO  TARGET  VELOCITY-TIME 
C***  MANEUVER  TABLE  IN  LAUNCH  COORDINATES.  - - 


I 


1 


Nt  P.'^IME  7h/7h  CPT  = 1 FTNH  ^.2  + 75075  07/16/7! 


— CALL  1 i'i!  C K P I I h A L r » T M A j y 0 T GM  j N T j N P X » Y D T b ^ N ERR  ) — 

CALL  lurZKP  (THALF,TMA,2'JTGM,NT,NPX,ZDTGG,NERR> 

calCulatc  the  Afioics  used  for  The  t-hatrix 

51  IF  (ZGTGt)  J)  MU1  = ATAN2  (Y0TG5,ZDTG5) 

IF  (ZjrOG.Eu.O.O.ANO.YDTGG.En.O.O)  MU1=C.0  — 

IF  ( ZOTG-G  .Ea.C.O.ANC.YQTGG.GT.O.O)  HUl  = 9a.*0TR 

— — IF  (Zl.'Tba  .t.U  . 0 . 0.  AND.  YUTG5.  LT.O  .0)  HUl  = 2 7 0 , *-DT  R 

52  IF  (XGTGv .GT.U. 0)  GO  TO  53 
IF(XuTG5.LT.O.O)  GO  TO  54 
.’■:U2  = 0. 0 

CO  TO  55 

53  IF(YOTG5.EQ.a.  J.AND.ZDTG5.EQ.C.0)  HU2=  27  0 . «-aTR 

—  ARb=SORr  (YUTG5**2  + Z0TG5<"*2)  

IF (YOTG5 .NE.n. Q .3R.20TG5 .NE.Q. 0)  M J 2 =36 0 , * OTR- A T AN  2 ( XD TG5, ARG) 

GO  TO  55  ' - — 

54  lFlYDTGb.Ea.C.G.ANn.ZOTG5.EQ.O.Q)  ,>)U2=90  . *DTR 

ARG  = S4RT  {YDTG5*»^2  + ZDTG5**-2 ) - 

IF  {YDTG5 .NE.O. O.OR.ZOTGS.NE.O. 0)  MU2=- ATAN2 (XDTG5, ARG) 

55  CONTINUE  ■ - 

DEFINE  THc  values  USED  IN  THE  T-MATRIX 

COGMUl  = COS(MUi)  — - 

SINHUI  = SIN(MUl) 

C0SMU2  = CCS(MU2)  ^ -- 

SINHU2  = SIN(MU2) 

ARG1=  ZOL!TG*SINMUi-YDDTG»COSMUi 

ARG2=XDUTG*CuSMU2+YDDTG*SINMU1*SINMU2+ZD0TG*C0SMU1»SINMU2 

- NU3  = Ar AN2 (ARGl , ARG2)  - - 

COSMU3  = C0S(MU3) 

SINRU5  = SIN(MU3)  - - - 

IFdCN.Mr.l.AND.III.Ea.l)  GO  TO  70 

—  IF(ISl.Ea.l)  GO  TO  70“  

PRINT  352,LEVEL,X0DTG,YDDTG,ZDDTGtX3TG5, YDTG5,Z0TG5 

— 70  CONTINUE  ■-  

CHANGE  radians  TO  DEGREES  FOR  PRINT  OUT 

NUl  = NUl  ♦ RTJ  

HU  2 = NU2  ♦ RTD 

—  NU3  = MU3  * RTO 

IF  (ISl.E’l.l)  GO  TO  90 

- IFdCN.EU.l.OR.III.NE.t)  PRINT  9Q1,MU1  tNU2,MU3  

90  CONTINUE 

362  FORMAT (1QX,53H  MANEUVER  NAS  IN  PROGRESS  AT  TIME  OF  DATA  COLLECTION 
♦./10X,SH  LEVEL  -=,I5,5X,8H  XDDTG  = » F 1 0. 4 , 5X  , 8H  YDDTG  =,F10.4,5X, 

— ♦AH  ZGOrc  -,F10.4/iex,8H  XDTG5-=,F10.4,5X-,8H-YDTG5'-  = ,F10,4,5X, 

*3H  ZDTG5  =,F10.4) 

901  FORMAT(1HO,10X»6H  HUl  = F8.2»5X,6H  MU2-=  F 8 . 2 , 5X » 6 H-MU3  = F8.2) 
RETURN 
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u 

C 

C 


oOi;.  ^jTINc  MISCOM  - 

THIS  £ ^k^hJJUTlNr.  CALCUL 
COOi^L'.:  .Ar£  i.Y3VEM 


CO  ..'./Ml'SP 
- CO:  -.ON/QEL 
CC'IMON  />;l 
CO.-^MON  /CD 
• ,X3Dl,Y0DE , 
COM.'lON/'1  ANE 


. (0)  ,YM  { 
V(5!])  , Y(5t\ 
0 S MJ 1 , S I 
■:o,YD,ZO 

/ITT  ,X- 


S THE  MISS  0I3TANCFS  IN  THE  TARGET 

,ZM(3) ,RM(3) ,NtICN 
Z ( 53) » VMXS<5Q) ,VMYS(5Q) ,VMZS(50)- 
1, COSMU2,SINMU2,COSHU3,SINMU3 
TG,YOTGtZDTG»XDEF, YDEF»ZDEF,ZZZ 

0,45) 


c 

ITT  ,X-  0 ,45)  , YD  EH  ( 3 0, 45) , ZDEM (3  0, 45)  , 

* XJTGH (3b)  .YuTGH  (30)  , ZDTGM (30) , X PRIME, YPRI ME , 

* ZPRIME,TIMeS (50) ,ITRAJ(6Q) ,ITAJ(60>,III,NAUJ 

COMMQN/GOMA/LE\/EL,IPTS,XXS  (50)  ,XDTGO  ,YOTGO,  Z0TG0,RL3,C0SE,  SPO,RI, 

*GAm,£OOT ,thetal,rn, 

* PPX{50),P?Y(50) ,PPZ(5Q) , TIME (50) , TMA S ( 3 0 ) , XO ' 


C 

C»** 


DDX= jX 
_ _ : J 

-XN  = FL0AT  (N) 

S,MX  = C. 

- ■:,MX2=0.  

NEXT  CARD  IS  TAPE  SPECIFIC  FOR  ITT=1  FOR-STATIONARY-TARGET 
IF(ITT.EQ.l)  GO  TO  19 

-FX  = C03MU2  ♦-COSMU3 ■ 

FY  = SINMUl  + SIf'lMUZ+C0SMU3-C0SMUl*SINMU3 

F2-SINMU1*SINMU3+C0SMU1*SINMU2*G0SMU3 

FXA  = FX 

FZA  = FZ 

- FX=CCSHU2-»SINMU3 

-FY=C0S.M;..1»C0SMU3+SINMU1»SINMU2+SINMU3 

f Z = COSMJ1  + SINMU2+SINMU3-SINMU1*COSMU3  ^ 

Fxn  = FX 

FZS  = FZ 

FY=3INMJ1+C0SMU2 

FZ  = COSMUl  ■'‘COSMU2  - 


I 


U INt  MlSCO.'t 


OPT  = l 
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FXC  =■  F K 
FYJ  = FY 
FZC  = FZ 
CONTINUE 
I F ( ( I T T . , 
IF{  (ITT.f 


Q. 1. OR. VVVV-LT. 338.0) . AND. ISl.NE. 1)  PRINT  67 
F. l.CR.VVVV.OE. 333. 0) . and. ISl .NE.l)  PRINT  66 


— 


I 


* N-cXT  CARD 
IF (ITT. E 

- IF(V\/VV. 

transform 

TARGET  CCJ 
REFERENCE 

— XT  ( I)  = 
YT(I)  = 

ZTd)  = 

CALCULATE 

AAA=1116 

VT=SORT( 

ARG=0.03 

PHO=0.00 

RM-VT/AA 

NPX'  = 2 

NX  = 3 

CALL  INT 
" NT  = 1853u 
ST=243,G 

AT=2. *WT 

IF ( lil.E 
I F ( I C M . E 
5 CUNTIi'iUE 
AT-AT'^OT 
xs-xv (I) 

— ys=yt (I) 
ZS-ZT  (I) 

TRANSFORM 
REFERENCE 
XT  1 1) =ZS 
YT (I) =-Y 

- ZT{I)=XS 
GO  TO  4 

TRANSFORM 
COORDINAT 
30  xT(I)=DX 
Y i ( I) =JX 
- ZT 1 I ) =-U 
4 CONTINUE 
IF( ISl.N 
15  F0RMAT(3 
SMX  = SMX4- 
1 SMX2=SMX 


IS  TAPE  SPECIFIC  FOR  ITT=1  FOR  STATIONARY  -TARGET 
T.l)  GO  TO  3 0 

.T.  3 3 8.0)  GO  TO  30  - - 

iMISSiLE  MINUS  TARGET)  LAUNCH  COORDINATES  TO  GENERALIZED 

‘DINATES  ■ ...  — 

) TO  ORIGIN  AT  TAILPIPE 

•XA»UX  <-  FYA^OY  <■  FZA*OZ -- 

■XB*OX  + FY3*OY  * F2B*DZ 

•XG*OX  + FYC+DY  + FZC^DZ  ^ 

ARGET  ANGLE  OF  ATTACK 

,80  + 0.  Q03894>'ZZZ  

;DCF**2+Y0EF*»2+ZDEF*»2) 

!C3*ZZZ  — - ” 

!37&92*EXP  (ARG) 


ERP  {QMtaMM,CLAA.NX,NPX,CLA,NERR) 

♦SCRT  (XDDE**2  + YODE*^^2+ZOOE**2)  / (PH0*\/T*VT*ST-»CLA*32.174) 
0.1)  GO  TO  5 

O.l.ANU.I.EQ.l)  PRINT  25, AT  ^ — 


GENERALIZED  TARGET  COORDINATES  TO  TARGET-FIXED  COORDINATES- 
U TO  ORIGIN  AT  TAILPIPE 

♦cost  AT) +XS»5IN( AT)  r 

S 

♦C0S<AT)-ZS»SIN(AT) - - 

MISSILE  launch  COORDINATES  TO  TARGET  EARTH  FIXED 

FS.  TARGET  IS  STATIONARY  OR  MOVING  AT  LESS  THAN  338  FT/SEC 

*S2-DY*SPL+DZ*S3  

*S4+uY+CPL+DZ»S5 

x-*-stl+dz*ctl 

£.1)  PRINT  15,XT(I),YT(I),ZT(I) - - 

( 3X ,E11.4)  ) 

XT(I)  

2 + XT (I ) •XT  (I ) 


DO  20  1= 
T(I)  = Y 


I 


OUTINE  MISCO:) 


/-♦/74  OPT  = l 
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20  CO  'i  i 1 Noe  — . - - - - 

CALCs^LATK  CC£FFICIt■^a3  FOR  MISSILE  PATH  PROJECTIONS  IN  TARGET  FIXED 
COORoI  -iAIES  - - 

R-Ft.RtNCE3  TO  ORIGIN  AT  TAILPIPE 


SMXI  = G . 

SMr  = SMT+  Td) 

2 SMAT^:s.!XI+Xr(U*T(I)  - 

D=XN'»SMX2-SMX*SMX 

!’=  <SMT*-SMX2-SMX^‘SKXT)/0  — 

A=  (XN»5MXT-SMT*SMX)/D 

II^  II  + 1 

99  CO  TGHia,ll),II 

YF-Axis  Intercept  and  slope  of  path  projection  in  xf-yf-plane 

R_FtRF.NC£0  TO  ORIGIN  AT  TAILPIPE 

-10  YM  ( i>  = 3 

YA  = A 

102  ^0  3 1=1, N 

— T(I)  = ZTd)  

3 CONTINUE 

GO  TO  21  

2F-AXIS  INTERCEPT  AND  SLOPE  Or  PATH  PROJECTION  IN  XF-ZF  PLANE 
REFCRtNCEO  TO  ORIGIN  AT  TAILPIPE  

11 

XF-AXIS  COORDINATE  OF  PATH-PIERCING-POINT-IN-XF-YF- PLANE 

ReFCRENCEU  TO  ORIGIN  AT  TAILPIPE 

XM  (2)=  -ZM(1)  / 2A  ■ - 

YF-AXIS  COORDINATE  OF  PATH  PIERCING  POINT  IN  XF-YF  PLANE 

Rc-FERENCtO  TO  ORIGIN  AT  TAILPIPE  — 

YM  (2) = YM(1)  + YA»XM (2 ) 


XF-AXIS  COORDINATE  OF  PATH  PIERCING  POINT  IN  XF-ZF  PLANE 

REFERENCED  TO  ORIGIN  AT  TAILPIPE  

XM  (3 )=  -YK(1)/YA 

- ZF-AXIS  COORDINATE  OF  PATH  PIERCING  POINT-IN  XF-ZF  PLANE 

referenced  TO  ORIGIN  AT  TAILPIPE 

ZM(3)-  ZM(1)-+  ZA-*XM{5)  - - --  

ZM ( 2 > = 0 . 

YM(J)=  - 

calculate  nearest  APPROACH  TO  ORIGIN  IN  YF-ZF , XF-YF, AND  XF-ZF  PLANES 

REFcRcNCED  TO  ORIGIN  AT  TAILPIPE 

DO  195  1=1,3 

155  RM(I)=  SORT (XM  (I)*XM{I)  + YH(I)*YM(I)  + ZM(I)*ZM{I))  - 

OX=ODX 

25  FORMAT (10X,25H  TARGET  ANGLE  OF  ATTACK  = ,F5.2,3X,4H  DEG) 

---66  format  {//i0X,50h  TARGET  FIXED  COORDINATES  OF  MISSILE  -REFERENCEO-TO 
• lOH  TAILPIPE.//) 

d7  FORMAT (//iax,43H  TARGET  EARTH  FIXED  COORDINATES  OF  MISSILE.//) 


3U3R0UTISt  MISCOM  74/74  OPT=l 


FTNH 


.2+75  07 


RETURN 


oooooo  ' j ol  j oooo 
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— — — SJti^OUTINE  McAN{AXM.AYMtAZH»ARM»NR^SXjSy^SZ»SS)  --  — ■ 


— THIS  SUBROUTINE  CALCULATES  THE  MEAN  ANO  STANDARD  DEi^IATlON  — 
OF  THE  MISS  DISTANCES 


COMMCN/COMP/IS1,I99,I88,CLMS{10)  .CLSKIO)  ,CLS2(1C)  ,NS AM, Q1  ,Q2 , 03 , 

— ♦04  1 05  f 06  i>si.W2'*W3tW4,W5,K6,IFT  ,IPLANE,IC0N,S2P1,  SI  P 1,  A 2P 1,  AlPl , — 
*S2P2,S1P2,A2P2,A1P2,S2P3,S1P3, A2P3, A1P3,S2P4,S1P4,A2P4,A1P4,S2P5, 

— ♦S1P5,A2P5,A1P5,S2P6,S1P6,A2P6,A1P6 

DIMENSION  AXM (2Q) , AY M ( 20 ) , AZM ( 20 ) , ARM ( 20 ) 

- XNR=FLaAr  <NR>  — • - - " - 

XNN=  SORT (XNR/ (XNR-1. ) ) 


■ 9 S3  = C.  - - - . - - 

svs=o. 

J=  J + 1 - - - 

DO  18  1=1, NR 

GO-TO  (11,12,13,14) ,J 

11  DUM=  AXM (I) 

GO  TO  15  - 

12  DUM=  AYM(I) 

-GO  TO  15  

13  0UM=  AZM (I) 

GO-TO  15  

14  DUM=  ARM  ( I)  ■■ 

15  S3=  SS+UUM  - - - - 

13  SVS=  SVS  + OUM^DUM 


MEAN 

SS=  SS/XNR 

- - SV3=  SVS/XNR  

-“STANUARO  OFVIATION - ' ■ 

—  DOD=S\/S-SS^SS 

U0J0=A35 (DOD) 

SVS  = 30Rr  (DODO)  *XNN - - 

GO  TO  (21,22,23,24) ,J 

— 21  3X=  SS - - - ' 

SVX=  SVS 

22  SY=  SS 
GC  TO  9 

—  23-SZ=  SS  - — - - 

SVZ=  SVS 

GO-TO  -9 — 

C 

i 24  PRINT  3 4 ■ 

PRINT  35,  SX,SY,SZ,SS 

PRINT  36,  SVX,SVY,SVZ,SVS  

IF  (I-jq.NE.l)  GO  TO  45 

—  IF(IPLANE.EQ.O)-GO  TO-45 

IF(IPLANE.EO.l)  GC  TO  40 
IF(I^LANE.£G.2)  GO  TO  41 


< 
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IF(IiJL'*N£.£a.3)  GO  TO  42- - - 

4 0 IF  (GV  Y.  Ly  r.S3Pl  .OR.SVY.  LT.31P1)  IFT  = 1 
IF  (GV2.GT.S2P2.0R.5\/Z.LT.S1P2)  IFT  = 1 
IF (SY.GT .A2P1. OR.SY.lt.  A IPl)  IFT  = 1 
IF(SZ.GT.A2O2.0R.SZ.LT.AlP2)  IFT=1 — 
GO  TO  45 

41-IF(S\/X.GT.S2P3.0R.SVX.LT.S1P3)  IFT-1 

IF  (SOY.GT.S2P4.0R.S\/Y.LT.S1P4)  IFT  = 1 

IF (SX.GT . A2P3.OR.SX.lt. A1P3)  IFT=1 

IF (SY.GT.A2P4. CR.SY.lt. A iP4)  IFT=1 
GO  TO  45 

42  IF  (SVX.G  r.S2P5.0R.S\/X.LT  .S1P5)  IFT  = 1 

IF  (SV2.G  T. S2P6.0R.SVZ.lt. S1P6)  IFT  = 1 

IF  (SX.GT. A2P5.0R.SX.lt. A 1P5)  IFT=1 
IF  (S2.GT.A2P6. OR.SZ.lt. A1P&) -I FT=1  ■ 

45  CONTlr(U£ 

RETURN  

C 

34-format (/ 15X2riXT 3X2 HYT6X2 HZ T8X2HRM) 

35  FORMAT(/5X4HMEAN4F10.3) 

36  FORMAT (3X6HST  0EV4F1C.4)  

ENU 
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SUriRu  J 1 
uESCaI? 
THj 
GAl 
F () 

ut; 

TH: 

UN.- 

VA;- 

INP'JT 

1 

^ - 


:Nh  lNTEPP(X,XT,YT,NX,NPX,Y,NeRR> 

ION 

, ROUTINE  USES  THE  LAGRANGE  INTERPOLATION  -ORMULA  TO  — 
:ULATE  T.HE  VALUE  Y = F(X)  FRCN  THE  FU.NCTION  TABLE  (YT(I)  = 
(I)).  THE  TABLE  OF  INOEPENjENT  VARIABLES  MUST  BE  IN 
’...ASING  ORDER  BUT  UO  NOT  HAVE  TO  3E  EVENLY  SPACED.  IF  THE 
;R£D  argument  X is  not  IN-THE  independent  VARIABLE-TABLE-*— 
ROUTINE  RILL  EXTRAPOLATE  USING  THE  NEAREST  POINTS  TO  THE 
lOhN  value.  THE  DEGREE  OF  ThE  I NT ERPOL AT  ION  FORMULA  IS  A 

:able  and  may  be  selected  by  the  user. 

THE  independent  VARIABLE  X. 

;T TABLE  OF  INDEPENDENT  X-VALUES.  MUST-BE-IN 

INCREASING  ORDER. 

T - TABLE  OF  OEPENDENT  Y VALUES.  YT ( I ) =F  ( XT  (I ) > 

IX  NUMBER  OF  POINTS  IN  XT. 

IPX  — number  of  points  USED  IN  THE  INTERPOLATING  FORMULA 
WITH  XT  AS  THE  INOEPENDE.NT  VARIABLE  TABLE.  (NPX-1> 
IS-THE'DEGREE  OF-THE-INTERPOLATION -FORMULA-USED. 


OUTPUT 


Y THE  INTERPOLATED  DEPENDENT  'VALUE.  Y = F(X> 

.NERR  NERR=0  IF  XT(1)  .LE.  X .LE.  XT{NX) 

- NERR  = 1 IF  X .LT.  XT(1)  

NERR  = 2 IF  X .GT.  XT(NX) 


DIMENSION  XT  (1)  ,YT  (1) 

* NERR  = C 

NP=NPX 

IFtNX  .LT.  NP)  — NP  =-  NX- 

NS  = (NX+loS)/26 

IH  = NP/2 

I = i 

:F(XT{I)  - X)30,20,10 
10  IH  = 0 
12  NE.RR  - 1 
GO  TO  7G 

13 -NERR  =2  — 

GO  TO  70 

20  Y = YT(H  

GO  TO  -19  9 
30  I = NX 

1F(XT(I)  - X)13,20,50 

-5C  -L  = IH  + 1 - 

IS  = NS+L 

IF (NX-IS)53 ,53,52  

52  DO  54  I=IS,NX,NS 

IF (XT (1)  - X) 54,20,56 

54  L = I 

•:,0  TO  5d  

56t.  = I-NS 

53  UO  60  I=L,NX  - 

IF(XT(I)  - X)6D,2G,70 

60  CONTINUE 
70  K = I - IH 

- N = < + NP  - 1 - — 

Y = 0.0 

IF  (N  - NX) 90,90,80 


— Js 


iia  CONTINUE 


120  Y = Y + YT{J)*P 

999  RETURN 
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SUdP'JUI  l.'iE 


analog  I NP  jT->(AOOG) 


DX ,0Y ,DZ 
X JOT 

rnoT  ■ 

zcor 

\\x 

Y Y Y 
LIL 


SCALED  XXX,yyY,Z2Z 
TIME 


X -NISSLE 
Y-MISSLE 
Z-MISSLE 
X-MISSLE 
Y-'tlSSLE 
Z-MISSLE 


MINUS  X-TARGET 
MINUS  Y-TARGET 
MINUS  Z-TARGET 
MINUS  X-TARGET 
MINUS  Y-TARGET 
MINUS  Z-TARGET 


VELOCITY( FT/SEC) 

V£LOCITY(FT/SEC)- 

VEL0CITY(FT/SEC) 

POSITION(FT) 

POSITION(FT) 

POSITION(FT) 


NOTE,  UX,UY,UZ  O'JERLOAD-UNTILL- lOCO-FT-TO  GO, HOWEVER,  XXX-YYY-Z2Z- 
SHO'JLD  NEVER  OVERLOAD 

analog  OUTPUTS (JACS) 

• X'JTGO.XCOMP  TARGET  VELOCITY  X-COMPONENT  (LAUNCH  COORD.)  

YOTGO.YCOMP  TARGET  VELOCITY  Y - COM PONEN T ( L A UNC H COORD.) 

— 2uTGO,ZCQMP  —TARGET  VELOCITY  Z-COMPONEN  T ( L A UNC  H COORD.)  

RlB  apparent  PLUME  LENGTH/BREAQT H RATIO 

^ COSE  - - - COS  OF  angle  BETWEEN  LOS  AND  CENTER  LINE  OF  TAR.^ 

5PO  SCALED  PLUME  ROTATION  ANGLE 

--  RI  INITIAL  RANGE  - - — 

gam  SCALED  T-PRIME  AT  1300  FT  TO  GO 

lDOT-  SCALED  GAMA/T-PRIME - 

ThETAL  INITIAL  ELEVATION  ANGLE  OF  TARGET  (SCALED) 

RN  - FOR  MICOM  HYBRID  - A RANDOM  NUMBER 

FOR  IR5S  APPLICATION  -IR3S  IRIS  RATIO  NO.  7 


uIuC.RETE  INPUTS 
IlN 


CDO-6600  INPUT-SENSE  LINES—  - 

SENSE  LINE  0 = C 0 0 0 00 0 0 0 0 0 0 0 00 1 

SENSE  LINE  1- = - 0 0 0 0 0 0 0 C 0 0 0 0 D 0 1 0 • 

SENSE  LINE  2 = 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 

SENSE  LINE  3 =-0000000030001000 

SENSE  LINE  4 = 00 C 0 0 0 0 0 3 0 0 1 0 0 0 0 

sense -LINE  5-=- 0000000300100000 

SENSE  LINE  6 = 0030000001000000 

SENSE  LINE  7 = C 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 

SENSE  LINE  8 = 0000000100000000 


-TR-22- 

TR-21 


uIuCRETE  OUTPUTS  (NOTE  BITS  ARE  NUMBERED  1-16  FROM  RIGHT  TO  LEFT) 

TOUT  CDC-6600  OUTPUT- CONTROL-LINES 

CNTRL  LINE  4 = 0 0 0 0 0 0 0 0 0 0 0 1 0 C 0 0 = 16  TR-33 

CNTRL  LINE  7 = 0000000010000000  = 128  TR-30 


OTHER  PROGRAM  VARIABLES  _ - . _ 

NPATri  NPATH=1,  WAITING  FOR  SENSE  LINE  6 HIGH 

— NPATH  = 2,  IN  REAL  TIME  LOOP 

IKRITE  IF  IWRITE  = 1,  WRITE  COMMENTS 

- - LEVEL  STATUS  OF  MANEUVER  - — - 

=-7  IMPLIES  NOT  IN  REAL  TIME 

__  =0  IMPLIES  IN  real  TIME 

+7  IMPLIES  TARGET  TRAJ  TABLE  EXCEEDED 

WMAN -C  IMPLIES  COLLECT  DATA,  =7  DONT  COLLECT- 

M«N  ARRAY  FOR  STORING  INTERPOLATED  VALUES 

''  MAX  MAXIMUM  POINTS  COLLECTED 


,5  O O O o O (.5  O U O O O O C3  O CJ  O o O O <3  U CJ  o o o o u o o o o o u o o o o o o 


r 


UTI\E  flight 


r,/7u 
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Ir'Tt, 


— MAXIMUM  number  OF  DATA  POINTS  THAT  CAN  BE 

ACTUALLY  COL LECTE D ( SEE  NADJ  IN  SUBROUTINE  INIT) 


:\'0 
Ir.JEX 
nC  .< 


NPTS 

IDA3 


LEN 


^ L ^ 

Slo' 

SL7 


INDEX  FOR  ASSIGNING  STORAGE  TO  XC OMP ♦ YCOM P , ZCOMP 
INDEX  OF  RANGE  TABLE  ENTRIES 
=-l,  TIME.LE. TPRIME.  =!»  TIME. GT, TPRIME 
CCJNTER,  INITIALLY  LENGTH  OF  ARRAY  TS  MINUS  ONE 

NUMBER  OF  POINTS -IN  RANGE  TABLE  - - — 

LOGICAL  VARIABLE,  =-l  IMPLIES  NOT  YET  IN  3-LOOP, 
=1  IMPLIES  FIRST  ^ASS  IN  3-LOOP  MADE.  — 

NUMBER  OF  TRUNKING  STATION  CONNECTOR  PAIRS  USED 
REAL  TIME  SIGNAL,  STATIC  CHECK  0K~ 

REAL  TIME  SIGNAL,  ICS  RECIVED 

-•REAL-TIME-SIGNAL, -RAMP  UP  READY — • - 

REAL  TIME  SIGNAL,  TERMINATION  OF  REAL  TIME 


UAl'IOS  PATCHING  REOJIREMENTS  (ONE  HOIS  AND  ONE  I3DIS) 


TRUNKING  

V- 5 u TC  W-53 
V-B2  rc  W-5U- 
V-SC  TO  w-51 
V-52  TO  W-51 


- — FORTRAN 

FOR  /IDIS2/1,IIDIS 
-FOR  /I'JIS2/l,:iDIS 
FOR  /IDIS2/2,IIDIS 
FOR  /IDIS2/2,II0IS- 


AD/4-PATCH 

TR0O-TR07  and  TR20-TR27 

TR40-TR47  AND  TR60-TR67 

TROO-TR07  AND  TR20-TR27 
TR40-TR47-AND  TR60-TR67 


V-51“TC  W-6C 
V “ 5 3 T C W - 6 0 
V-51  TO  W-61 
V-53  TO  W-ol 


F0R-/0QIS2/1,ICUIS TRi0-TR17- AND-TR30-TR37' 

FOR  /ODIS2/1,IODIS  TR50-TR57  AND  TR70-TR77 

FOR  /0UIS2/2,I0DIS  TR10-TR17  AND  TR30-TR37 

FOR  /ODIS2/2,IODIS  TR50-TR57  AND  TR70-TR77 


THE  HIGH  ORDER  CDC-560C  BIT  CORRESPONDS  TO  TROX,  WHERE  X=0,2,4,6 


DADIC3  ANALOG  PATCHING  REQUIREMENTS 


W-G3  TO  V-06 
W-t!C  TO  V-00- 


FOR  /*ADCl/49, ADC 
F0R-/*ADC1/ADC 


TR310-317  AND  TR330-337 
• TR 01 0-0 17-ANO-TR030-0 37 


-W-13-TC  V-07  - 
W-10  TO  V-Ql 
W-11  TO  V-D3 


F0R-/^0AC1/49,DAC 
FOR  /♦DACl/DAC 
FOR  /*DAC1/17,DAC 


TR35  0-3  57-AND-TR370-3  77- 
TR05Q-057  AND  TR070-077 
TR150-157  AND  TR170-177 


THE  FIRST  ADC  IN  EACH  GROUP  OF  16  CORRESPONDS  TO  AD/4-TR-10. 
THE  FIRST  DAC  IN  EACH  GROUP  OF  16  CORRESPONDS  TO  AD/4  TR-50. 


— INTEGER  PIN(B) ,P0UT(3)  - 

PC 


RE*.L  MAN  (20  3)  , MI  SSED  ( 7 ) , MISS  , L AUNG  H 
INTEGER  WMAN  — - 


XMISSd)) 

STATUS 

SPO,RI, 


o o o ^ ! o f ; o o 


luTINE  FLIGHT 


/ H / 74 
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♦C?L,Al,VT:,'<h(iG)  ,YE<30)  , Z£  ( 3 0 ) ,2AL  T , N E P.S  , CLA  , NPX  , NX  , CLA  A { 1 0 ) , 

+ PHO,As!G,  AAAtSCALtl  ,TREAL,TMA  (30)  , X T A , YT  A , ZT  A , SC  ALE  V , QMM  ( 1 0 ) , QM 

*,SA,CA,\/'IX(5.'))  ,\/MY(5l})  .\/MZ(50)  - — 

*,£C0->  (17  ),HR7EE{17),FI7T(9),FCI77(9) 

CJNH.0N/+ADC1/AQIN  (10  ) 

COHMCN/'^ jAGl/LAJ,NCH(ll> - 

COMMON/*  IL)IS2/IIN 

JAM, ON  

COMMON/* OCIS2/ICUT 

J A M , 0 N TEST 

C0M.M0N/*l)AC1/17  , A JiMTST  ( 10  ) TEST— 

jAil , OFF  TEST 

NAMELIST/TEST/LEVELtIPTS.XXS  , XDT S 0 ^ Y DTGO » ZOTG 0 i RLB » COSE, SPO# RI f 

*GAM,EDOT,TH£TAL,RN, 

* PPX,PPY,PPZ,TIME,TMAS,XDTGMS  - ■ - 

* ,YOTGMS,ZOTGMS,XKAN,XMISS,NT 

....  ,XCOMP,YCOMP,  ZCOMP,TAMA  - ,OELTAR—  ,VM  — '.G-.GGG 

*,Xj3,YC0,200,IjXG,DYG,DZG,S2,S3,S4,S5,XDM  ,YDM  ,ZDM  , 

♦RL15''C,SCALEP,F1,F2,F3,G1,G2,G3,XC,YC»ZC,S1,RRR,SR,^PL,CTL,STL, 
*CPL,A1,VTI,X£  ,YE  ,Z£  , ZALT , NERR, CLA , NPX, NX , CLA A , 

*FH0,ARG, AAA,SCAL£T,TREAL,TMA  ,xta,yta,zta,scalev,qmm  ,qm — — 

*,SA,CA,yMX  , VMY  ,VMZ 

*,ECO'J ,HR7EE ,FI7T ,FCI77- , MAN  , MISSED, -i;AUNCH,TS 

* ,IT1,KCK,ICR2,I0A3,IND, INDEX, MAX 


initialization -OF  FLAGS- ANO-COUNTERS  - 

L£y£L  = -7-  - 

IT1=1 

ICR2=29 

IDA3  = -1 

INw  = 0 

— — IN;;£X  = 0-'--  

IC jT=0 

IWRITC  --1 

MAX  = 5 0 ^ 

CHECK  PATCHING-STATUS 

PIN(1)=3RV50 

-- - PIN  (2)  =3RV51  

PIN(3)  =3P,V52 

PIN (4 ) =3RV53  

PIN(5) =3RV07 

PIN<6>=3RV06  — 

CALL  PATSTAT(PIN,POUT,LEN) 

00  10  1=1, LEN  - 

10  WRITE(6,  .3COC)PIN(I),POUT  (I) 


RESERVc  HYBRID  EOUI^MENT 


OOO  OOO  OOO  ' !oOC>  OOO  I looo'ooo 


IOijT  = <i^*lb-l  

TEST 


JAM  OUT  PR£-REAL-TIME  TEST  FUNCTIONS TEST- 

TEST 

ADiNTSTt  i)  = -1.0 

AD1NTST(2)  = l1.0 

ADINTSK  3)  = 0.  0 

AOINTST(4)=-1.0 

AOlNTST(5)=a.O 

AUINIST  ifj)  = 0.  0 

AOINTSTJ 7)=0.0 

AOIKTSK  -i)=-l.a 

ADINTSr(U)=0.0 

AOINTST ( 10) =Q.  0 
PAUSE  1 

JAM  ALL  CONTROL  LINES  LOW  IN  PREPARATION  FOR  REAL-TIME 


TEST 

TEST- 

TEST 

TEST 

TEST 

TEST- 

TEST 

TEST- 

TEST 


" -WAIT  FOR  STATIC  CHECK-COMPLETE (SENSE  LINE  5 = BIT-6)-  - 

40  CALL  INrCNVT(LARA,l) 

1-LS  = IIN.  AN0.32 

IF(SL5.NE.O)GO  TO  50 

GO  TO  4C 

- 50  -COfiTINU"  - 

WRirE(6,S000) 


PRINT  CONTENTS  OF  COMMON  BLOCK  COMA 


WRITE (6,  TEST) 


ENTER  RcAL  TIME 


call  SIMRUN (ISTAT ) 

WRITE (6,200  0 ) ISTAT 

call  REMAR;<(17H  job  in  REAL  TIME) 
CALL  OHOLC  - 


PRINT  REmAINOING  RCT  ANO  AND  REAL  TIME  -ERRORS 


r 


FLI&HV  74/74  QPT=1  FTNH  4.2+75075  07/i6/7 


- - IF  (PCTOVr.P.GT  .0  .U  ) W^IT£(6,9002)  - 

9002  FO!-lP.AT(’^  RCT  OVERUN*) 

IFtACCOVcR.GT.G.0)WRIT£(6,9003)  . 

9003  FORNAFC*  OVER  VOLTAGE  ON  ADCS+> 

IF  ( j ACO V ER . G T . u . 0 ) WRITE  (&»9004)  

9004  FORMAT{»  OVER  VOLTAGE  ON  UACS*) 

IF<AD0UV:iR.GT.G.C)WRITE(6.9D05)  

9005  FORMATf*  ADDRESS  OUT  OF  RANGE  IN  REAL  TIME*) 

KRITii  (o,  i EST)  — - — - - 

RETURN 

1 000  FORMAT(*l  RESERVATION  ERROR  CODE=*  ,2XtO20) - 

2000  F0RMAT(13H  REAL  TIME  STATuS=,O20) 

8 000*FOkMAT(/*0TRJNK— LINE  PATCHING'  l5''*tR3»*'— TO  - *^» R3 > - — 

3000  F0RMAT(22H  STATIC  CHECK  COMPLETE) 


c>  C)  c c>  o ' n a o o o o i ! i o o o 


;nE  REAL! 
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— • iOrfROUriNfc  REAL!  — 

REAL  MA,\  (c'JQ)  , M I SSED  ( 7 ) , M I S3  , L AUNC  H 
INTEGER  RKAN 
DIMENSION  T3130) 

.^QUH/ALENCE  (TS(1),  TM  AS  ( 1 ) ) , ( MAN  ( 1 ) t XM  AN  ( 1 ) ) » ( M 15  SED  ( 1)  , XMISS  ( 1 ) ) 
C0MM0N/r.l-EE£/ITl,KC<,ICR2,IDA3,IN0,  INO£XtMAX,IWRITE,NPATH,STATUS 

—  COMMON/COMA/LEVEL, IPTS ,XXS (5  0 ) ,X DT GO , Y DTGO , Z DTGO , RL3 tC OSE , SPO» R I , 

•GAM.EOOT ,TH£TAL,RN, 

♦ PPX (50 ) ,PPY (50) , PP2 (50) , TIME (50 ) , TMA S ( 3 0 ) » XO TGMS ( 30) 

• ,yDTGMS(30) ,ZOTGMS(30) ,XMAN(4,50) ,XMISS(7),NT 

* ,XCQMP,YCOMP,ZCOMP,TAMA (30)  ,D£LTAR  (30 ) tVM  (30)  ,G,GGG 
+ .XDO.YOD ,2DG,OXG,DYG,D2GtS2,S3  tS4,S5,XDM (30) , YDH{30) ,ZOM(30) , 

—  ♦RLtJ<,SC.\LIP,El,F2,F3,Gl,G2,G3,XG,YC,2C,Sl.RRRfSRoPL,CTL,STL, 

*CPL,Ai,VTI,XE (30)  ,Y£ (3  0)  ,ZE(3D ) , Z A LT , NERR , CLA , NPX , NX , CLA A ( 1 0 ) , 

♦kmO,ARG, AAA ,SCAL£T,TREAL,TMA(3  0) , X T A , Y T A . ZTA , SC ALE  V , QMM( 10 ) ,QM  

*,SA,CA,VMX(5G),VMY(50)»\/MZ(50) 

♦,ECOS(17)  ,H;R7EE(17),FI7r(9),FCI77(9) - - - - 

COMMON/^ ADC 1/ADIN ( 10 ) 

—  COMMJN/^^DACl/LAUNCHdl) 

COMMJN/+TDIS2/IIN 

— jam, ON  -■  

COMMON/^'OLISZ/IOUT 

JAM, OFF  

C 

G0‘'T0(1,2),NPATH  - 

TRANSFER  ICS  TO  AD4  - — 

LA'JNCH(i)=XOTGO  ^.  102375 

LA-JNCH(2>=YUIG0-*.  102375 

LAJNCH(J)=ZJTG0  +.102375 

LACNCH (4 ) =RLS  +.102375  — — 

LAJNCH(5)=COSe  +.102375 

LAjnCH(j)=SPO  +.102375  — 

LAUNCH{7)=RI  +.102375 

LA'JNCH(0  )=GAM  - + . 102375  

LAUNCH(P) =EDJT  +.102375 
- LAL;NCH(lG)  = THErAL+. 102375 
LAUNCH(11)=RN  +.102375 

SEND  STATUS  *3IT  TO  AO/4  INDICATING  ICS  SENKSENSE  LINE  4 = BIT  5) 


I0UT=16 

WAIT  FOR  RAM?  UP  SIGNAL  FROM  AD/4(SENSE  LINE  6 = BIT  7) 
SL6=IIN.AM0.64 

IF  (SL6.NE.fl)NPArH  = 2 

CALL  SIMIDLE 

2 CONTINUE  --  


READ 


READ  AMJ  SCALE  ADCS 


«1 


t 


bUki'i 


uTlN£  R£ALI 


7-*/7k 


0PT=1 


FTNH  2+75075 


07/16/7 


DX  = AUI^J(  1) 
jY^AJINC?) 

02  = ADIN(:<) 

— OT  = Aii  I N ( ) 

XDJT=ACIN(5) 

—  VC0T=AaiN(6) 

• ZtJ0T=A^lrU7) 

XX\'  = AJIN(8) 

YY Y=Aaru ( 9) 

ZZZ=AD1N(10) 

C 

— — C STOxt  ax,DY,OZ,XOOT,YOOT,ANL)-ZDOT  -FOR-  DIAGNOSTIC-PURPOSES 

C 

--  JXG  = -iX  + SCAL£P  

CYG=DY*3CALEP 

—  DZG  = DZ*SCALEP  

xoq=xgot 


- YD0=YUGT 

ZDO=ZDOT 


c- 

c 

COMPARE  ACTUAL  RANGE  WITH  RANGE  TABLE  ENTRY 

' c 

IF  UXS  (I.NUEX+1)  .GT.DX)  GO  TO  5 

c ■ 

PPX  (IN0EX+1)=DX 

PPY 

ppz  (i.\a.:x+i)=Dz 

TINEi INJLX+i) -DT 

VMX  ( lNJEX+1) =XDOT 

— — - 

\/HY  ( INouX  + i)  - VOOT 

VMZ  ( IN0EX+1> =ZDOT 

WMAN=7 


C 

■ -I0A3=1  

" l! 

il 

C 

TP?<Inii  - 

1 

C 

5 CONTIUUE 

j 

C 

’ *1  Kc  AL  = UT  AL c I ' 

0 o 

1 

, 

IF ( TRt AL . Lt • G ) GO  TOlO 

Ci<Ml  

KCK=1  

S<<  = GGG+  1.0 

Mioi-  UX+XlCT  + UY»YD0T+DZ'^2D0T 

IF  C-IISS.  LE.  0. 0)  GO  TO  10 

MISSc.Q(i>=DX  - 

XlbSGDCZ )=0Y 

NIjS£C(S)=u2  

'•ISSEa{'.  ) =XDOT 

M I GCE  rj  < 5 ) = Y 00 T 

MIS5tD(6)=Z00T 

:iI::>Sc.D  (7  ) =MISS  


ki 


/1.02G75 

/I, 02 37 5 — 

/I. 02375 

/I.  0 2 37  5 — - 

/I. 02375 

/I.  02375  

/1.G2375 

/1. 02375  

/i. 02375 

/1. 02375  - 


0!j  ooo  C.JOO  ooo  ooo:  pob  i ! ! !o6o  ooo 


3aoti-  cof^riNut. 


SYSTEM  hDLO(SENSE  LINE  7~=  BIT  3) 

I0UT=iE3 

CHECK  REAL  TIME  STATUS  WORD 
STATUSES  TAT (777777777777777777  778) 


RETURN  VO  BATCH  JOB 
CALL  SIMSTOP 
CKTIMc 


la  CONTIt^UE 

IF (LEVEL. CT,0)GO  TO  1000 

IF  (DT.LT.TSd)  )G0  TO  1000 

L£V£L=0 

load  7 

20  CONTINUE ; 

IF(DT.LE.TS(31-ICR2) )G0  TO  30 

INC 

IT1=ITH-1 

IF  (IT1.E0.NT)-  GO'  TO  25  

IF(ICR2.LE.1)G0  to  25 

ICR2=ICR2-1 

GO  TO  20 

■25  CONTINUE 

LEVEL=7 

-GO-TO'-IOOO  

INT INTERPOLATE  -XCOMP,  YCOMP,-ZCOMPyXCYYC  ,-ZC 

■-■30-COfiTINUE — ■ ■ - 

IT2  = IT1*-1 

IF(WMAN.EQ.O)GO  TO-100 

INJ=IN0+ 1 

IF  (MAX.LTjINO)  GO-  TO-100 T 

MAN(4-»IN0-3)=DT 

100  CONTINUE  -- 

CALC  - INTERPOLATE  XCOMP  , YCOMP  ,ZC0HP,XC , YC , ZC ) 

niV=TS(1 12) -TS (ITl) — 

RAIIO=(OT-TS(ITl) )/DIV 

XC0MP=X0TGMS(IT1) ♦RATIO* (XDTGMS{IT2)-XDTGMS(IT1)  ) 

XC=XDM  <ITl)+RATIO* (X  DM ( I T2 ) -XDM ( I T1 J J 

YC0MP=YDTGMS(1T1) ♦RATI0*(YDTGMS(IT2)-YDTGMS(IT1)) 

, I YG=YCM(IT1) ♦RATIO* (YOM (IT2) -YDM(ITl) ) 

I—'  — ZCCMP  = ZJTGMS(IT1) ♦RATIO* (ZOTGM S (IT  2 > -ZDTGMS ( ITl) > 

1 ZC=Z0M(iri)+RATI0*(ZDM(IT2)-2DM(ITl)) 

i I IF (KCK.lT.0)SKK=1.Q*GGG*TR£AL/G  


DUTINE  REALT  7U/7k  0PT=1  FTNH  h.Z*7S^75  07/16/7 


C CAT 

IFdDAC.  GT.O)GO  TO  200 

- - AAA=li.uu.5694*2ALT+1116.89  

VTI=SQRT {XC^XC+YC*YC+ZC*ZC) 

- IF(VTI.Lr.333.0)GO-  TO  2 0 0 

ARG^C. 00  0C3*2ALT 

PHO  = 0.GD237692*£XP{ARG) - 

QM  = Vr I/A  AA 

C INTERPOLATE  FOR  

NPX  = 2 

CALL  INT£RP(aM,QMM,CLAA»NX ,NPX , CLA , NERR ) 

CALL~INTCRP(TR£AL,TMA,XE,NT,NPX,XTA,N£RR) 

CALL  INTERP(TREAL,TMA,Y£,NT,NPX,YTA ,N£RR) 

-  CALL-INTEKP(TReAL-,TMA,ZE-,NT,NPX,2TA,N£RR) 

C 

- C CALCULATE  COSEtRLB^Al  ' ' - 

C 

- A1  = XTA*X  TA  + YTA-^YTA  + ZTA»ZTA  --  - — 

Al= 0.01795329*4. 6 370 3A242*SQRT|A1)/(PH0*VTI*VTI*CLA> 

SA  = SIN(A1) — 

CA=COS(A 1) 

- XXX  = XXX*2Q475. 0/SKK  - ■ 

YYY=YYy*2C475. C/SKK 

2ZZ=ZZZ*20475.0/SKK  

RRR=XXX*XXX+YYY*YYY+ZZ2*ZZZ 

: RF  L £ T = RR  R 

' RRR=VTI*3QRT {RRR> 

F1=S2*XXX-YYY*SPL  + S3*ZZZ 

F2=XXX*S4+YYY*CPL*ZZZ*S5 

F3  = GTL*ZZZ-STL*XXX 

SR=SaRT{XC*XC+YC*YC) 

S1-CA+ (SA/SR)*-2C 

G1=S1*XC 

G3=ZC*CA-SA*SR 

Elll=<Fi*Gl+F2*G2  + F3*G3)/RRR ^ 

COC£=l.Q-Elli*Elli 

!<LO=GaRT  (COSE) : 

RL6=RL3*RLBK/SCALET 

- - COS£=£lil/l. 02375 

C 

- C COMPUTE -IRSS- IRIS-RATIO- NO. -7 (PART-OF  PLUME-IRRADIANCE  COMP.)  - - 

C 

NPX  = 2 

CALL  INTERP(Elil,EC3S,HR7EE,NX ,N?X.HR7E.NERR) ^ 

HC7FT=12C9.675*HR7E*EXP(-2.341*AL0G (RFEET/22965.83 1) ) 

P1IRSS=( 2 .0*ATAn( 0. 04834004884/RLB) ) - 

P2IRSS=230, 34375*RL8/RFEET 

^77  =(  (P2IRSS**2)*SIN(P1IRSS))/(2,0*<COS(P1-IRSS/2.0))»-*2) 

pjrui=j. 1*EXP(1,0C3258*AL0G(AT7) ) 

H7  = MC7FT/PJTU1 


oooo,  I I I ! ; OOOj  ; I OOO 


r 


1 
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!MX  = 9 

MPX  = 2 

CALL  INVERP  ( F 17 , F 17 T , FC 1 77 , NX , NPX , FCI 7 , NERR)  - 

RN=FCI7 

WRlTE(6,CH£Ci<3)  - SOFT- 

COMPUTE-SPOjThE-SCALEO-PLUME- ROTATION- ANGLE 

RCl  = S£;Rr  (Fl*-Fi+F2*F2)  ■ 

RC1=F2/RC1 

RC0=AC0S  IRCl)  - - 

RC  X = ~ RC 1 

RCT  = -SIN (RCB)  

IF (Fl.GE, 0. 0)RCY=-RCY 

F11=(F2<^G2-F3»G2) /VTI— 

F22= (G1»F5-G3*F1) /VTI 

F33=tG2*r 1-G1*F2) /VTI 

FCR=SQRT (F11*F11+F22*F22+F33*F33) 

G11  = F11»RCX-  - 

G22=F22»RCY 

Clll= (G11+G22) /FCR  = 

Tlll^ACOS  (Clli) 

IF (F33.GE.Q.0) GO  TO  155 

IF (Till. L£.1.570796326)GO  TO  156 

TRP  = T111-4V71233398 — 

GO  TO  159 

155  -TRP=1. 570796326-1111 

GO  TO  159 

-155-  TRP=1.57Q79632&  + T111 

159  SP0=TRP/SCALET 

T^DO'CONTIN  UE'- 

IF  tWIlAN.FG.OlGO  TO  200G 

IF <MAX.L£.IN0) GO  TO  2000 

rAN<4*Ixr)-2)=XC0MP 

MA'Hh^IN  J-D-YCO.'IP 

■ {',»INU3  = ZC0MP 

IN3EX  = lSiJEX  + l 

- - IF  (iNOtX  .GT.IPTS)  GO  TO  3000  --  -- - 

200Q  CO.,Ti:jJ£ 


6YPASG  - UPDATE  THE  DACS 

LAUNCH  (1  ) =XCOPP  '‘.102375 

LAJNCH<2)  =YCO*iP  ♦.102375  

launch  (3  ) =ZCullP  ». 102375 

LA  JNCrt(4  ) ^RLB -*.102375 

LAwNCH(5)=COJ£  *.102375 

lAJ.NCH(6)=''P3  *.102375  

LAurtCH  (7  ) =KI  *.102375 

LAoNCH (8 ) =GAM *.102375  - 

LAjr:CH(9)  =EOOT  *.102375 

L ACr;CM  { 15  ) = THE  T At  *.102375 

LAUNCH(11)=RN  *.102375 


J 
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APPENDIX  G 

SUMMARY  OF  STINGER  DIGITAL  COMPUTER  FUNCTIONS 


DYNAMICS.  INC 


PRE-REAL  TIME  PART 


1.  Manual  Mode  Operation 

• Select  target  crossing  angle,  altitude,  and  range  from  stored 

table.  j 

• Select  target  maneuver  from  stored  table  ; 

• Select  target  maneuver  point  at  missile  launch  j 

• Adjust  maneuver  table  for  selected  point  at  launch  \ 

• Scale  maneuver  table  for  analog  compatibility 

• Determine  type  of  target  from  maneuver  table 

• Calculate  initial  target  elevation,  azimuth,  and  position 

« Calculate  euler  angles  for  earth  fixed- to-laimch  coordinate 
transformation 

• Transform  initial  target  velocity  into  launch  coordinate  system 

• Calculate  parameters  for  the  analog  variable  scale  factor 
equation 

• Generate  table  of  missile  range-to-go  relative  to  target  for  use  in 
controlling  the  A/D  data  collection  from  the  analog  computer 

• Scale  initial  launch  conditions  for  analog  compatibility 

• Transform  target  maneuver  table  into  launch  coordinates  and  scale 

• Clear  storage  area  where  data  is  to  be  accumulated 

2.  Automatic  Mode  Operation  Has  Same  Sequence  As  Manual  Mode 

Except;  j 

• Target  crossing  angle,  altitude,  range,  maneuver,  and  maneuver  j 

point  are  selected  from  a stored  simulation  plan  | 

• There  are  no  manual  time  delays  associated  with  selection  or 
with  problem  setup  and  execution  of  subsequent  simulations 
stored  in  the  simulation  plan 

G.  1 


f 


REAL  TIME  PART 


• Convert  initial  conditions  into  form  compatible  with  D/A 
converter 

• Wait  for  AD/ 4 analog  computer  to  signal  readiness  to  receive 
initial  conditions 

• Send  initial  conditions  D/A  to  the  AD/4  computer  and  start 
analog  real  time  simulation 

• Read  A/D  relative  position,  velocity,  and  time  from  analog 

• Convert  A/D  information  into  digital  computer  notation 

• Compare  A/D  relative  position  information  (range-to-go)  with 
previously  generated  table 

• A/D  information  is  collected  and  stored  as  data  whenever  the 
above  range-to-go  comparison  is  close 

• Real  time  interpolated  target  maneuver  is  stored  at  same  time 
as  A/D  data  collection  above 

• If  A/D  information  indicates  missile  has  passed  the  target,  the 
real  time  part  is  terminated  and  control  is  given  to  the  post- 
real  time  part 

• Calculate  target  angle  of  attack  in  real  time 

• Calculate  target  jet  plume  apparent  length- to- breadth  ratio  in 
real  time 

• Send  target  velocity  components  and  length- to-breadth  ratio  D/A 
to  the  analog  computer 

• Collection  of  last  desired  A/D  data  terminates  the  real  time  part 
and  gives  control  to  the  post  real  time  part 


G.2 


POST- REAL  TIME  PART 


I 

• Scale  collected  A/D  data  for  digital  compatibility 

• If  nonstationary  target,  calculate  Euler  angles  for  transformation 
from  launch-to- target  fixed  coordinates 

• Calculate  the  least  squares  missile  approach  trajectory  from 
collected  A/D  data  and  project  miss  distance  information 

• If  the  target  is  a high  speed  jet,  perform  a lethality  test  with  the 
probability-of-kill  model 

• Repeat  the  simulated  flight  20  times  and  obtain  average  PK,  mean, 
and  standard  deviation  of  miss  distance  information 

• Give  control  to  pre-real  time  part  for  next  simulation 


G.3 


o o n 


1 

UTIME  OACS  /‘*/74  OPT  = l FTNH  4. 24-75075  07 


SU'or'iCUUNE  ClACS  — 

niMc:NS10N  ADIN(IO) 

JAM, ON 

COM;-lO,'!/*aACl/17,  AOI.NTST  (10) 

J A M , C F r 

C0MM0N/EEEtE/ITl,KC<,lCR2, ID A3 , I NQ , I NDEX ,M AX , I WR I I E , NPAT H 

IF  {NFATH.NE.2)  TIM  = 0.0  — 

IF  (NPATH  .tQ.2)  TIM  = TIM4-.0  01 

GONPuT-  oIMPLIFItO  ADCS 


AOIN(1)=-100.0/1JO.O 

IF  (TIM.oF,4.G) AOIN (1) =100, 0*( TIM-5. U) /lOU.O  -• 

Al)IN(.')  = 0.0 

AuIN(J)=0.0  — 

C NOT::  TIME  HAS  AN  ADDED  FACTOR  OF  10.0 

A0IN(4)=TIM  /lU.O 
ADIN(5)=(11.0*riM-5.0)/100.0 

AOlN(b)= (50.D-10.Q*TIM)/10Q.fl  — 

AL1IN(7)=(5G.Q-10.0*TIM)/100.0 
Aj:n(3)=  (liJ.J*IIM-50.0)/100.U 
A JIN (9) =0.0 

ADIN{ia)-0.0  — 

DO  10  1=1,10 

10  ADINTSTC I)=ADIN{I) 

IF  (TIM. GT. 15.0)  CALL  SIMSTOP 

CALL  SIMIDLE 

END 


I 


i 


DISCLAIMER  NOTICE 


THIS  DOCUMENT  IS  BEST  QUALITY 
PRACTICABLE.  THE  COPY  FURNISHED 
TO  DDC  CONTAINED  A SIGNIFICANT 
NUMBER  OF  PAGES  WHICH  DO  NOT 
REPRODUCE  LEGIBLY. 
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